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ABSTRACT 


This  document  provides  the  guidelines,  limitations,  and  modifications 
required  to  perform  a  structural,  fracture  analysis  using  Grif f ith-Irwin 
fracture  mechanics  principles.  It  serves  as  an  introduction  to  fracture 
mechanics  for  those  personnel  who  are  concerned  with  fracture  strength  esti¬ 
mates  for  aerospace  structural  applications.  Illustrations  and  hypothetical 
examples  are  included  which  show  how  engineering  solutions  for  critical  crack 
size  and  fracture  stress  may  be  made.  The  critical  stress  intensity  (fracture 
toughness)  concept  is  used  as  a  basic  factor  for  the  fracture  analysis  of 
materials.  For  most  crack  situations,  a  stress  intensity  factor  can  be  com¬ 
puted  which  can  be  related  to  critical  conditions  and  estimates  made  of  criti¬ 
cal  crack  lengths,  stresses,  and  crack  propagation  behaviors.  To  provide  a 
complete  and  accurate  fracture  analysis,  the  user  is  encouraged  to  become 
familiar  with  all  aspects  of  the  analysis  and  its  limitations. 

This  document  has  been  approved  for  public  release  and  sale;  its 
distribution  is  unlimited. 
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Symbol 


SYMBOLS  AND  NORMAL  UNITS* 


Units 


o 

2a 

2a 


ec 

Aa,  da 

2a 

W 


A  «A„fTVr4 
str.  STR 

B 

b,  c 


C 

2c 

E 

E 


>EC 


str.  ’  STR 


F 

f 

G 

A 

A 


Crack  length  (usually  V2  total  length)  or 
crack  depth  (part- through  crack) 

Initial  crack  length  (usually  V2  total  length) 

Total  crack  length  (tip  to  tip) 

Critical  crack  length  (tip  to  tip) 

Critical  crack  length 

Equivalent  Griffith  crack  length 

Critical,  equivalent  Griffith  crack  length 

Incremental  crack  extension 

Crack  aspect  ratio,  (crack  length/panel  width) 
Area  or  constant 

Cross  sectional  area  of  stringer 
Material  thickness 

Length  from  center  line  to  point  of 
concentrated  force 

A  constant 

Crack  length  (part- through  crack) 

Young’s  Modulus  or  Modulus  of  plate  psi 

Modulus  of  stringer  psi 

(See  Figure  18) 

Force  per  unit  thickness  or  length 
A  function  of 
Shear  Modulus 

Strain  energy  release  rate  inch- 

Critical  for  instability  inch- 

Moment  of  inertia  of  stiffener  about  its 
neutral  axis  perpendicular  to  the  plane  of 
the  plate 


*  Except  where  noted  in  text. 
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inch 

inch 

inch 

inch 

inch 

inch 

inch 

inch 

inches^ 

inches^* 

inch 

inch 

inch 

3 

or  ksi  x  10 

3 

or  ksi  x  10 
inch 

lbs/ inch 

psi 

•lbs/ inch^ 

2 

*  lbs/inch 
inch^ 


Symbols 

Units 

Kr  ii*  hi 

Stress  intensity  factors  for  cracking  Modes 
of  Figure  5. 

ks  Winch* 

kt 

Ic 

Critical  stress  intensity  factor  for  Mode  I 
cracks  (Plane  Strain  Fracture  Toughness) 

ks  Winch' 

K 

c 

Critical  stress  intensity  factor  for  an 
infinite  sheet  for  Mode  I  cracks  (Plane 

Stress  Fracture  Toughness) 

ks  Winch1 

K 

c 

Average  critical  value  of  plane  stress, 

Mode  I  stress  intensity 

ksiN/inch ' 

KIscc 

Critical  stress  intensity  for  Mode  I  in 
environmental  (corrosive)  conditions 

ks  Winch' 

Kb 

Bulging  stress  intensity  factor 

ks  Winch 

KIi 

Initial  value  of  stress  intensity  factor 

ks  Winch* 

“h 

Hoop  membrane  stress  intensity  factor 

ks  Winch' 

Kth 

Threshold  value  of  stress  intensity  factor 

ks  Winch' 

kq 

Tentative  value  of  plane  strain  fracture 
toughness 

ks  Winch 

k 

Elastic  constant  for  plane  strain  or  plane 
stress 

M 
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I  INTRODUCTION 


This  report  contains  guidelines  by  which  fracture  mechanics  principles 
may  be  used  in  the  analysis  of  the  fracture  strength  of  aircraft  structures. 

It  serves  as  an  introduction  to  the  use  of  this  discipline  in  structural, 
fracture  analysis.  The  data  contained  herein  is  not  all  inclusive,  but  supple¬ 
mentary  to  other  information  available  in  the  field  of  fracture  mechanics. 

Fracture  mechanics  has  been  employed  to  evaluate  material  fracture  behav¬ 
ior  with  a  good  degree  of  success.  If  materials  can  be  ranked  by  their  criti¬ 
cal  stress  intensity  factor  (K^itical) ,  the  basic  unit  of  fracture  mechanics, 
then  this  unit  can  by  applied  to  the  analysis  of  aerospace  structure.  Such  a 
critical  parameter  has  been  determined  for  thick  section  fracture  and  provides 
the  proper  trends  in  determination  of  critical  fracture  stress  for  plane  strain 
behavior . 

Thus,  a  critical  K  can  be  used  in  fracture  analysis--provided  the  limita¬ 
tions  of  fracture  mechanics  are  fully  understood.  It  is  the  purpose  of  this 
report  to  show  how  fracture  mechanics  can  be  used  during  the  vehicle  design 
stages  for  fracture  stress  and  critical  crack  length  estimations  as  well  as 
optimum  materials  selection.  To  provide  a  basis  for  the  use  of  the  stress  in¬ 
tensity  factor  in  analysis,  it  is  necessary  to  review  its  historical  evolution, 
derivation,  terminology,  and  relationship  to  accepted  design  philosophies. 

This  is  reviewed  in  Sections  II  through  V.  In  many  structural  applications 
plasticity  is  not  limited  and  geometric  effects  predominate.  These  parameters 
must  be  accounted  for  and  lead  to  modified  solutions  to  the  critical  K.  concept 
(Sections  VI  and  VII).  Of  necessity,  other  influencing  factors  may  become 
dominant,  such  as  environmental  effects  and  require  additional  modification  of 
the  basic  critical  K  concept  (Section  XI). 

With  these  basic  concepts  as  background,  and  knowledge  of  the  limitations 
of  fracture  mechanics  which  require  modification  to  these  concepts,  it  is  then 
possible  to  compare  and  analyze  typical  structural  fracture  problems.  These 
solutions  are  presented  in  Sections  VIII  and  IX  and  are  obtained  from  the 
following  basic  premise  that  there  is  a  K  solution  for  the  crack  problem  of  in¬ 
terest  which  can  approach  critical  conditions.  Following  this  premise,  it  is 
then  possible  to  make  engineering  estimates  and  recommendations  based  on  frac¬ 
ture  mechanics  analysis  of  structure  which  may  become  cracked  by  any  means  dur¬ 
ing  its  development  or  service  lifetime.  The  accuracy  of  these  estimates  is 
dependent  on  the  choice  of  critical  stress  intensity  factor. 

It  is  the  purpose  of  this  report  to  make  the  user  of  fracture  mechanics 
analysis  aware  of  the  basic  concepts,  limitations,  and  usage  in  design.  This 
can  only  be  accomplished  in  a  step-by-step  review  of  all  sections  of  this 
document.  In  this  manner,  the  user  will  become  familiar  with,  and  proficient 
in,  the  use  of  fracture  mechanics  analysis. 
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II  HISTORY  OF  FRACTURE  MECHANICS 


To  understand  the  motivating  force  behind  the  rapid  rise  of  this  new 
discipline,  one  must  go  back  to  the  large  number  of  fracture  problems  which 
have  occurred  over  the  years.  In  1919,  a  2-million  gallon  hot  molasses  tank 
failed  at  a  rivet  line  with  complete  disaster.  Recently,  another  such  tank 
fractured,  resulting  in  an  estimated  clean-up  requiring  2-3  weeks. 

Fracture  mechanics,  from  its  unimpressive  beginning  with  an  English 
scientist  studying  the  failure  mechanisms  in  glass,  to  its  present  employment 
in  analyzing  rocket  motor,  tankage  structure,  and  aircraft  primary  structure, 
has  come  a  long  way.  During  the  war  and  post-war  era,  the  Liberty  Ship  inci¬ 
dents  which  saw  1450  welded  plate  failures,  19  of  which  were  total  ship  losses, 
were  a  prime  motivation  in  extending  the  theory  of  fracture  cf  brittle  solids 
to  more  realistic  structural  materials.  Thus,  the  theory  of  Griffith,  the 
English  scientist,  was  modified  and  updated  to  account  for  semi-ductile  behavior. 
There  have  been  other  fracture  problems,  of  course,  such  as  rashes  of  steam 
generator  turbine  and  motor  failures,  Polaris  motor  cases,  the  260-inch  solid 
rocket  engine  case,  to  name  a  few,  which  also  helped  to  establish  fracture 
mechanics  as  an  analysis  method,  as  well  as  a  material  evaluation  parameter. 

The  cost  of  replacement,  loss  of  hardware,  etc.,  have  dictated  the  necessity  of 
utilizing  the  tools  of  fracture  mechanics.  This  fact  is  even  more  important 
now  as  the  use  of  higher  strength  steels,  titanium  alloys,  and  exotic  materials 
become  more  evident  in  high  speed  structures. 

II. 1  THE  GRIFFITH/IRWIN  ERA 

The  first  in-depth  research  of  sharp-crack  fracture  was  performed  by 
A.  A.  Griffith(l)  of  the  Royal  Aircraft  Establishment,  England.  In  1920  he 
published  his  results  on  the  effect  of  surface  scratches  on  the  mechanical 
strength  of  solids.  His  work  was  truly  the  initial  step  in  the  foundation  of 
fracture  science.  Working  with  brittle  material  (glass),  Griffith  postulated 
that  an  existing  crack  or  flaw  will  propagate  if  the  total  system  energy  is 
lowered  in  the  presence  of  a  crack  or  flaw.  Therefore,  Griffith  provided  a 
means  of  estimating  the  theoretical  strength  of  solids  and  gave  the  correct 
relationship  between  fracture  strength  and  defect  size  for  brittle  material . 

Based  on  G.  R.  Irwin's  research  at  the  Naval  Research  Laboratory,  the 
Griffith  theory  of  fracture  was  ultimately  shown  to  be  not  only  strain  and 
surface  energy  dependent,  as  shown  by  Griffith,  but  also  highly  dependent  on 
the  work  of  plastic  deformat ion( 2 ) #  This  led  to  the  conclusion,  supported 
independently  by  E.  Orowan(^),  that  for  relatively  ductile  materials  the  work 
of  plastic  deformation  is  much  larger  than  that  of  surface  tension  for  engineer¬ 
ing  materials.  Subsequent  papers  by  Irwin  indicated  that  the  energy  approach 
to  fracture  is  equivalent  to  a  critical  stress  distribution.  From  this  he 
developed  the  concept  of  fracture  toughness . 

A  quantative  relationshop  could  now  be  placed  on  the  fracture  process 
which  included  the  dimensions  of  the  crack  or  flaw,  the  nominal  stress  field 
near  the  crack,  and  a  property  of  the  material  which  governed  the  energy 
balance  of  the  material  in  the  presence  of  a  crack  under  stress  or  the  so- 
called  fracture  toughness.  Thus,  for  the  first  time,  a  characterization  of 


2 


the  fracture  process  for  materials  of  limited  ductility  was  available  to  the 
stress  engineer  and  a  new  method  of  analysis  was  developed  called  fracture 
mechanics , 


The  Griffith/Irwin  fracture  mechanics  approach  involves  an  energy 
balance;  that  is,  the  energy  required  for  crack  growth  and  available  system 
energy  (stored  elastic  energy  released  by  the  material  when  the  crack  grows). 
When  the  energy  from  the  elastically  strained  material  surrounding  the  crack 
is  equal  to  or  exceeds  the  energy  required  to  support  crack  growth,  crack 
extension  will  occur  without  additional  increase  in  load  or  stored  elastic 
energy.  Griffith’s  research  in  brittle  glass^^  indicated  that  for  a  homo¬ 
geneous,  isotropic  material  containing  a  crack  of  length  2a,  having  a  thick¬ 
ness  B  undergoing  a  uniform  tensile  stress  perpendicular  to  the  crack,  the 
crack  would  extend  an  incremental  distance  (da)  on  either  end.  The  energy 
required  for  this  incremental  crack  extension,  which  created  new  fracture 
surfaces  2B  da,  was  the  energy  required  to  overcome  the  solid  state  surface 
tension  of  the  material,  Usurface.  (Due  to  symmetry,  only  one  half  of  the 
crack  is  being  analyzed.) 


dW 


req. 


U 


surface 


2B  da  * 


Thus,  the  surface  energy  required  per  incremental  crack  extension  is 


l 

B 


dW  dW 

■  recl:  =  -  re(l?  =  2U 

da  dA  surface 


(II-D 


where  or  ^  is  the  crack  surface  area.  Therefore,  Griffith  indicated  that 

the  surface  energy  is  a  constant  for  brittle  material. 


Griffith  then  showed  that  the  energy  available  (dWavai^#)  from  the  re¬ 
lease  of  stored  elastic  energy  as  the  crack  grew  was 


or 


dW  . .  =  (Stored  Elastic  Energy )x(Affected  Volume) 

avail.  ' 

dWavail.  IT  *  (  2  v  Ba  da] 


This  equation,  written  in  the  form  of  energy  available  per  incremental  crack 
extension,  is 


avail . 
B  da 


dW 

avail. 


dA 


7T0-2a 

E 


(H-2) 


where  E  is  the  elastic  modulus  of  the  material.  Griffith  then  formulated  the 
critical  energy  balance  equation  which  states 


avail.  surface 

or  from  Eqs.  II-l  and  II-2 
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(  dW  .,  \ 

avail.  I 

_  I  ( 

<  dW  \ 

req- 

B  ' 

i  da  / 

B  ' 

*  da  / 

*The  symbol  B  in  fracture  mechanics  is  standard  terminology  for  thickness. 
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and 


( II- 3 ) 


2U 


2 

TTCT 

surface  E 


or  written  in  the  familiar  Griffith  form,  equation  II-3  is 


cr 


critical 


2E  U 


surface 
7T  a 


(II-4) 


for  the  critical  energy  case. 

critical  *-s  t^ie  critical  stress  for  unstable  crack  propagation  (usually 
written  crc).  The  concept  stated  by  equation  II-4  is  that  the  critical 
stress  is  inversely  proportional  to  the  square  root  of  the  crack  length  and 
directly  proportional  to  the  square  root  of  the  surface  tension.  This  equa¬ 
tion  forms  the  basis  for  most  of  the  modern  theories  of  fracture  of  solids. 

Orowan  in  the  early  1950's  and  Irwin,  during  the  same  period,  showed  the 
influence  of  local  plastic  deformation  near  the  crack  tip  on  the  fracture 
control  process.  0rowan(3)  stated  that  for  local  plasticity,  the  work  of 
plastic  deformation  can  be  treated  in  a  manner  similar  to  Griffith's.  He 
also  postulated  that  if  the  region  of  crack  tip  plasticity  was  dimensionally 
small  compared  to  other  body  dimensions,  elastic-stress  analysis  could  be 
used  successfully.  Thus,  Orowan' s  modified  Griffith  equation  accounting  for 
plasticity  effects  was 


Orowan  also  estimated  that  the  work  or  energy  of  plastic  deformation  (Up^ast^c) 


was  1000  times  as  large  as  ^surface*  therefore,  the  latter  term  could  be  neglected 
in  the  analysis  of  materials  of  slight  ductility.  In  this  case,  equation  II-5 
becomes  _ 


2E  U 


plastic 

7ra 


cr 


c 


or 


U 


(II-6) 


This  analysis,  independently  supported  by  the  work  of  G.  R.  Irwin,  was 
the  first  recognition  of  the  role  of  plasticity  in  fracture. 


It  became  apparent  with  additional  research  by  Irwin(2)  and  his  asso¬ 


ciates,  that  the  energy  of  plastic  deformation  (Up^ast£c)  was  not,  however, 
a  constant,  but  a  variable  dependent  on  length  of  crack;  i.e.,  the  larger  the 
crack  the  greater  the  plastic  energy.  To  illustrate  this  finding,  refer  to 


Figure  1.  A  crack  of  original  half  length  aQ  under  increasing  load  will 


grow  slowly  and  shows  an  increase  in  plastic  energy  dissipation,  2Upiastic 
(dashed  curve).  The  shape  of  this  curve  will  be  dependent  on  specimen 
geometry  (width,  thickness,  etc.)  and  crack  speed.  During  the  rising  portion 
of  the  plastic  energy  dissipation  curve  (crack  growth  resistance,  or  R  curve 
is  the  now  accepted  terminology),  resistance  to  crack  growth  occurs  with 
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INCREASING  ENGERGY  INPUT 


INCREASING  HALF  CRACK  LENGTH 


FIGURE  I.  GRAPHIC  REPRESENTATION  OF  THE  FRACTURE  PROCESS  FOR  IDEALLY 
BRITTLE  AND  OTHER  THAN  BRITTLE  BEHAVIOR 
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increased  plastic  deformation.  Additional  input  energy  is  required  (increase 
in  external  load)  until  the  point  of  tangency  of  the  R  curve  and  the  avail¬ 
able  energy  curve  is  reached  (see  point  A,  Figure  1  ).  At  this  point,  un¬ 
stable  crack  growth  occurs,  i.e.,  the  energy  available  for  crack  extension  is 
greater  than  the  dissipated  plastic  energy.  Therefore,  the  crack  growth 
process  is  now  unlimited  with  no  further  increase  in  load  (input  energy)  for 
crack  extension,  and  the  crack  grows  catastrophically.  One  can  visualize 
this  process  by  mentally  walking  up  the  R  curve.  Extra  energy  is  needed  to 
get  up  the  hill,  but  once  the  crest  is  reached  (point  A)  momentum  will  pre¬ 
vail  and  less  effort  is  required  to  proceed. 


Point  D  and  its  associated  resistance  curve  represents  typical  behavior  of 
an  ideal,  brittle  material.  It  will  be  noted  that  negligible  crack  extension 
occurs  prior  to  instability. 


As  shown  by  Irwin,  the  energy  available  per  crack  extension  is  really 
the  elastic  strain-energy  release  rate  (j£?)  associated  with  increasing  load. 
Thus,  at  the  point  of  tangency,  (A,  Figure  1  )  the  following  conditions 

prevail : 


R  = 


critical 


and 

dsi?  _  dR 
dR  da 


From  these  relationships,  it  follows 


or  from  Eq.  II-6 


critical 


2U 


plastic 


(II-7) 


A  quantative  relationship  between  critical  half  crack  size  (ac),  stress 
(°"c)>  and  fracture  toughness  is  shown  by  equation  II-7.  The  fracture 
toughness  is  represented  by  the  critical  material  strain-energy  release  rate 
(*/c).  This  equation  indicates  material  fracture  behavior  on  a  gross  scale. 

The  control  of  fracture  toughness  will  not  be  a  constant,  but  dependent  on 
crack  length,  (a).  To  date,  an  analytical  means  is  not  available  to  deter¬ 
mine  the  plastic  work  associated  with  the  Griffith/Irwin  instability  equation. 
However,  the  problem  was  solved  experimentally  by  Irwin  and  his  co-workers. 

The  parameters  needed  to  determine  a  critical  strain-energy  release  rate  in 
equation  II-7  are  the  load  and  crack  length  associated  with  critical  crack 
extension.  The  load,  and  hence,  critical  stress  can  be  determined  directly 
from  any  universal  testing  machine.  The  question  remained: how  could  the  assoc¬ 
iated  critical  crack  length  (aQ  or  2ac)  be  determined  accurately?  This  problem 
was  solved  by  Irwin^using  an  ink  stain  technique.  Using  a  centrally  slotted 
panel,  he  placed  a  small  quantity  of  staining  ink  at  the  notch  roots.  The 
specimen  was  then  uniaxially  loaded  perpendicular  to  the  slot.  During  the 
rising  portion  of  the  R  curve  (see  Figure  1),  the  ink  could  slowly  follow  the 
crack  extension  by  capillary  action.  Once  the  critical  point  of  crack  exten¬ 
sion  was  reached,  the  ink  could  not  keep  up  with  the  rapidly  advancing  crack; 
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therefore,  an  indication  of  critical  crack  length  would  be  outlined  by  the 
distance  the  ink  traveled  during  slow  crack  extension.  However,  consistent 
calculations  oiJPc  could  not  be  made.  The  ink  stain  method,  due  to  interpre¬ 
tation  and  experimental  difficulties,  has  since  been  replaced  by  more  sophis¬ 
ticated  techniques,  which  will  be  discussed  later  in  this  report. 

Thus,  for  the  past  50  years,  fracture  research  has  shown  the  evolution 
of  "a  new  discipline  termed  f fracture  mechanics*,  which  provides  a  rational 
basis  for  analysis  and  control  of  flaw  related  fracturing  of  structures. "(6) 

II. 2  ASTM  TASK  FORCE 

In  1958,  the  Department  of  Defense  initiated  a  study  effort  within  the 
American  Society  for  Testing  and  Materials  (ASTM)  motivated  by  several  frac¬ 
ture  problems  in  rocket  propellent  tanks.  As  the  result  of  this  contact,  a 
special  ASTM  committee  was  established  on  Fracture  Testing  of  High-Strength 
Sheet  Materials.  (-*)  The  importance  of  the  structures  fracture  problem  can 
be  seen  by  the  appointment  of  6  representatives  from  the  aerospace  industry 
and  government  to  this  special  10-man  committee.  Much  will  be  said  later  in 
this  report  about  the  contributions  of  this  committee  to  the  standardization 
of  fracture  testing  methods,  terminology,  and  analysis. 

The  relationship  of  the  analysis  by  fracture  mechanics  of  aircraft  and 
aerospace  structures  in  preventing  fracture  in  areas  which  must  be  fracture 
or  damage  tolerant  will  become  evident  during  the  development  of  this  report. 


II. 3  SPECIMENS 


One  condition  which  must  prevail  in  order  to  determine  a  vqlid  fracture 
toughness  value  is  that  limited  plasticity  accompany  the  fracture  process.  It 
will  be  recalled  that  the  solution  of  Griffith,  Orowan,  and  Irwin  are  baaed  on 
linear  elastic  approximations;  therefore,  the  crack  tip  plasticity  must  be 
confined  to  a  small  zone  compared  to  other  specimen  dimensions.  For  this 
reason,  many  different  specimen  designs  have  been  introduced  and  evaluated 
since  the  first  fracture  toughness  specimen.  As  more  specimen  types  were 
introduced  (notched  rounds,  edge  notch,  etc.),  it  became  obvious  that  stan¬ 
dardization  of  test  specimens  and  methods  was  required  to  provide  valid  frac¬ 
ture  toughness  measurements.  The  efforts  of  experimentalists  and  analysts 
within  The  American  Society  for  Testing  and  Materials  have  been  instrumental 
in  evaluating  these  design  and  test  methods.  The  primary  contributors  to  this 
effort  have  been  associated  with  the  Naval  Research  Laboratory  (G.  R.  Irwin, 

J.  M,  Krafft,  A.  M.  Sullivan,  J.  A.  Kies,  R.  W.  Boyle,  and  their  co-workers) 
and  NASA  Lewis  Research  Center  (J.  E.  Srawley  and  W.  F.  Brown,  Jr.,  and  their 
co-workers).  (See,  for  example,  Ref.  7  and  8.)  Illustration  of  the  various 
test  specimens  in  use  today  will  be  shown  later  in  Section  XII  of  this  report. 
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Ill  DESIGN  PHILOSOPHIES 


Fatigue  or  structural  damage  can  be  expected  during  the  operational  life 
of  aircraft  or  spacecraft  structures.  To  account  for  this  in  modern  aerospace 
vehicle  structural  design,  two  philosophies  are  frequently  considered  to  ensure 
attainment  of  a  satisfactory  operational  life.  One  of  these  is  the  "fail  safe" 
philosophy,  which  is  also  associated  with  the  damage  tolerant  concept.  The 
other  is  the  "safe  life"  concept,  which  is  really  a  given  life  or  specified 
lifetime  philosophy.  In  the  following,  these  two  philosophies  will  be  further 
defined  and  their  applications  explored. 


III.l  FAIL  SAFE 


The  fail  safe  concept  has  been  applied  to  everything  from  electronic 
hardware  to  load-carrying  structure  in  aircraft  systems.  In  structural  design, 
its  philosophy  states  that  during  the  service  life  of  a  structure,  fatigue 
cracks  or  damage  will  not  progress  to  a  catastrophic  condition  prior  to  detec¬ 
tion  during  regular  inspection  periods.  Therefore,  the  fail  safe  philosophy 
assumes  that  readily  detectable  damage  due  to  fatigue  or  accidental  occurrence 
will  remain  in  a  somewhat  stable  condition  upon  subsequent  service  load  his¬ 
tory.  In  aircraft,  this  concept  covers  all  of  the  primary  structure,  so  that 
loss  of  one  conponent  does  not  reduce  the  total  strength  of  an  assembly  to 
dangerous  levels.  Most  of  these  components  consist  of  two  or  more  separate 
segments  constructed  in  a  manner  that  a  crack  can  propagate  completely  through 
one  section  without  propagating  into  a  second  area;  or  perhaps  a  structure  is 
designed  or  materials  selected  to  take  advantage  of  slower  crack  growth  rates, 
thus  effectively  increasing  the  probability  of  crack  detection  during  inspec¬ 
tion. 


FAA  certification  for  transport  aircraft  neccessitates  the  fail  safe  con¬ 
cept  and  dictates  that  particular  structural  areas  be  analyzed  and  tested  to 
meet  fail  safe  requirements.  In  commercial  design,  concern  is  for  the  unde¬ 
tected  crack,  or  puncture  removing  a  structural  panel  between  stringers  or 
frames.  It  also  applies  to  the  fatigue  or  damage  of  frames  and  stringers. 
However,  this  concept  cannot  be  overlooked  in  military  aircraft,  either.  In 
military  aircraft,  similar  damage  can  be  caused  by  projectiles  or  shrapnel, 
and  fail  safe  concepts  apply. 

The  fail  safe  or  damage  tolerant  philosophy  often  utilizes  stringers  as 
built-in  crack  or  tear  stoppers.  Tear  straps  are  now  an  integral  part  of  al¬ 
most  all  new  commercial  designs.  The  function  of  these  straps  as  crack  stop¬ 
pers  will  be  discussed  further  in  other  sections  of  this  report.  Using  exist¬ 
ing  fracture  mechanics  principles,  a  fail  safe  analysis  can  be  made  in  which 
critical  stresses  can  be  determined  for  planned  designs.  Because  the  ultimate 
substantiation  of  these  assemblies  requires  proof  testing,  a  good  comparison 
can  be  made  with  the  fracture  analysis.  However,  no  indication  is  usually  pro¬ 
vided  as  to  structural  fatigue  resistance.  Improvements  in  design  are  usually 
made  by  varying  testing  parameters  and  relying  on  fracture  mechanics  principles 
for  estimating  fail  safe  stresses. 
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III. 2  SAFE  LIFE 


The  safe  life  concept  assumes  that  a  given  structure  will  not  develop 
fatigue  cracks  during  its  service  life.  Once  cracks  do  occur  in  the  safe  life 
structure,  its  life  is  assumed  complete* 

It  is  the  aim  of  this  concept  that  during  the  design  phase,  materials, 
stress  raisers,  and  other  fatigue-effecting  parameters  be  thoroughly  considered 
to  prevent  the  development  of  catastrophic  cracks  under  service  loading  condi¬ 
tions.  To  achieve  this  goal,  fatigue-sensitive  structure  is  subjected  to  simu¬ 
lated  service  testing  during  the  aircraft  development  stage.  Thus,  a  conserva¬ 
tive  or  realistic  test  program  provides  an  estimate  of  service  life  and  esta¬ 
blishes  inspection  and  maintenance  schedules. 

It  can  be  seen  that  in  many  cases  this  concept  can  lead  to  costly  retro¬ 
fitting  when  thoroughly  analyzed  safe  life  structure  develops  fatigue  damage 
during  a  testing  phase.  It  is  quite  possible  that  this  may  take  place  some 
time  after  production  items  have  been  delivered.  Problems  also  arise  when 
service  loads  are  not  accurately  simulated,  or  a  particular  design  is  pressed 
into  a  different  service  load  history  than  that  for  which  it  had  been  safe 
life  designed. 

Thus,  the  principles  of  safe  life  design  depends  on  a  fairly  accurate 
estimate  of  service  loading,  usually  involving  a  random  spectrum.  All  designs 
must  then  be  analyzed  for  fatigue,  keeping  in  mind  the  starting  points  for 
damage  such  as  stress  raisers  (welds,  cutouts,  rivet  holes,  notches,  etc.), 
and  the  material  selection  should  be  optimized  for  suppression  of  fatigue. 

It  can  now  be  seen  that  both  philosophies  of  structural  design  have  cer¬ 
tain  limitations.  However,  both  concepts  can  and  should  be  used  in  aircraft 
design  to  provide  the  optimum  f,safeff  structure  with  fracture  mechanics  as  an 
analysis  tool  and  proof  testing  as  a  back-up  requirement  for  critical  structure. 
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IV  PLANE  STRAIN  VS.  PLANE  STRESS  IN  FRACTURE  ANALYSIS 


IV. I  "ELASTIC  OR  BRITTLE"  FRACTURE 


The  terms  brittle  or  elastic  fracture,  in  general,  refer  to  an  elasti¬ 
cally  stressed  body,  excluding  that  region  surrounding  the  crack  as  the  crack 
grows  during  the  tearing  process.  Through  the  terms  plane  strain  and  stress, 
we  can  tie  together  the  fracture  process  with  fracture  mechanics. 


IV. 2  PLANE  STRAIN 


In  a  cracked  body,  a  zone  of  plasticity  (see  Section  VI)  occurs  at  the 
crack  tip  during  increase  of  a  remotely  applied  force.  In  a  thick  body,  this 
plastic  action  is  suppressed  by  interior  constraint  due  to  thickness.  This 
is  analogous  to  the  definition  of  plane  strain  as  given  in  Reference  9;  that 
is,  if  the  z  dimension  of  a  stressed  body  (thickness)  is  large,  then  all  z 
direction  strains  and  displacements  are  zero.  This  definition  could  well 
describe  the  behavior  of  a  thick  section  fracture,  except  for  that  region 
near  where  the  crack  meets  the  free  surface  faces.  In  that  region  the  z 
direction  strains  and  displacements  are  not  zero,  and  lead  to  the  development 
of  shear  deformation  or  "shear  lips"  near  the  free  surfaces.  Therefore,  in 
a  thick  section,  plane  strain  conditions  are  governed  by  the  material  thick¬ 
ness. 


One  may  visualize  an  analogous  situation  with  a  common  tensile  test. 
Shown  in  Figure  2a  is  a  representative  fracture  for  a  material  of  limited 
ductility,  elastic,  except  for  the  region  surrounding  the  fracture.  A  negli¬ 
gible  amount  of  local  plasticity  occurs,  though  the  thickness  and  fracture 
occurs  as  indicated  by  the  load  deformation  curve. 

IV. 3  PLANE  STRESS 


In  a  thin,  cracked  body  subjected  to  in-plane  loads,  there  is  essentially 
no  constraint  in  the  z  (thickness)  direction  (3-dimensional  strain  state),  and 
plane  stress  conditions  prevail.  This  lack  of  constraint  is  contrasted  with 
the  plane  strain  condition  and  leads  to  relatively  large  plane  stress  yield 
zones  which  are  on  the  order  of  the  material  thickness  itself. 

Using  the  tensile  test  analogy  once  again,  Figure  2b  indicates  a  fracture 
in  a  ductile  material  which  is  elastic,  except  for  the  fracture  region.  Gross 
cross  section  deformation  has  taken  place  and  typical  shear  failure  predomi¬ 
nates.  At  this  point  it  must  be  realized  that  the  tensile  test  analogy  of 
Figure  2  is  for  illustrative  purposes  only  and  does  not  depict  plane  strain  or 
stress  (fracture  mechanics)  behavior. 

The  material  thickness  enters  into  the  modes  of  failure  and  is  the  con¬ 
trolling  factor.  Thus,  for  a  given  material,  plane  strain,  plane  stress,  or 
mixtures  of  the  two  can  be  encountered  by  changing  the  material  thickness,  and, 
therefore,  the  amount  of  plastic  constraint.  Examples  of  the  usual  fracture 
appearance  of  a  material  with  plane  strain,  plane  stress,  and  mixed  mode  type 
failures  are  shown  in  Figure  3. 


*The  general  cross-sectional  shape  is  material  and  thickness  dependent. 
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FIGURE  2a.  PLANE  STRAIN  BEHAVIOR  (PICTORIAL) 
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FIGURE  2b.  PLANE  STRESS  BEHAVIOR  (PICTORIAL) 
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PLANE  STRESS  (SLANT)  FRACTURE 


PLANE  STRAIN  (FLAT)  AND  MIXED  MODE  FRACTURE 


FIGURE  3.  THROUGH-THE -THICKNESS  VIEWS  OF  PLANE  STRESS 
&  PLANE  STRAIN  &  MIXED  MODE  FRACTURES 


PLANE  STRESS  (V-SLANT)  FRACTURE 
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A  familiar  residual  strength  diagram  is  shown  in  Figure  4.  These  frac¬ 
ture  envelopes  would  be  representative  of  an  aluminum  alloy  containing  central 
cracks  which  were  small  compared  to  specimen  width.  For  long  cracks,  plane 
strain  behavior  predominates  where  the  nominal  stress  is  small  compared  to 
the  material  yield,  and  thickness  is  sufficient  to  suppress  gross  crack  tip 
plastic  action.  In  this  region,  the  fracture  criterion  follows  that  proposed 
by  Griffith  (see  Eq.  II-4).  As  the  crack  length  or  thickness  decreases,  a 
region  of  mixed  mode  failure  (plane  stress  and  plane  strain)  predominates  due 
to  increased  plastic  deformation  until  a  condition  of  plane  stress  predomi¬ 
nates.  The  elastic  based  theory  of  fracture  must  then  be  modified  to  account 
for  increased  plastic  action. 

IV. 4  FRACTURE  APPEARANCE 


A  post  mortem  investigation  of  a  fracture  face  can  tell  much  about  the 
nature  of  the  stress  state,  as  well  as  identify  the  possible  failure  mechan¬ 
isms.  In  fatigue,  the  appearance  of  the  fracture  surface  can  provide  clues 
to  the  investigator  as  to  the  prior  stress  history.  In  fracture,  a  macro¬ 
scopic  view  of  the  fracture  surface  can  supply  information  as  to  the  prior 
state  of  stress  also.  In  this  way,  a  decision  can  then  be  made  as  to  the 
type  of  fracture  analysis  which  is  appl icable--plane  stress  or  plane  strain. 

If  a  through- the- thickness  crack  has  developed,  the  preference  for  crack 
progression  during  fracture  can  be  associated  with  the  plane  stress,  plane 
strain,  or  mixed  mode  regions  of  Figure  4,  as  shown  with  representative 
fracture  surfaces  of  each  region.  Therefore,  fractures  under  plane  stress 
conditions  usually  lead  to  through- the- thickness ,  45-degree,  slant  or  V-slant 
cross  sections.  The  development  of  the  slant  fractures  at  the  fracture  faces 
are  referred  to  as  Mshear  lips,f  and  provide  one  of  the  many  indications  as 
to  the  validity  of  plane  strain  fracture  toughness  tests.  The  significance 
of  the  shear  lip  development  will  be  discussed  in  Section  VI  of  this  report 
in  more  detail. 

The  square  appearance  of  flat  fracture  indicative  of  plane  strain  frac¬ 
ture  is  also  shown  in  Figure  4.  In  this  case,  plastic  deformation  (through- 
the- thickness )  at  the  crack  tip  is  minimized  and  shear  lip  development  negli¬ 
gible.  It  is  within  this  region  that  fracture  mechanics  works  with  a  high 
degree  of  confidence,  because  the  equations  of  fracture  are  based  on  elastic 
analysis . 

The  plane  stress  and  mixed  mode  regions,  although  under  intensive  study, 
require  modification  to  the  elastic-analysis  based  theory  in  order  to  use 
fracture  mechanics  effectively.  At  the  present  time,  it  is  these  modifications, 
many  of  which  are  empirically  based,  which  produce  the  greatest  limitations 
on  the  use  of  fracture  mechanics  in  design  of  aircraft  structures.  However, 
the  basic  concepts  are  sound,  and,  through  an  understanding  of  the  limitations 
of  fracture  mechanics,  a  working  knowledge  of  the  extension  of  the  analysis  to 
more  ductile  behavior  can  be  used  with  confidence. 

With  this  background  in  theory,  history,  and  basic  failure  modes  of  frac¬ 
ture,  we  can  now  approach  the  analysis  of  the  stress  state  in  the  presence  of 
cracks  and  then  a  major  concept  of  fracture  mechanics-- the  stress  intensity 
factor. 


13 


■°  * SS3HXS  TVNIHON 


03 

Csl 

£ 

o 

2 

W 

o 

2 

o 


14 


FIGURE  4.  RESIDUAL  STRENGTH  DIAGRAM  SHOWING  REGIONS  OF  FRACTURE  BEHAVIOR  AND  GENERAL  TRENDS  IN 
FRACTURE  APPEARANCE  FOR  PLANE  STRAIN,  PLANE  STRESS,  AND  MIXED  MODES 


V  LINEAR  ELASTIC  STRESS  ANALYSIS  OF  CRACKED  PLATES 


Through  the  use  of  linear  elasticity,  information  can  be  gained  as  to 
the  effect  of  cracks  on  the  stress  distribution  of  an  essentially  elastic 
solid.  Griffith(l)  made  use  of  the  Inglis  stress  solution(lO)  for  an  ellip¬ 
tical  hole  in  a  uniform  tension  field.  He  then  degenerated  the  elliptical 
hole  into  a  crack  solution  on  which  his  energy  rate  theory  is  based. 
Sneddon(H)  in  ^946  gave  the  first  solution  for  stress  field  expansion  near 
the  crack  tip.  However,  it  was  not  until  ten  years  later  that  Irwin^l^)  and 
M.  L.  Williams^*^  observed  the  general  applicability  of  these  solutions  and 
provided  the  extensions  which  made  it  possible  to  apply  them  to  a  cracked 
isotropic  elastic  body. 


V.l  MODES  OF  CRACK  SURFACE  DISPLACEMENT 


In  fracture  mechanics,  three  displacement  modes  of  the  crack  surfaces  are 
considered.  These  modes  are  shown  in  Figure  5.  Indicated  are  the  opening 
(I),  sliding  (II),  and  tearing  modes  (III).  The  sliding  mode  may  also  be 
referred  to  as  the  edge-sliding  mode.  It  can  be  seen  from  Figure  5  that  any 
crack  loading  problem  could  be  solved  by  including  combinations  of  two  or 
three  modes  for  various  loading  situations. 

V . 2  STRESS  INTENSITY  FACTOR 


Since  the  crack  surfaces  are  stress  free  boundaries,  that  portion  near 
the  tip  of  the  crack  predominantly  influences  the  local  stress  field.  There¬ 
fore,  remote  forces  and  boundaries  only  affect  the  intensity  of  this  local 
stress  field.  This  intensity  is  defined  as  the  crack  tip  stress  intensity 
factor,  Kj,  which  will  be  developed  in  succeeding  paragraphs.* 

Consider  the  opening  mode  (I)  crack  displacement  of  Figure  5.  Near  the 
leading  edge  of  the  crack,  the  free  body  diagram  would  appear  as  shown  in 
Figure  6.  The  stress  and  displacement  field  can  be  analyzed  as  a  plane  strain 
(extensional )  problem  in  elasticity.  Irwin  using  the  method  of  Wester- 

gaard(l^),  analyzed  this  problem  and  indicated  the  following  stress  fields 
and  displacements  for  Mode  I  at  a  distance,  r,  close  to  the  crack  tip. 


Stress  Equations  (Mode  I) 
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*A  physical  interpretation  of  the  stress  intensity  factor 
reflects  the  redistribution  of  stress  in  an  elastic  body 
tion  of  a  crack  and  which  reflects  mode  and  magnitude  of 
through  the  tip  region. 
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OPENING  MODE  I 


FIGURE  5.  MODES  OF  CRACK  SURFACE  DISPLACEMENT 
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FIGURE  6.  COORDINATES  AND  STRESS  COMPONENTS  OF  CRACK  TIP  FIELD 
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cr  =  0  for  plane  stress  conditions 
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or  cr  =  v(  cr  +  cr  )  for  plane  strain  conditions 
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Displacement  Equations  (Mode  I) 


and  w  =  0  for  plane  strain  conditions. 

V. 3  THE  WESTERGAARD  METHOD  (STRESS) 


(V-2) 


One  of  the  several  methods  of  analyzing  cracked  plates  to  obtain  stress 
intensity  solutions  will  be  presented  here.  As  mentioned  previously,  there 
are  several  methods  available,  and  a  good  review  is  contained  in  Reference  15. 
The  Westergaard  method  will  be  summarized  for  Mode  I  cracks  of  geometry  simi¬ 
lar  to  Figure  7.  Using  the  Westergaard  method,  K  can  be  computed  from  the 
(Westergaard)  stress  function  Z.  Obtaining  the  stress  function  involves 
guessing  for  any  particular  configuration;  however,  for  those  well  versed  in 
elasticity,  it  is  not  difficult.  For  the  Griffith  crack  configuration  of 
Figure  7,  with  the  coordinates  of  Figure  6  and  elastic  material,  the  plane 
extension,  elastic  equilibrium  equations  are  given  as: 
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(V-3) 


and  the  strain/displacement  relationships  and  Hookean  conditions  lead  to  the 
compatibility  equation, 
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FIGURE  7  GEOMETRY  OF  CENTRAL  CRACKED* INFINITE  PLATE 
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(V-4) 
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By  defining  an  Airy  stress  function,  *  ,  the  equilibrium  equations  V-3 
are  satisfied  in  terms  of  stresses 
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(V-5) 


Substitute  equations  V-5  into  V-4 


v4*  -  v2  (v2*)  -  0 


(V-6) 


To  solve  the  crack  problem,  the  stress  function,  ♦  ,  must  satisfy 
equation  V-6  and  its  boundary  conditions.  The  stress  function,  ♦  ,  is  chosen 
as, 
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(V-7) 
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which  will  satisfy  equation  V-6  if  ^!s  are  each  harmonic,  i.e«,  V  “i  *  0, 

Define  a  complex  variable  z  as  z  *  x  +  iy  and  then  functions  of  this 
variable  will  be  Z(z),  and  its  derivatives  are 
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(V-8) 


which  have  harmonic  real  and  imaginary  parts  if  the  function  is  analytic, 
i.e.,  Z  =  Real  Z  +  i  Imaginary  Z.  Hence, 


V2(Re  Z)  -  V2(Im  Z)  -  0 


(V-9  ) 


20 


This  solution  is  the  result  of  Cauchy -Riemann  conditions  which  specify, 
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(V-10) 
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By  using  equations  V-10,  the  functions  "Z  through  Z  can  be  differentiated. 


V. 3(a)  Mode  I 

(14) 

It  is  convenient  to  use  a  specific  form  of  the  Westergaard  stress 

function  for  Mode  I  cracks.  Westergaard  defined  an  Airy  stress  function,  <f> 

<f>  =  Re  1I  +  y  Im  J  (V-ll) 

where  Zj  and  Zj  are  related  through  equations  V-8.  (This  is  true  for  the 
Mode  I  crack  except  at  singular  points,  i.e.,  at  concentrated  load  points.) 
Using  equation  V-10  to  differentiate  equation  V-ll  to  form  the  stresses 
according  to  equations  V-5,  we  have, 


a  =  Re  Z,  —  y  Im  ZT 1 

x  I  J  I 

a  -  Re  ZT  +  y  Im  ZT 1 

y  I  J  I 

T  =-  y  Re  Z  1 
xy  J  I 


(V-12) 


At  this  point  in  the  analysis,  any  function  Z  (z),  which  is  analytic  will 
give  stresses,  by  equations  V-12,  which  automatically  satisfy  the  equations 
of  elasticity.  It  only  remains  to  find  those  functions  of  Z  (z)  which  will 
solve  the  crack  problems,  and  satisfy  boundary  conditions. 


Near  the  crack  tip  the  stress  free  crack  surfaces  are  the  boundaries 
which  dictate  Z  (z).  With  the  coordinates  taken  at  the  right  hand  end  of  a 
crack  parallel  to  the  x-axis,  Z  has  the  form 


£  (z) 

sf* 


(V-13) 


where  j(z)  is  a  well-behaved  function  and  approaches  a  real  constant  at  the 
origin.  This  results  in  T  Xy  and  Oy  approaching  zero  at  the  crack  surfaces 
because  the  surfaces  are  stress  free.  The  character  of  $  (z)  away  from  the 
crack  tip  is  then  unspecified  and  can  be  adjusted  to  solve  many  Mode  I  crack 
configurations  (symmetrical  extension) . 
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Close  to  the  crack  tip  region,  or  when  |z|  - *"0,  f  (z)  may  be  replaced 

by  a  real  constant  Kj,  and  equation  V-13  becomes 


K, 


(V- 14 ) 


|  z| — *0  sflTTz 


The  constant  Kj  may  be  thought  of  as  representing  the  constant  term  in 
the  series  of  the  MacLaurin  expansion  of  f  (z).  By  taking  polar  coordinates 
at  the  origin  (see  Figure  6) 


z  =  r  e 


ie 


( V- 15 ) 


and  using  equations  V-14  and  V-15,  the  stresses  can  be  computed  according  to 
equations  V-12  as 


K, 


v/  2  7T 
Kt 


xy 


sJTW 

h 

sfl  V 


9 

COS  — 

2  L 


9 

cos  — 
2 


9 


3  9 


sin  —  sin  — 

2  2  J 


9 


3  9 

1  +  sin  —  sin  — 
2  2 


(V- 16 ) 


9  9  36 

sin  —  cos  —  cos  — 

2  2  2 


which  are  the  "exact”  elastic  crack  tip  stress  field  equations  as  r — -0,  and 
in  that  general  vicinity  (see  also  equations  V-l).  Equations  V-16  neglect 
only  the  higher  order  terms  in  r  (beyond  the  constant  terms  in  MacLaurin’ s 
expansion  of  f  (z)  ).  The  factor  Kj  is  called  the  crack  tip  stress  intensity 
factor  and  has  a  constant  value  in  the  vicinity  of  the  crack  tip  through  the 
unspecified  character  of  £  (z) .  Kj  is  seen  to  depend  on  the  mode  of  loading. 

V.4  THE  WESTERGAARD  METHOD  (DISPLACEMENT) 

The  y  direction  strain,  by  Hooke’s  law,  is 

a 


e  =  - -  —  ( cr  +  <r  ) 

y  d  y  E  e  v  x  7.’ 

For  plane  strain  ( €  =  0)  conditions,  Hooke's  law  states 

cr  =  v  (cr  +  c r  ) 
z  x  y 


( V— 1 7 ) 


(V-18) 


Substituting  equations  V-12  and  V-18  into  equation  V-17,  and  integrating, 


gives 


=  L±ii  ^2(1  -  V  )  Im  Zj-yRe  ZjJ 


(V-19 ) 
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and 


u 


1  +  V 
E 


(V-20) 


[a- 


IV)  Re 


ZI~ 


y  Im  z. 


Noting  that  E  =  2G(1  +  V  )  and  substituting  equation  V-14  in  its  polar  form 
into  the  general  displacement  equations  V-19  and  V-20,  u  and  v  become, 


u 


V 


—  v/r/2  IT  cos  — 
G  2 


1 


0  fJ  ,  .  2  6 
2  V  +  sin  — 

2 


—  >/r/ 2  7T  sin  — 
G  2 


9  , ,  2  0 

2  1/  —  cos  — 

2 


( V- 2 1 ) 


and  w  =  0  for  plane  strain  conditions.  Equations  V-16  and  V-21  are  the  same 
as  the  stress  and  displacement  solutions  for  Mode  I  through  cracks  developed 
by  Irwin^^)  using  the  Westergaard  (stress)  method  (see  equations  V-l  and 
V-2.  Similar  solutions  can  be  obtained  for  Modes  II  and  III,  and  are  given 
in  Ref.  15. 

It  must  be  remembered  that  this  is  by  no  means  the  only  method  to  obtain 
stress  solutions  for  crack  problems.  However,  it  presents  one  of  the  simplest 
means  to  analyze  cracks. 

This  two  dimensional  analysis  of  Mode  I  cracks  indicates  that  the  magni¬ 
tude  or  intensity  of  this  distribution  is  dependent  on  Kj.  Thus,  the  rela¬ 
tive  intensity  of  the  stress  field  at  any  point,  obtained  by  basic  stress 
analysis  of  the  crack  tip  region,  is  the  crack  tip  stress  intensity  factor, 

K.* 

V. 5  DIMENSIONAL  CONSIDERATIONS 


It  now  becomes  necessary  to  discuss  the  influence  of  crack  size,  geometry, 
or  applied  load  on  the  overall  fracture  problem.  Consider  the  geometry  and 
stress  conditions  of  Figure  7.  A  crack  of  half  length  "a"  under  a  uniform 
gross  tensile  stress,  O'  ,  is  introduced  into  an  infinite  plate.  Mode  I  dis¬ 
placement  is  evident,  and  two  parameters  are  known,  cr  ,  and  "a”.  Due  to 
crack  symmetry,  only  one  crack  tip  stress  intensity  need  by  considered.  Using 
dimensional  analysis,  P.  C.  Paris'^)  solved  this  problem  through  the  use  of 
equations  V-l  and  found 


*  When  the  subscript  I  does  not  accompany  the  stress  intensity  factor,  open¬ 
ing  mode  is  assumed.  The  units  are  normally  ksi  yinch,  but  may  be  seen 
as  lbs- inches-^/^ . 
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Kj  =  Const.  Or  yj a~ 


.  •  Mode  I 


( V- 22 ) 


where 


K 


II 


=  K 


III 


=  0 


Similarly,  Modes  II  and  III  may  be  analyzed  by  dimensional  analysis. 

KII  =  Const#  Ty/IT 


Mode  II . ( V-23 ) 


where 


ki  ■  Km  ■  ° 


for  uniform  in-plane  shear,  and 

=  Const.  Tyfa~ 


where 


Ki  -  Ku  - 0 


•  *  Mode  III 


(V- 24) 


for  out-of-plane  shear  at  infinity. 

The  value  of  the  constant  (Const.)  is  found  as  yflT for  all  three  modes. 
Therefore,  the  stress  intensity  factors  for  an  infinite  plate  with  the  crack 
geometry  of  Figure  7  are 


K  =  Or  ^  7Ta 


Kn  * 


(V- 25) 


Km  =  r/iT 

using  both  dimensional  and  elastic  stress  solutions. 


The  stresses  in  linear  elasticity  depend  linearly  upon  applied  load, 
hence  Equations  V-16  imply  that  the  stress  intensity  factors  contain  load  as  a 
linear  factor.  Dimensional  examination  of  the  same  equations  show  that  K 
must  also  contain  a  characteristic  length  parameter  (lineal  dimension,  etc.) 
for  a  particular  body  configuration,  including  crack  length.  In  the  Griffith 
crack  configuration  (see  Section  II)  the  only  characteristic  length  is,  a, 
the  crack  length.  The  similarity  between  Equations  V-25  and  the  Griffith 
energy  criteria  for  crack  stability  (see  Equation  II-7)  indicate  that  instability 
occur  at  a  constant  value  of  K. 


There  are  many  other  means  of  determining  stress  intensity  factors  (i.e., 
complex  stress  functions,  numerical  collocation,  inner-outer  expansion)  which 
will  not  be  included  here,  but  lead  to  the  same  dependency  on  crack  length  and 
stress  as  eqs.  V-25.  An  excellent  summary  of  stress  intensity  factor  solutions 
is  contained  in  Reference  15. 
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V.6  STRESS  INTENSITY  -  STRESS  CONCENTRATION 


Similar  terms  and  symbols--dissimilar  parameters  but,  the  stress  inten¬ 
sity  factor  (K)  may  be  obtained  as  a  limiting  case  of  the  stress  concentration 
factor  (Kt ) ,  i.e.,  when  the  notch  radius,  p  ,  approaches  zero. 

For  a  crack  of  finite  tip  radius,  p  ,  the  surrounding  stress  state  is 
governed  by  the  stress  intensity  factor.  A  proportionality  then  exists  be¬ 
tween  Kj  and  O  m,  where  CTm  is  the  maximum  stress  near  the  crack  tip  due  to 
the  presence  of  the  stress  raiser  (crack).  In  order  to  make  this  propor¬ 
tionality  dimensionally  correct,  the  notch  radius  must  be  introduced;  hence, 
the  complete  solution  is: 


K_  =  Const.  O’  fp~ 
I  m  v  r 


( V-26 ) 


It  can  be  recognized  that  this  stress  intensity  solution  is  the  same  as  that 
for  a  Mode  I  crack  (see  equations  V-25)  where  "a"  is  replaced  by  the  charac¬ 
teristic  crack  tip  or  notch  radius,  P  .  The  Kt  for  an  elliptical  notch 
(half  crack  conditions)  can  be  used  as  the  limiting  case  for  a  crack,  or  when 
the  ellipse  collapses,  i.e.,  flattened  and  approaches  a  crack  shape. 


Kt(for  elliptical  notch) 


1  + 


( V- 2  7  ) 


Solving  equation  V-27  for  and  substituting  into  the  equated  Kj  solutions 
of  equations  V-26  and  V-25,  when  P  approaches  zero  (crack  situation),  equa¬ 
tion  V-26  becomes 

K  =  o^Lo  °m~  ^  (V-28 ) 

I  p  — ►  U  m  2 

and  the  constant  (Const.)  is  found  to  be  — — ■  for  P  =  0.  By  this  procedure, 
the  stress  concentration  analysis  may  be  used  to  determine  Kj  mode  stress 
intensities  for  various  shaped  flaws  in  finite  structure  as  long  as  the  stress 
concentration  factors  are  available  in  a  form  similar  to  Equation  V-27. 


Another  interesting  feature  of  this  analysis  is  that  the  constant  (Const, 
was  found  to  be  \ 111  .  This  factor,  when  combined  with  the  finite  width  correc 

tion  (usually  signified  by  the  symbol  Y*)  is  dependent  on  both  specimen  and 
crack  geometry.  The  finite  geometry  corrections  will  be  discussed  later  for 
various  crack  geometries. 

V. 7  THE  RELATIONSHIP  BETWEEN  STRAIN  ENERGY  RELEASE  RATE  Q#)  AND  STRESS 

INTENSITY 

It  can  be  observed  that  the  relationship  between  strain  energy  release 
rate  (jfl )  and  stress  intensity  (K)  as  given  by  Equations  II-7  and  V-25  are 
quite  similar.  Both  are  measures  of  the  elastic  stress  field  effects  near 
the  crack  tip.  It  has  been  shown(^)  that  the  relationship  for  displacement 
Modes  I,  II,  and  III  is 

*The  correction  factor  (Y)  usually  includes  the  \ term  in  its  calcu¬ 
lation.  In  this  report,  (A)  will  be  used  to  denote  finite  width  correction 
which  does  not  include  the 'JW  . 
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£t.=  2 


II 


II 


A 


=  (1  +  v  )  2 

III  E  III 


(Plane  Strain) 


(Plane  Strain) 


.  (V-29) 


and  for  plane  stress  conditions 


(Plane 


Stress ) 


(V-30) 


where  V  and  E  are  Poisson's  ratio  and  modulus  of  elasticity,  respectively. 
Although  most  of  the  modern  terminology  of  fracture  mechanics  uses  the  stress 
intensity  terminology,  the  relationships  of  equations  V-29  and  V-30  can  be 
used  to  convert  strain  energy  release  rate  data  to  K  values  directly. 


V . 8  SUPERPOSITION  OF  K  SOLUTIONS 


In  many  cases  for  complicated  loading  problems,  it  may  be  expedient  to 
add  the  stress  intensity  solutions  in  Mode  I  only,  to  solve  for  a  combined 
loading  problem.  For  example,  in  linear  elasticity  the  effects  of  two  or 
more  loads  may  be  superimposed  by  adding  stresses  and  displacements  (not 
energy)  at  each  point.  In  fracture  mechanics,  a  problem  may  arise  such  as 
shown  in  Figure  8  ,  which  might  represent  a  crack  in  a  structure  under  uni¬ 
axial  stress  with  superimposed  compressive  rivet  forces.*  The  stress  inten¬ 
sity  factor  for  the  total  problem  is  then  the  algebraic  summation  of  the  two 
simpler  problems. 

Superposition  principles  also  apply  to  the  geometric  corrections  (  ATs) 
and  the  principles  involved  will  be  developed  in  later  sections  of  this 
report. 

V . 9  CRITICAL  K  CONCEPT 


The  stress  equations  of  eq.  V-l  indicate  that  the  stress  conditions  near 
the  crack  tip  are  of  similar  distributions  (1/ J r  )  and  vary  only  in  stress 
intensity  from  one  case  to  another.  Therefore,  it  follows  that  unstable 
crack  extension  will  take  place  when  the  stress  intensity  factor,  ,  reaches 
a  value  which  is  critical,  that  is  when  any  combination  of  load,  crack  geometry, 
and  configuration  becomes  critical.  This  situation  occurs  when  K-j-  approaches 
K^c,  where  the  subscript  (c)  denotes  a  critical  value  in  fracture  mechanics 
terminology.  This  critical  value, K^c , is  the  plane  strain  fracture  toughness. 

For  plane  stress  fracture,  the  symbol  Kc  is  used. 

*Since  Mode  II  and  III  distribution  are  entirely  different  from  Mode  I,  K 
factors  cannot  be  added  for  different  modes  (i.e.,  Kx  +  Kn  f  KtoTAl)* 
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FIGURE  8  EXAMPLE  OF  STRESS  INTENSITY  SUPERPOSITION 


In  the  case  of  plane  strain  fracture  toughness  (Kxc)*  the  following 
assumptions  apply: 

1.  A  small  zone  of  plasticity  occurs  at  the  crack  tip. 

2.  Griffith  conditions  (brittle  behavior)  predominate. 

However,  in  a  given  material  the  elastic  stress  field  is  always  changed 
due  to  the  non-linearity  at  the  crack  tip;  but  always  in  the  same  way.  There¬ 
fore,  fracture  will  occur  at  what  is  termed  an  apparent  stress  intensity,  Kj;c. 
Critical  conditions  still  apply,  however,  (K^  'Kxc)  for  materials  which  show 
limited  ductile  behavior. 

V.10  SUMMARY 


In  this  section  the  general  forms  of  the  stress  and  displacement  fields 
surrounding  a  crack  were  considered.  An  introduction  to  the  form  of  the  stress 
intensity  factor  K  for  the  various  crack  displacement  modes  and  the  formula¬ 
tion  of  K  factors  from  stress  concentration  factors  was  also  presented. 

Up  to  now,  all  crack  and  fracture  criterion  have  been  based  on  linear 
elastic  theory;  that  is,  plastic  constraint  confined  to  a  small  area  near  the 
crack  tip.  With  these  fundamentals,  we  can  now  approach  the  real  life,  duc¬ 
tile  material  situation  where  plasticity  takes  place  and  requires  modification 
to  the  elastic  solutions. 


28 


VI  CRACK  TIP  PLASTICITY 


VI. I  ESTIMATES  OF  PLASTIC  ZONE  SIZE 


In  most  real  engineering  materials,  there  exists  near  the  crack  tip, 
under  increasing  load,  a  zone  of  plasticity  (plastic  zone)  which  varies  in 
size  and  shape  for  states  of  stress.  The  distribution  and  size  of  this  zone 
has  been  studied  by  several  investigators .  ( 17-20)  Due  to  highly  complex 

analysis  required  to  solve  the  plasticity  equations  in  Modes  I  and  II,  a 
solution  has  only  been  obtained  for  Mode  III. (21)  An  estimate  can  be  made, 
however,  for  Mode  I  cracks  if  it  is  assumed  that  the  elastic  stress  field 
equations  V-l  approach  yield  conditions  as  the  elastic-plastic  boundary  is 
reached.*  Using  Octahedral  Shear  theory,  the  relative  size  of  the  plane 
stress  and  plane  strain  plastic  zone  (r) is  found  to  be 


where 


r 


r 

y 

the  plastic 


2  $  ,i  -  .20, 

cos  —  (1  +  3  sin  —  ) 

2  0  [,  ,  .2  ,  ,  .20 

cos  —  (1  -  2  V)  +3  sin  — 

zone  radius,  ry,  is  given  by 


(Plane  Stress) 


(Plane  Strain) 


(VI-1) 


r 

y 


(Plane  Stress) 


(VI-2) 


and  the  actual  plastic  zone  diameter  is  2ry 

The  plane  stress  equation  VI-1  indicates  the  following: 

1.  The  distance  ry  is  the  apparent  crack  tip  location  directly  ahead  of 
the  actual  crack  tip  and  prescribes  the  elastic-plastic  boundary 
(plane  stress  plastic  zone  radius). 

2.  The  plastic  zone  is  larger  than  that  given  for  plane  strain. 

Liu, (2^0  in  1961,  proposed  a  model  of  the  plastic  zone  which  is  depicted 
in  Figure  9.  Experimental  evidence  has  shown  that  this  model  is  a  good 
representation  of  the  zones  for  plane  stress  and  plane  strain.  The  relative 
sizes  of  the  plastic  zones  are  as  indicated. 

Dugdale^^  analytically  described  the  shape  of  the  plastic  zone  on  the 
surface  as  "wedge”  shaped  or  strip  model,  and  these  types  of  zones  have  been 
observed  experimental ly ( 1 7 > 19 ) •  others  propose  a  "dumb-bell"  shape  during 
the  early  stages  of  deformation.  Many  feel  that  the  Dugdale  equation  esti¬ 
mates  the  plastic  zone  size  more  accurately  for  certain  materials  than 


*  This  method  becomes  "exact"  as  the  material  stress-strain  curve  becomes 
linear. 
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FIGURE  9  PICTORIAL  REPRESENTATION  OF  THROUGH-THE-THICKNESS  PLASTIC 
ZONE  DEVELOPMENT  AS  FIRST  VISUALIZED  BY  LIU  (Ref.  20) 
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that  given  by  equation  VI-2.  This  equation  assumes 
length  can  be  described  mathematically  as 


2r  = 

y 


a  wedge  shaped  zone  whose 


(VI-3 ) 


which  in  terms  of  the  stress  intensity  factor  is  (when  ry  <  <  a) 


77  Kt 


2r  =  •? 

y  2 


2  a 


ys 


(VI-4) 


The  reason  for  these  models  will  be  discussed  in  the  next  section. 
However,  what  is  important  is  that  at  the  crack  tip  a  plastic  zone  develops 
under  increasing  stress  and  affects  the  stress  intensity  solution. 


VI. 2  LIMITATIONS  ON  PLASTIC  ZONE  DIMENSIONS 

The  plastic  zone  must  remain  small  compared  to  other  dimensions  of  the 
material  (i.e.,  crack  length,  width),  or  elastic  fracture  mechanics  is  in¬ 
validated. 


If  equation  VI-2  is 
Mode  I  (see  eqs.  V-25), 
the  applied  gross  stress 


combined  with  the  stress  intensity  equation  for 
the  plastic  zone  size  can  be  seen  as  a  function  of 
and  crack  length.  (Also  see  equation  VI-3.) 


(VI-5) 


The  plastic  zone  size  given  by  equations  VI-5  or  VI-3  will  be  small  com¬ 
pared  to  half-crack  length,  "a",  as  long  as  the  gross  stress  is  well  below  the 
material  yield  strength.  This  relationship  places  some  limitation  on  the  ratio 
of  O/o  ,  which  was  reported  in  the  fifth  ASTM  Special  Committee  Report(23) 
as  yS  a 

net  <  0.8  (VI-6) 


a 


ys 


where  (7  t  is  the  stress  on  the  uncracked  section  at  fracture.  Thus,  when 
the  size  of  plastic  zone  is  a  large  portion  of  the  uncracked  net  section, 
linear  elasticity  does  not  represent  the  stress  field  with  adequate  accuracy. 
The  relationship  of  equation  VI-6  was  determined  empirically,  and  its  validity 
in  the  high  stress  (small  crack  size)  range  must  be  verified  by  testing  ex¬ 
perience.  This  involves  a  computation  directly  from  equation  VI-2  to  insure 
that  the  plastic  zone  size  is  small  compared  to  other  specimen  dimensions. 


VI. 3  INFLUENCE  OF  WORK  HARDENING  ON  THE  PLASTIC  ZONE 


Work  hardening  affects  the  crack  tip  plasticity  in  the  following  manner: 

Large  Work  Hardening  --  Equations  VI-1  apply  and  "dumbell"  or  circular 

zones  are  formed  (see  Figure  9  ). 

Small  Work  Hardening  --  Concentration  of  slip  ahead  of  the  crack  forms 

shear  bands  and  the  Dugdale^2)  »»WedgeM  model 
applies  (see  Figure  10  and  equation  VI-4). 

Figure  10  shows  the  shapes  of  experimentally  observed  plane  stress  plastic 
zones  for  several  materials  of  various  work  hardening  coefficients,  n,  as 
summarized  by  Hahn  and  Rosenf ield (1 7) .  Reference  17  contains  an  excellent 
summary  of  the  influence  of  strain  hardening  on  fracture. 
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FIGURE  10.  EXPERIMENTALLY  OBSERVED  SURFACE, 
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PLANE  STRESS  PLASTIC  ZONES 
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Another  illustration  of  the  plane  stress  plastic  zone  is  shown  in 
Figure  11.  It  can  be  observed  that  the  so-called  V-slant  fracture  surface 
will  occur  due  to  this  type  of  plastic  zone.  Slip  within  the  plastic  zone 
(at  45°  to  crack  tip  plane)  is  one  of  the  causes  of  this  type  of  fracture. 

It  will  be  remembered  that  the  state  of  plane  stress  or  plane  strain  also 
influences  the  formation  of  slant,  V-slant,  and  flat  fracture  (see  Figure  3). 


VI. 4  THICKNESS  EFFECTS  ON  FRACTURE 


To  differentiate  between  the  terms  thick  and  thin,  the  radius  of  the 
plane  stress  plastic  zone  (eq.  VI-2)  can  provide  a  relative  measure  when  com¬ 
pared  to  material  thickness,  B.  Test  results  have  indicated  that  the  frac¬ 
ture  appearance  (through  the  thickness)  is  related  to  a  normalized  plastic 
zone  size,  as  shown  in  Figure  12.  This  behavior  is  typical  and  indicates 
the  transition  zones  as  a  function  of  the  crack  tip  parameter,  .  The 
shape  of  this  curve  will  depend  on  many  parameters,  and  the  fracture  mode 
transition  on  the  interplay  between  material  properties,  geometry,  and  force. 
However,  as  a  general  rule,  the  cut-offs  of  0.2  for  plane  strain  and  1.0  for 
plane  stress  are  generally  accepted.  However,  there  is  one  material,  beryl¬ 
lium,  which  does  not  show  a  transition  from  flat  to  shear. 

This  brings  us  to  the  effect  of  material  thickness  (B)  on  the  critical 
stress  intensity  factor,  Kc.  Figure  13  shows  these  trends  for  titanium, 
steel,  and  aluminum  alloys.  The  form  is  similar  in  many  cases  to  the  curve 
of  Figure  12.  For  example,  the  nominal*  value  of  stress  intensity  for  7075- 
T6  has  a  minimum  at  36  ksi^/inch  for  >  0.8  inch  thick  material.  This  limiting 
value  represents  the  plane  strain  fracture  toughness,  Kjc,  for  7075-T6.  A 
mixed  mode  region  occurs  for  0.5  >B  >0.15  inches ,  and  the  plane  stress  region 
occurs  when  B<0.15  inches  and  Kc^85  ksi^inch.  The  other  curves  of  Figure 
13  show  similar  trends,  some  more  gradual  (Ti-6A1-4V)  and  others  more  drastic 
(Ti -B120VCA) . 


There  have  been  several  indexes  proposed  to  measure  the  tendency  toward 
plane  strain  or  plane  stress  behavior.  In  a  plate,  the  plastic  zone  radius 
at  the  surface  is  given  by  the  plane  stress  equation  VI-2  or,  VI-4;  and  in  the 
center  of  the  plate  by  the  plane  strain  equation  VI-7,  which  is  1/ (2^2)  or 
approximately  1/3  the  size  of  the  plane  stress  zone  equation. 


r 

y 


Plane  Strain 


(VI-7) 


The  critical  conditions  for  plane  strain  and  plane  stress  can  be  estimated  by 
comparing  the  surface  plastic  zone  radius,  eq.  VI-2,  with  thickness,  or  when 


plane  strain  conditions  prevail.  And  when 


*CAUTION--The  critical  stress  intensity  values  shown  in  Figure  13  must  be 
considered  as  average. 
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cr 
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FIGURE  II.  V-SLANT,  PLANE  STRESS  PLASTIC  ZONE  AS  OBSERVED  BY 
HAHN  6c  ROSENFIELD(l7) 
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NORMAL  FRACTURE  APPEARANCE 


PLANE  STRAIN 


r 

Pm 

t 


S 

M 

s 


0  0.2  0.4  0.6  0.8  1.0 

RADIUS  OF  PLASTIC  ZONE.  ^ 

THICKNESS  B 


FIGURE  12.  TREND  IN  FRACTURE  MODE  APPEARANCE  VS. 

CRACK  TIP  PLASTIC  ZONE  PARAMETER 
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CRITICAL  STRESS  INTENSITY  FACTOR,  K  (ksi  ^ inch ) 


240n 


0 


0 


0.2 


“T““  -1— 

0.4  0.6  0.8 

THICKNESS,  B  (inches) 


1.0 


FIGURE  13.  NOMINAL*  CRITICAL  STRESS  INTENSITIES  FOR  SEVERAL 
MATERIALS  AS  A  FUNCTION  OF  THICKNESS  (Data  from 
Refs.  24  and  23. ) 


*CAUTION  -  Values  should  be  considered  as  average . 
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B  <  «  2 


2 


plane  stress  conditions  predominate. 


Relationships  such  as  the  above  can  be  useful  in  setting  a  limit  on 
specimen  size  for  determining  Kxc.  The  American  Society  for  Testing  and 
Materials(®)  has  set  a  requirement  that 


B  and  a  — 


(VI-8) 


in  order  to  obtain  valid  plane  strain  fracture  toughness  values,  under  current 
accepted  practice.  However,  since  these  are  empirical  estimates,  the  exact 
method  of  determining  Kjc  for  a  particular  material  is  to  increase  specimen 
thickness  until  critical  K  values  reach  a  limiting  value  (see,  e.g.,  7075-T6, 
Figure  13). 


VI. 5  INCLUSION  OF  PLASTICITY  AND  GEOMETRIC  CORRECTIONS  IN  K  SOLUTION 

It  was  shown  that  the  radius  of  plastic  zone,  ry,  affects  the  material 
critical  stress  intensity  factor.  Therefore,  it  must  be  accounted  for  in 
some  manner  in  the  basic  K  equations  V-22,  23,  and  24.  It  may  be  thought  of 
as  adding  an  incremental  length  to  the  existing  crack.  For  example,  equation 
V-22  becomes 

Kj  *  X(Jy/  a  +  ~f  (r  )  (VI-9) 

/  K.  \2 

where  the  plastic  zone  radius,  ry,  is  a  function  of  I  j ■  )  from  eq.  VI-2  or  7, 

and  X  is  the  correction  for  finite  geometry.  \  / 

In  a  previous  section,  equations  V-25  gave  the  K  solutions  for  an  infinite 
panel  containing  a  central  crack  of  half-length  "a".  For  centrally  cracked 
finite  panels.  Irwin'**'  suggested  the  following  equation  for  stress  intensity 
from  the  work  of  Westergaard,  including  geometric  correction,  to  account  for 
finite  specimen  widths,  (W)  and  crack  geometry. 

K  =  oj  W  tan  ^2.  (VI-10) 

This  is  the  so-called  "Irwin  Tangent  Formula"  where  the  factor  W  tan 
includes  the  geometric  correction  factor  (X)  discussed  previously.  The  form 
of  equation  V-25  accounting  both  for  plasticity  and  finite  geometry  using  the 
"Irwin  Tangent  Formula"  becomes 

Kj  =  O  s/  W  tan  +  f  (r  ))  (VI-11) 

where  the  plastic  zone  size  is  determined  from  equations  VI-2,3,  or  7  depending 
on  stress  state  and  expected  zone  shape. 
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The  normal  form  of  the  opening  mode  stress  intensity,  including  finite  width 
correction  is  usually  written  as 


or 


(VI-12) 


for  centrally  cracked  panels.  It  has  been  determined  that  the  nIrwin  Tangent 
Correction"  factor  has  an  estimated  accuracy  of  ±5%  for  crack  length  to  panel 
width  ratios  (— j''* **less  than  0.4.  Other  correction  factors  have  been  developed 
for  the  center W  crack  geometry (26-28) #  t he  expression  in  use  today  is  by  Isida, 
which  has  an  accuracy  of  ±0.5%  for  1*  <0.7.  The  forms  of  the  various  correc¬ 
tion  factors  which  account  for  relative  size  of  specimen  to  crack  dimensions, 
will  be  discussed  fully  in  Section  VII. 6. 


VI. 6  SUMMARY 


The  plastic  zone  at  the  crack  tip  has  been  shown  to  affect  the  state  of 
stress,  fracture  mode  appearance,  and  most  importantly,  critical  stress  in¬ 
tensity  factor.  It  is  itself,  in  turn,  a  function  of  thickness  and  material 
properties.  Estimates  can  be  made  of  the  plastic  zone  size  and  result  in  an 
apparent  crack  length  which  can  be  used  to  determine  stress  intensities  for 
semi-ductile  fracture.  The  complete  form  of  stress  intensity  equation  will 
include  corrections  for  plate/crack  geometry  and  for  plastic  zone. 

In  the  following  section,  the  stress  intensity  solutions  will  be  pre¬ 
sented  for  various  crack  loading  situations. 


*  Plastic  zone  correction  not  included. 

**Also  referred  to  as  "crack  aspect  ratio", 


/TOTAL  CRACK  LENGTH  \ 

l  PANEL  WIDTH  J 
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VII  BASIC  STRESS  INTENSITY  FACTORS 


This  section  presents  solutions  of  the  basic  stress  intensity  factor  K 
for  crack  geometries  and  loading  modes  most  encountered  in  aircraft  applica¬ 
tions.  These  expressions  have  been  abstracted  from  Reference  15  which  con¬ 
tains  an  excellent  review  of  available  K  solutions.  The  original  sources 
are  also  given  in  Reference  15.  Using  superposition  principles,  it  will  be 
possible  to  combine  many  of  the  results.  In  most  cases  the  opening,  Mode  I, 
tensile  stress  intensity  will  be  emphasized  since  this  is  the  most  common 
case  for  aircraft  application. 

NOTE:  ALL  K  EQUATIONS  IN  THIS  SECTION  ARE  BASED  ON  ELASTIC  SOLUTIONS  AND 
ARE  FOR  INFINITE  AND  SEMI -INFINITE  GEOMETRY  (i.e.,  without 
geometric  (finite  width)  corrections).* 

VII. 1  UNIFORM  LOADING 

Uniform,  Tensile  Stress  (Normal  to  Crack) 


Kj  =  a  yira"  (vii-i ) 


(Mode  I) 


Uniform  Shear  Stress 


with 


Kn  =  r/^T  (VII  -  1A ) 


(Mode  II) 


K 


I 


K 


III 


=  0 


*Except  Where  Noted. 
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Uniform,  Transverse  Shear  Stress  (Parallel  to  Crack) 


Uniform,  Twisting  Moments 


F 


Uniform,  Shear 


(Although  the  K  solutions  are  consistent 
with  Kt 


K 


Kfi  the  standard  notation  is 


shear’  lending  for  flexure‘> 


F  =  Pure  Twisting  Moment,  (per  unit  length) 


6  F  JH a~ 

^shear - <VII-2> 

E 


with  K,  =  0 

bending 


F  =  Uniform  Shear  (per  unit  length) 

NOTE:  Results  are  independent  of  F* 


K 


shear 


(VII-3) 


with  K,  =  0 

bending 
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Uniform,  Bending 


M  =  Uniform  Bending  Moment 

NOTE:  Results  are  independent  of  M* 
(Equilibrium  Moments) 


ending 


(VII-4) 


with  K  .  =  0 

shear 


Uniform  Tensile  Stress  -  Crack  Arrays 


(Mode  I) 


J 


P  =  Force  per  unit  thickness 
Kj  @  points  A  = 


with  KXI  =  0 


Uniform,  In  Plane  Tensile  Stress  -  Inclined  Crack 
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Uniform  Tensile  Stress  -  Corner  Crack 


(Mode  I) 

Kj  =  MR  a /FT  (VI1-7 ) 


e 

90° 

0 

0.705 

0.2 

0.650 

0.4 

0.625 

0.5 

0.620 

0.6 

0.625 

0.8 

0.650 

1.0 

0.705 

where  values  of  the  elastic  magnification  factor  M,  for  values  of  angular 

9 

orientation  qqo  are  as  above  (Reference  25,  see  also  Section  IX). 
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In-Plane  Bending 


KI  = 


(Mode  I ) 
6M 


BW 


or  K  =  °  NOM  yWa" 


M  =  bending  moment 
°"nOM  =  nominal  stress 


with 


ku  -  Km  -  ° 


Axial  Tension  (Double  Crack) 


(Mode  I ) 
Kx  =  CT/il 


with  Ku  =  KnI  ,  0 


(VII-8 ) 


(VII-9) 
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STRESS  INTENSITY  FACTOR  COEFFICIENT,  ^(f ) 


Uniaxial  or  Biaxial  Stressed  -  Cracks  at  Holes 
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Uniform  Stresses  on  Part  of  Crack  Surface 


(Mode  I) 

O  =  Uniform  stress  on  one  surface 
T  —  Uniform  Shear  on  crack  faces 


Kj.  @  A 


_a_  /r 

2  V  77 


4  -1  b 

—■sin  — 
a 


.  -1  c 
sin  — 
a 


1  +  V 


.1 1  _,/£ 

\ 2 

V 

/ 

(' 

,2 

+ 1 

v/1- (-■ 

1 

4_  . 

T  (c-b) 

=m 

2  yir 

\  k+1  / 

4  V 

(Plane 

Strain) 

V 

V 

(Plane 

Stress ) 

(VII-  II) 


VII. 2  CONCENTRATED  FORCES 

Concentrated  Force  (F),  on  a  Crack  Surface 

with  Arbitrary  Inclination 

(Mode  I) 


F  =  Force  per  unit  thickness 


(VII-12) 


p 

ii 

2 

y^ 

V  k+1  ) 

q 

/a+b 

+  2 

y^a 

y  a-b 

3 

-  4  y 

(Plane  Strain) 

3 

1 

-  V 

+  v 

(PI 

ane  Stress) 
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Concentrated  Forces  on  Crack  Surface 


(Mode  I) 

P  =  Force  per  unit  thickness 
Q  =  Force  per  unit  thickness 
K 


(Mode  II) 


Concentrated  Forces  on  Crack  Surface 

(Colinear  Cracks)  vMpjie_ 12. 


P  =  Force  per  unit  thickness 
Q  =  Force  per  unit  thickness 
(Mode  I) 


ri - T 

^  p  VC  -  a 

K_  @  A  *  r  , 7 - r 

1  2  Jlfa  (C  -  a) 


(Mode  II) 


J2  2 

Vc  -  a 


KH0A=2^(c-  a) 


(Mode_I ) 


n 2 

p  vc  —  a 

KJB - — - 

1  2  y/TTa  (c  +  a) 


(Mode  II  ) 

^  q  Vc2  - 

2  y/TTa  (c  +  a)y 


@  A 


(a  B 


(VII-13) 


(VII-14) 


(VII-15) 
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Splitine  F orco  ol  rr~.« 


:  ;  .ngi;  l  .*r  bee  t  ion  (a>>2c) 

(Mode  I ) 

V  Concentrated  Force  per  unit  thickness 


K 


I 


3.464  Pa 


(VII- 16 ) 


with  K  =  0 


VII.  3  THERMAL  I'TL'CLD  "TRECEE5 


Thermal  Load  on  Crack  Surface  (Upper  L  Lower 
Crack  Surfaces  at  Same  Temperature) 


(Mode  I ) 


k  3  -  4^  (Plane  Stress) 
i-  (3  —  v)  (Plane  Strain) 

a  +v) 


(Vll-17 ) 


&  —  Coefficient  of  Expansion 

Lt  Thermal  Conductivity 

q  Equivalent  Heat  Flow  Required  to 

Produce  Constant  Tempera. ure  on 
the  Crack  Surfaces 

E  Young’s  Modulus 
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VII. 4  STRESS  INTENSITIES  FOR  KNOWN  STRESS  CONCENTRATIONS 


cr 


tt  t  t 


1111 


cr 


(Mode  I) 


with  K 


K 


III 


=  0 


(from  Eq.  V-28) 
(VII-18 ) 


(Mode  II) 


KTT  =  lim  &  n JW  Jp 
II  p — ►  o  m  v 

with  Kj  =  K  =  0 


(VII-18A) 


(Mode  III) 


K 


III 


lim 
P—  O 


m 


JF  JP 


(VII-18B) 


with  Kj.  -  K  =  0 

where  (Tm  is  the  maximum  normal  stress 

at  the  root  of  the  notch  caused  by  axial 
tensile  stress,  O,  and  r  is  the  maximum 
shear  stress  adjacent  to  trie  notch. 

(See  above  for  loading  conditions.) 


48 


VII.  5  PART  THROUGH  C.RACK 


Uniform  Tensile  Stress  -  Semielliptical 


rr 


/  /  / 


(T 


(Mode  I) 


with  Kjj  =  0 

Other  examples  of  part  through  cracks  will  be  presented  in  Section  IX. 


49 


VII.  6  FINITE  WIDTH  CORRECTION 


Since  the  stress  intensity  solutions  previously  presented  in  this  section 
were  based  on  infinite  geometry,  their  usage  for  finite  panel  problems  must  be 
modified  by  finite  width  correction  factors  (see,  e.g.,  eq.  VI-9  to  VI-12). 

With  the  use  of  these  factors,  fracture  toughness  values  can  be  used  to  provide 
design  criteria. 

The  width  correction  factor  will  be  signified  by  the  symbol  X  (see  Section 
V).  For  a  center  cracked  sheet,  seven  solutions  have  been  in  use  throughout 
the  years.  The  ones  we  will  concern  ourselves  with  will  be  those  currently  in 
accepted  usage.  Finite  width  correction  factors  being  analytic,  lend  them¬ 
selves  to  graphic  presentations.  Presented  in  Figures  14  through  18  are  the 
X  solutions  for  various  crack  geometries  as  a  function  of  crack  length  to 
panel  width  ratio  (crack  aspect  ratio).  NOTE:  These  corrections  have  been 
proven  accurate  as  long  as  the  stress  O  is  remotely  applied,  i.e.,  panel  length 
>  >  W.  In  a  center  cracked  panel,  length  >  3W. 

Geometric  corrections  (or  K  calibrations)  for  current,  standard  fracture 
mechanics  specimens  are  presented  in  polynominal  form  in  Reference  8.  These 
may  be  used  to  supplement  the  graphical  presentations  for  greater  accuracy 
for  various  crack  loading  situations. 


VII. 7  SLOW  TEAR  IN  CRACKED  STRUCTURE 


During  static  stressing  of  cracked,  ductile  materials,  the  crack 
length  increases  from  an  initial  length  aQ,  until  a  critical  length  ac  is 
reached  (see  Figure  1).  Many  times  in  design  problems  it  is  necessary  to  deal 
with  the  initial  crack  size  and  estimates  must  be  made  of  critical  crack  lengths* 
Several  empirical  estimates  of  slow  crack  growth  to  a  critical  value,  ac,  have 
been  presented  (Ref.  31,  32)  which  have  the  form 


ac  =  (1  +  const.)  aQ  (VIl-20) 

where  the  constant  (const.)  provides  the  best  data  fit.  Newman^^  proposed 
that,  as  a  best  data  fit,  the  following  equation  applies: 

(const. ) 

ac  =  1  +  (const. )x |  l! 


(VII-21) 


where  the  two  constants  were  determined  from  experimental  data.  One  shortcoming 
of  equation  VII-20  is  that  when  aQ  »  0,  or  for  very  small  crack  lengths, 

ac - -0.  In  ductile  materials  one  would  expect  some  slow  tear  to  occur  from 

even  the  smallest  discontinuity.  For  long  initial  crack  lengths  >  0.2), 
equations  VII-20  or  VII-21  can  be  used  with  confidence  in  predicting  the 
critical  crack  lengths.  It  is  recommended  that  equation  VII-21  be  used  for 
greater  accuracy  for  predicting  critical  cracks  of  engineering  size. 

The  slow  tear  behavior  of  thin  gage  material,  where  the  plane  stress  fracture 
mode  prevails  is  currently  under  study  by  a  special  task  force  of  the  American 
Society  for  Testing  and  Materials  (Committee  E-24)  on  Fracture  Testing  of  Metals. 
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FIGURE  14.  FINITE  WIDTH  CORRECTION-CENTER  FIGURE  15.  FINITE  WIDTH  CORRECTION-DOUBLE 

CRACKS  (TENSION)  (REF.  29)  EDGE  CRACKS  (TENSION)  (REF.  15) 


FINITE  WIDTH  CORRECTION 


6.O-1 


5.0- 


4.0- 


3.0- 


2.0- 


1.0- 


=  cr/^T  X 


Mil 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7 


CRACK  ASPECT  RATIO,  ft 


FIGURE  16.  FINITE  WIDTH  CORRECTION-SINGLE  EDGE 
CRACKS  (TENSION)  (REF.  30) 
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FINITE  WIDTH  CORRECTION, 


FIGURE  17.  FINITE  WIDTH  CORRECTION-SINGLE  EDGE 
CRACK  (BENDING)  (REF.  8) 
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FINITE  WIDTH  CORRECTION 


FIGURE  18  FINITE  WIDTH  CORRECTION-ECCENTRIC 
CRACK  (TENSION)  (REF.  29) 
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VIII  METHODS  OF  RESIDUAL  STRENGTH  ANALYSIS 


In  the  overall  view  of  structural  analysis  for  fracture  prevention,  ini¬ 
tial  efforts  should  be  made  to  analyze  fracture  or  fatigue-sensitive  areas 
at  an  early  stage  of  design.  Figure  19  shows  a  logical  progression  from 
conceptual  design  to  fracture  analysis.  The  interplay  between  material  and 
fracture  parameters  can  be  seen  by  the  discussion  thus  far. 

Many  methods  of  residual  strength  determination  are  available.  Through- 
cracked,  uniaxial  loaded  sheet  specimens  and  structure  can  be  analyzed  using 
what  are  essentially  modifications  of  notch  analysis.  All  analytical  methods 
require  determination  of  a  "material  constant"  as  a  basic  parameter.  In 
this  regard,  they  are  in  basic  agreement  with  the  fracture  mechanics  approach 
The  degree  to  which  these  methods  agree  for  the  through-crack,  thin  sheet, 
residual  strength  problem  will  be  investigated  in  this  section. 

In  fracture  analysis  there  are  several  choices  of  available  concepts. 

In  fail  safe  design  it  will  be  remembered  that  concern  is  for  the  remaining 
strength  in  the  presence  of  an  existing  crack  or  defect.  Residual  strength 
methods  are  based  on  the  occurrence  or  presence  of  such  cracks.  Therefore, 
in  the  early  part  of  this  section,  the  nature  of  the  overall  problem  will  be 
given  which  will  be  solved  by  fracture  analysis  methods.  The  methods  of 
analysis  for  the  simple  case  of  a  center  cracked  sheet  will  then  be  presented 
and  finally  a  detailed  problem  will  be  solved  which  will  encompass  and  com¬ 
pare  all  pertinent  fracture  analysis  methods.  In  this  manner,  the  applica¬ 
bility  and  simplicity  of  the  fracture  mechanics  approach  will  become  evident. 


VIII. 1  FRACTURE  ANALYSIS  CONCEPTS 

VIII. 1(a)  Fracture  Mechanics  (Critical  K  Concepts) 

As  we  have  seen  so  far,  fracture  mechanics  is  based  on  the  stress  analy¬ 
sis  in  the  presence  of  an  existing  flaw  or  crack.  To  analyze  a  center 
cracked  panel,  the  basic  equation  (VI-12)  gives  the  interrelationship  be¬ 
tween  stress,  crack  length,  and  toughness,  which  can  be  solved  for  critical 
conditions;  i.e.,  when  the  stress  intensity  factor  K  reaches  a  critical  value 
the  operating  stress  cr  — ►  (Tc  and  the  crack  length  a — ►ac.  If  the  material 
fracture  toughness  is  known,  then  an  analysis  can  be  made  to  determine  frac¬ 
ture  stress  ( cr c )  for  a  tolerable  crack  size  (ac). 

VIII. 1(b)  Emperically  Derived  Crack  Strength  Plots 

In  many  design  applications,  it  may  suffice  to  use  what  are  termed  resi¬ 
dual  strength  comparison  plots.  Fracture  data  from  many  tests  on  panels  of 
various  widths  are  plotted  in  Figure  20  for  center  cracked  7075-T6  clad 
aluminum  panels  of  a  given  thickness.  One  fact  evident  from  such  plots  is 
the  change  in  fracture  stress  with  panel  width.  One  drawback  to  using  this 
type  of  data  in  analysis  is  the  expense  encountered  in  obtaining  such  data. 
This  method  is  usually  employed  in  material  comparisons  and  does  not  lend 
itself  to  design  application. 
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FIGURE  19.  FRACTURE  ANALYSIS  PROGRESSION  IN  DESIGN 


CRITICAL  FRACTURE  STRESS _ 

MATERIAL  ULTIMATE  TENSILE  STRENGTH 


1.0 


O  3 
b  b 


0.6 


CRACK  BUCKLING  NOT  RESTRAINED 


6  INCHES 
12  INCHES 

24  INCHES 


— I- 

0.6 


a 


CRACK  ASPECT  RATIO,  — 

w 


(DATA  ABSTRACTED  FROM  "THE  EFFECT  OF  FINITE  SPECIMEN  WIDTH  ON 
THE  RESIDUAL  STRENGTH  OF  LIGHT  ALLOY  SHEET,"  BY  D.  BROEK, 
NLR-TR  M. 2152  NATIONAL  AEROSPACE  LABORATORY,  NETHERLANDS, 
SEPT.  1965.) 


FIGURE  20.  RESIDUAL  STRENGTH  AS  A  FUNCTION  OF  PANEL 
WIDTH  -  0.08  INCH  GAGE  7075-T6  CLAD 
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Also,  care  must  be  exercised  in  using  these  data  in  analysis,  since  crack 
buckling  can  occur  (out  of  plane  of  panel)  due  to  compressive  stresses  parallel 
to  the  crack  induced  by  Poisson’s  ratio  effects  and  the  longitudinal  tension 
loading.  These  compressive  stresses  with  subsequent  buckling  can  cause  a  20% 
to  30%  reduction  in  fracture  stress  from  those  panels  where  buckling  was  re¬ 
strained.  Caution  must  be  exercised  as  many  times  this  testing  parameter  is 
neglected  in  the  reported  data. 

VIII. 1(c)  Method  A 

In  1961,  CrichlowOA)  f  irst  presented  an  analysis  method  for  cracked 
panels  which  included  an  effective  width  parameter  that  empirically  corrected 
for  crack  buckling,  Crichlow^^  has  since  employed  the  Notch  Strength 
Analysis  (NSA)  parameter  of  Neuber(36)?  modified  for  his  analysis.  It  will  be 
seen  later  that  several  of  the  residual  strength  analysis  methods  employ  a 
modification  of  the  NSA.  As  with  the  fracture  mechanics  concept,  a  ’’material 
constant”  (Neuber  constant)  must  be  found  from  these;  and  with  knowledge  of 
this  constant,  critical  fracture  stress  can  be  determined. 

VIII. 1(d)  Method  B 

In  an  effort  to  improve  on  the  effective  width  concept  of  Crichlow, 

Chr is tensen-Denke ^  have  proposed  that  the  stress  distribution  in  a  uni- 
axially  loaded,  cracked  sheet  is  equal  to  the  tensile  ultimate  strength  in  the 
effective  width  region  and  is  a  Wes tergaard ^ ^  distribution  beyond  this  region. 
By  integration  and  correction  for  finite  width,  the  formula  for  gross  area 
stress  is 


where  CL  is  an  empirically  determined  ’’material  constant”  which  is  a  function 
of  the  material  plastic  zone  size. 

VIII. 1(e)  Method  C 

Since  this  method(38)  is  the  basis  for  Royal  Aeronautical  Society  Data 
Sheets,  it  will  not  be  repeated  here.  It  is  essentially  based  on  material 
property  (stress-strain  derived)  and  effective  fatigue  crack  tip  radius  para¬ 
meters,  and  predicts  crack  strengths  based  on  these  properties* 

VIII, 1(f)  Method  D 

This  method^^),  which  also  employs  an  effective  crack  radius,  is  the 
predecessor  to  the  Notch  Strength  Analysis  proposed  by  Kuhn,  which  is  outlined 
below. 
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VIII. 1(g)  Method  E,  Notch  Strength  (NSA)  and  Crack  Strength  Analysis  (CSA) 


These  methods  utilize  a  modified  form  of  the  Neuber  theoretical  stress 
concentration  factor.  The  CSA  me thod^^)  ,  first  introduced  by  Kuhn,  uses  this 
modified  stress  concentration  factor  at  fracture  to  calculate  the  net  area 
stress.  A  ’'material  constant”  is  employed  in  the  basic  stress  concentration 
equation  and  can  be  determined  from  a  single  fracture  test.  Also,  the  effects 
of  crack  buckling  can  be  handled  by  the  introduction  of  a  correction  factor. 

In  the  NSA  method,  the  "material  constant"  is  determined  from  stress-strain 
data  which  are  needed  as  a  basis  for  making  crack  strength  predictions. 

VIII. 1(h)  Method  F  (Par t- through  Cracks) 

The  stress  intensity  factor  for  a  crack  which  is  part  way  through  the 
thickness  is  given  by  Irwin(^l)as 


K  =  1.1 


a  /t:  ^ 
v  Q 


which  is  reasonable  accurate  for  crack  depths  <.  50%  of  the  material  thickness. 
The  flaw  shape  parameter,  Q,  is  a  function  of  both  crack  or  flaw  depth  to 
length  ratio  and  operating  stress  to  material  yield  stress  ratio.  Assumptions 
are  made  that  an  elliptical  crack  front  is  formed,  which  usually  occurs  in 
service.  Tiffany,  et  al.,(^2)  have  presented  an  analysis  which  can  be  used 
for  thinner  gage  designs  with  deeper  cracks  in  both  crack  propagation  and 
fracture.  These  are  presented  later  in  a  special  sub-section  of  Section  IX 
devoted  to  the  part-through  crack  problem. 


VIII. 2  SAMPLE  PROBLEM  USING  VARIOUS  FRACTURE  ANALYSIS  METHODS 

MATERIAL  Material:  7075-T6  Aluminum 

DATA  Panel  width:  W  =  20  inches 

Thickness:  B  =  0.100  inch 
Material  yield  strength:  CFyS  =  77  ksi 

Material  ultimate  strength:  O u  =  84  ksi 

Based  on  an  established  routine  inspection  schedule,  it  is  assumed  that 
a  simple  rectangular  aluminum  panel,  uniformly  loaded  in  tension  across  the 
width,  will  develop  a  central  through- the- thickness  crack  (normal  to  loading) 
4  inches  long  before  detection. 

Question:  Will  the  panel  fail  at  the  given  45  ksi  operating  stress  level 

before  the  crack  is  detected  at  the  routine  inspection  interval? 


NOTE:  At  this  point  a  decision  is  made  to  analyze  the  panel  using  a  plane 

stress  toughness  value.  This  decision  is  predicated  on  the  material 
gage  (0.100  inch)  which  should  produce  a  plane  stress  slant  fracture  in 
this  alloy.  Increased  proficiency  and  familiarization  with  fracture 
mechanics  will  make  this  decision  easier.  For  those  materials  which 
have  been  investigated  thoroughly,  this  is  a  simple  task;  however,  other 
materials  may  require  searching  for  valid  fracture  toughness  values. 
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VIII*2(a)  Fracture  Mechanics  Analysis 


From  fracture  toughness  data  such  as  Figure  13,  it  can  be  seen  that  the 
plane  stress  fracture  toughness,  Kc,  is  approximately  83  ksi\/inch  for  7075-T6 
of  this  thickness.  The  20-inch  width  also  dictates  that  valid  Kq  values  would 
probably  be  obtained  for  this  crack  geometry,  (i.e.,  infinite  panel  behavior 
and 

The  infinite  plate  solution  for  stress  intensity,  equation  VII-1,  is 


K 


I 


oJTa 


The  finite  width  corrected  form  from  equation  VI-12  is 


Step  1.  Solve  equation  VI-12  for  gross  stress, 


a  --r±r 

and  let  critical  conditions  apply;  i.e.,  K^. — ►  K^,  a 
Therefore , 


a  ,  and  a — ►a  . 
c’  c 


(VIII-1) 


Step  2.  Determine  the  finite  width  correction  factor,  X  .  For  the  center 
crack  geometry  the  crack  aspect  ratio,  —  =  total  crack  length  £s 
*  7  *  W  specimen  width 


2a  _  4  inches  _  n  ?n 

W  20  inches  # 

and  from  Figure  14,  X  for  this  center  crack  aspect  ratio  is  1.03. 

Step  3.  Solve  equation  VIII- 1  for  the  given  conditions,  assuming  that  a  4-inch 
crack  is  the  critical  condition. 


a 

c 


85  ksiv/inch 

incwrT.03 


33  ksi 


NOTE:  The  critical  stress  value  can  be  reduced  by  207.  to  307.  if  crack  buckl¬ 
ing  takes  place.  The  Kj.  value  reported  here  assumes  an  unbuckled 
panel  configuration. 
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Thus  for  specified  conditions,  the  sheet  can  sustain  a  gross  area  stress 
which  is  less  than  33  ksi  in  the  presence  of  a  4  inch  long  crack.  As  this  is 
277,  below  the  required  operating  stress  (45  ksi),  the  requirements  are  not  met 

From  this  illustration,  a  change  in  inspection  schedule  to  detect  a  crack 
smaller  than  4  inches  or  a  reduction  of  operating  stress  level  or  both  is 
indicated . 

VIII. 2(b)  Crack  Strength  Plot 

For  comparison,  the  crack  strength  plot  of  Figure  20  can  be  used  to  obtain 
a  very  rough  estimate  of  critical  failure  stress  for  this  sample  panel  geo¬ 
metry.  A  critical  stress  of  0.3  <JU  or  25  ksi,  which  is  447,  lower  than  the 
anticipated  operating  stress,  would  be  predicted  using  these  plots.  The  data 
of  Figure  20  are  for  crack-buckled  sheet;  hence,  a  207,  reduction  in  fracture 
stress  over  that  obtained  by  fracture  mechanics  analysis  would  be  expected 
when  no  crack  buckling  effects  are  considered. 


VIII. 2(c)  Method  A  Analysis 

The  equation  given  in  Ref.  35  will  be  used  to  analyze  the  sample  problem, 


O 

c 


„  / 1  2a 

O  l  1  "  TT 
_ u  \  W 

,  s/T  >'  JIT 
4T 


(VII-2) 


where  A’ 

4P7 

e 

and  E 


is  the  Dixon^^)  finite  width  correction  factor 

the  empirically  determined  Neuber  ,!material  constant’1 

ultimate  tensile  elongation 

the  Young's  Modulus  in  ksi  x  10^. 


Step  1. 


Find  the 


finite  width  correction 


factor : ( ^6 ) 


A’ 


=  0.82 


NOTE:  The  critical  stress  value  can  be  reduced  by  207,  to  307,  if  crack  buckl¬ 
ing  takes  place.  The  Kc  value  reported  here  assumes  an  unbuckled 
panel  configuration. 
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Step  2.  Given  the  modulus,  ultimate  elongation,  and  Neuber  constant  for 
7075-T6  of  this  thickness  and  rolling  direction,  solve  equation  VIII-2  for 
gross  fracture  stress: 

E  =  10.3  ksi  x  103 


e  =  137.  =  0.13 

and  J  P  '  from  published  data(43)  =  0. 125  J  inch 


a 

c 


84  ksi 


4  inch 
20  inch 


1  + 


*FT  0.82  inch 

0.125/I^r  (i  +  °-8L<0-"><‘;:3  ksl  x  I<M)J) 


a  =32  ksi 
c 

The  value  of  critical  fracture  stress  predicted  by  this  method  is  approxi¬ 
mately  297o  below  the  operating  stress  level.  If  crack  buckling  is  expected 
to  take  place,  this  value  of  critical  stress  is  reduced  in  the  Method  A  analy¬ 
sis  by  the  empirical  factor  developed  by  Kuhn  and  Figge:^^^  1  -  Y  ^  , 
where  Y  is  another  material  constant  which  has  been  found  to  be  0.001  for 
aluminum  alloys. (35)  This  reduces  the  unbuckled  critical  stress  value  (32 
ksi)  by  47>  and  predicts  a  crack  buckled  critical  stress  of  30.7  ksi,  which  is 
327»  below  the  operating  stress  level. 

VIII. 2(d)  Method  B  Analysis 

In  using  this  method  the  "material  constant",  OL  ,  must  be  determined 
from  a  fracture  test, 


a  = 


2&r 


1  ^ 


2a.  XT 

w 


-  1 


+  1 


vau 

For  7075-T6,  where  O'  =  0.75  inch,  CFc  is  (Ref.  37): 


G  = 
c 


1  +  2ar 


(VIII-3 ) 


84  ksi  | 

4  inch  \ 

20  inch  / 

yj  1+4  inch 

’  i  i 

1 

2  ^ 

(0.75  inch) 

O’  =  32.4  ksi 
c 


for  an  unbuckled  crack  of  the  geometry  of  the  sample  problem,  or  287,  below 
the  expected  operating  stress  level. 
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VIII. 2(e)  Method  E,  Crack  Strength  Analysis  (NSA)  &  (CSA) 

The  net  section  stress  at  fracture  is  given  in  the  NSA^^)  method  as 


<xtl 


a 


net 


1  +  a  Af  Jli 


and  in  the  CSA  method 


(40) 


as 


O 

net 


1  +  a'  A'  JT 


(VIII-4) 


(VIII-5 ) 


The  Oi  term  in  equation  VIII-4  is  a  material  constant  which  assumes  no  notch 
strengthening.  To  obtain  a  better  correlating  equation,  equation  VIII-5  was 
introduced.  It  contains  a  modified  "materials  constant"  (  Oi  f )  and  ultimate 
strength  (O  u  ),  determined  from  material  on  which  a  crack  strength  test  has 
been  run.  We  will  utilize  equation  VIII-4  to  determine  the  net  failure  stress. 
The  value  of  Oi  for  7075-T6  is  1 . 4^nch  ^  from  published  data(^)  for  this 
thickness,  and  X1  is  the  Dixon  finite  width  correction  term  as  utilized  in 
the  Crichlow  analysis.  From  equation  VIII-4, 


O 


net 


84  ksi 

T 


1  +  l^y/inch  (0.82)  yfl inch 


O  .  =  32  ksi  . 
net 

The  critical  gross  stress,  Oc,  for  a  4  inch  crack  would  be  25.6  ksi  or  437® 
lower  than  the  operating  stress  for  an  unbuckled  central  crack.  The  buckled 
crack  geometry  reduces  this  value  again  by  47o  for  a  critical  stress  of  24.6 
ksi. 


A  comparison  of  the  results  from  these  various  analyses  is  provided  in 
Table  I.  The  predicted  fracture  stress  for  both  the  buckled  and  unbuckled 
crack  condition  in  the  sample  geometry  is  indicated  as  a  percentage  of  the 
original  assumed  operating  stress,  45  ksi.  This  comparison  shows  that  the 
agreement  between  methods  is  quite  good  in  the  unrestrained  condition.  The 
NSA  analysis  with  buckling  restrained  would  imply  lower  fracture  stress. 

In  any  case,  the  answer  to  the  question  posed  in  this  sample  problem  is  that 
the  7075-T6  panel  will  not  be  able  to  withstand  an  operating  stress  of  45  ksi 
for  the  crack  geometry  in  question. 
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TABLE  I 

COMPARISON  OF  FRACTURE  DATA  FOR  VARIOUS  ANALYSIS  METHODS 


ANALYSIS  METHOD 

Critical  Stress  and  Percent  of  Operating 
Stress  (45  ksi) 

With  Crack  Buckling 

No  Crack  I 

iuckling 

.  Otsi) 

(%1 

(ksi) 

1  (%) 

Fracture  Mechanics 

26.4-23.1 

(50-58) 

33 

(73) 

Strength  Plot 

25 

(56) 

— 

— 

Method  A  Analysis 

30.7 

(68) 

32 

(71) 

Method  B  Analysis 

— 

— 

32.4 

(72) 

Notch  Strength  Analysis  (NSA) 

24.6 

(55) 

25.6 

(57) 

VIII. 2(f)  Common  Grounds 

In  analyzing  this  simple  case  of  a  cracked  plate  in  tension,  it  becomes 
apparent  that  the  following  conditions  are  common  for  all  the  available  thin 
sheet  analysis  methods: 

1.  A  "constant"  or  "constants"  related  to  some  material  property  must 
be  known  or  determined. 

2.  Finite  width  corrections  apply. 

3.  In-plane  buckling  of  the  sheet  at  the  crack  tip  reduces  the  critical 
fracture  stress. 

4.  All  methods  except  the  fracture  mechanics  analysis  are  based  on 
stress  reduction  from  the  material  ultimate  tensile  strength. 

It  should  be  remembered  that  the  calculations  of  fracture  stress  are 
fundamentally  based  on  the  "material  constant",  be  it  fracture  toughness  or 
notch  factor;  hence,  the  accuracy  will  depend  on  how  this  parameter  is 
determined.  In  fracture  mechanics,  it  has  been  shown  to  be  both  width  and 
thickness  dependent  in  the  plane  stress  fracture  mode.  However,  as  more  Kc 
data  becomes  available,  and  as  a  better  understanding  of  the  factors  involved 
in  general  yielding  ( (7  c  >  (TyS)  becomes  known,  a  better  base  of  confidence 
in  the  fracture  mechanics  analysis  of  thin  sections  will  be  provided.  Similar 
statements  apply  to  the  other  analysis  methods,  which  have  resulted  in  new 
or  modified  "material  constants"  to  better  fit  the  available  data.  The  intro¬ 
duction  of  tougher  new  materials  have  dictated  changes  in  the  fracture  analysis 
and  more  modifications  are  inevitable. 
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IX  INFLUENCE  OF  FRACTURE  IN  AIRCRAFT  DESIGN 


In  conventional  aircraft  structural 
(stress),  "ultimate  load"  (stress),  and 
employed.  For  purposes  of  this  report, 
follows : 


design,  the  terms  "limit  load" 
"factor  of  safety"  are  extensively 
the  terms  may  be  briefly  defined  as 


Limit  load  (stress)  is  the  maximum  load  (stress)  expected  in  operational 
service  (usually  set  by  or  derived  from  structural  design  criteria  speci¬ 
fications,  and  represents  that  value  for  which  the  probability  of  exceeding 
that  load  is  very  small). 

Ultimate  load  (stress)  is  the  limit  load  times  the  factor  of  safety 
(specifications  usually  require  that  the  structure  be  capable  of  sustaining 
ultimate  load,  as  defined;  structure  collapse  just  beyond  ultimate  is 
accepted ) . 

Factor  of  safety  is  the  multiplying  factor  (usually  set  by  specifications) 
to  convert  limit  load  (stress)  to  ultimate  load  (stress).  Customary  USAF 
practice  is  to  employ  a  factor  of  safety  (ultimate  of  1.5). 


The  examples  presented  in  this  section  assume  that  the  factor  of  safety 
is  1.5  and  that  the  tension  stress  encountered  at  limit  load  (limit  design 
stress)  is  2/3  of  the  ultimate  tensile  strength  of  the  material  (2/3  <ru). 

(This  simplification  is  for  illustrative  purposes  herein  only;  it  is  recognized 
that  other  requirements,  e.g.,  fatigue,  stiffness,  etc.,  often  dictate  lower 
values  of  stress  at  limit  and  ultimate  loads.) 

NOTE:  The  examples  presented  in  this  section  are  for  illustrative  purposes 

only.  In  practice  each  problem  will  have  its  unique  solution. 

IX. 1  MATERIAL  TRADE-OFFS 


Quality  control  dictates  that  cracks  must  be  detectable  within  certain 
limits.  As  an  example,  in  a  center  cracked  panel*,  a  1/2-inch  long  crack  can 
be  detected  in  routine  inspections.  The  material  thickness  is  planned  to  be 
0.100  inch.  There  is  a  choice  of  three  materials  which  possess  approximately 
the  same  ultimate  tensile  strength  properties.  The  design  stresses,  mechani¬ 
cal  properties,  and  fracture  toughness  are  listed  in  Table  IX-II.  Infinite 
panel  behavior  is  assumed  in  this  design  (  <0.1)  and  the  finite  width 

correction,  X  ,  from  Figure  14,  is  1.00. 

From  the  tabular  data,  the  2024-T3  alloy  appears  to  be  the  most  suitable 
material  based  on  average  plane  stress  fracture  toughness  (Kc).  The  critical 
gross  fracture  stress,  crc ,  is: 

c  95  ksi  \f  inch _ 

0.25  inch  (1.00) 

(T  -  107  ksi 

c 

The  predicted  critical  fracture  stress  is  2.3  times  the  guaranteed  minimum 
yield  for  this  material.  This  indicates  that  fracture  will  occur  with  gross 
plastic  deformation  and  yielding  on  the  net  section  rather  than  directly 


*Although  this  crack  geometry  occurs  rather  infrequently  in  service,  it  pro¬ 
vides  a  good  illustrative  basis  for  material  fracture  evaluation.  As  such, 
it  will  be  used  here  for  just  such  a  comparison. 
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from  the  influence  of  the  1/2-inch  crack.  It  would  be  acceptable  under  these 
circumstances  to  choose  a  less  frac ture- tough,  but  stronger  alloy  with  a  higher 
limit  design  stress. 

TABLE  II 

DATA  FOR  MATERIAL  TRADE-OFF  STUDY  PROBLEM 


Alloy 

Designation 

(ksi ) 

<r  a) 
<&> 

Kc(2,3,4) 

^lim/5^ 

Limit  Design 
Stress  (ksi) 

°c 

( ks  .1 ) 

O’c 

°Y. 

O’c 

(ksiv/inch) 

O’u 

°~ys 

°"lim. 

2024-T81 

CLAD 

67 

59 

65 

45 

73 

1.1 

1.2 

1.6 

7075-T6 

CLAD 

73 

64 

85 

49 

96 

1.3 

1.5 

2.0 

2024-T3 

CLAD 

63 

46 

95 

42 

107 

1.7 

2.3 

2.5 

NOTES : 


(1)  MINIMUM  VALUES  FROM  MIL-HDBK-5 

(2)  AVERAGE  TEST  VALUES  -  PARTICULAR  DESIGN  WOULD  CALL  FOR  FRACTURE 
TOUGHNESS  VALUES  FOR  0.100 -INCH  MATERIAL  IN  THE  WIDTH  OF  INTEREST. 

(3)  PLASTIC  ZONE  ADJUSTMENTS  HAVE  BEEN  INCLUDED  IN  THESE  VALUES,  AND 
NO  CRACK  BUCKLING  IS  ASSUMED. 

(4)  ENVIRONMENTAL  EFFECTS  HAVE  BEEN  NEGLECTED. 

(5)  BASED  ON  A  F.S.  OF  1.5. 


The  other  possibilities  listed  in  Table  II  are  2024-T81  anc  7075-T6.  The 
critical  fracture  stresses  were  computed  as  for  the  2024-T3  alloy.  The  ratios 
of  critical  stress  to  yield  and  ultimate  are  similar,  with  7075-T6  having 
slightly  higher  values.  Therefore,  since  the  limit  design  stress  is  10  percent 
higher  for  the  7075  alloy  and  the  densities  are  comparable,  the  7075-T6  material 
would  be  chosen  over  the  2024-T81  alloy.  It  would  also  provide  higher  strength 
to  weight  capability  as  well  as  sufficient  fracture  toughness  to  cope  with  the 
inspectible  crack,* 

IX. 1(a)  Material  Substitutions 

If,  during  the  material  selection  phase,  or  perhaps  the  design  stage,  it 
may  be  desirable  to  consider  evaluating  different  materials,  how  can  this  be 
handled?  For  example,  it  may  be  planned  to  substitute  steel  sheet  in  place  of 
the  7075-T6  alloy  for  possible  weight  savings  and  higher  allowable  stresses. 
Mechanical  property  and  fracture  data  are  listed  in  Table  III  for  a  4340 
(0.035-inch  gage)  steel. 

To  determine  the  smallest  possible  crack  which  can  be  tolerated  in  the 
structure  at  the  limit  design  stress,  we  have,  when  Kj — Kc ,  2a/W — ►O, 

\ - ^1,00,  and  ^ 

a  =  — 
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and  for  critical  conditions, 


*  Environmental  effects  have  not  been  considered.  If  they  are  present,  7075-T6 
could  be  a  poor  choice  (see  Section  XI), 
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2 

-r  =  1.92  inches 
7 T 

for  the  7075-T6  alloy,  and  0.77  inches  for  the  4340  steel.* 

NOTE:  These  estimated  values  of  critical  crack  length  include  crack  tip  plas¬ 

tic  zone,  ry  length  due  to  the  use  of  plasticity  corrected  Kc  values; 
i.e.,  2ac  ^  2ac  by  visual  inspection,  but  2ac  =  2(ac  +  ry)  (see  Equa. 
VI-9,  Section  VI). 

With  limit  strength-density  ratios  about  equal  and  the  greater  crack 
tolerance  of  the  aluminum  alloy  at  the  design  stress,  the  7075  alloy  remains 
a  better  choice.  Other  operational  considerations  such  as  inspection  for 
small  cracks  and  stress-corrosion  cracking  lead  to  lower  threshold  Kc  values 
which  should,  for  comple teness ,  be  considered  in  this  substitution  study. 


“ 

2 

t - 

K 

85  ksi  */ inch 

c 

b  =  2 

7 T 

49  ksi  (1.00) 

It  has  been  determined  that  environmental  heating  of  the  structure  due  to 
performance  changes  in  the  prototype  design  makes  the  choice  in  aluminum 
marginal.  Titanium  alloy  Ti-6A1-4V  sheet  is  being  considered  as  a  replacement 
material.  Pertinent  material  and  fracture  data  for  the  Ti-bAl-4V  alloys  are 
shown  in  Table  III. 


TABLE  III 


DATA  FOR  MATERIAL  SUBSTITUTION  STUDY 


A1  loy 
Desig¬ 

cr 

uu 

cr  ( 1  ^ ) 

y 

Density 

P  ( lb/ 

K  (2, 3,4) 
c 

<rLiM.<5) 
Limit  Design 

p 

Critical 
Crack  Length 

nation 

(ksi) 

(ksi) 

cu . in.  ) 

(ksi^inch ) 

Stress  (ksi) 

(in  xlO^) 

2ac ( Inches ) 

7075-T6 

CLAD 

73 

64 

0.101 

85 

49 

485 

1.92 

4340 

200 

176 

0.283 

154 

140 

49  5 

0.77 

Ti -6A1 - 
4V 

134 

126 

0.160 

160 

89 

555 

2.06 

NOTES : 


(1)  MINIMUM  VALUES  FROM  MIL-HDBK-5. 

(2)  AVERAGE  TEST  VALUES  -  PARTICULAR  DESIGN  WOULD  CALL  FOR  FRACTURE 
TOUGHNESS  VALUES  FOR  THE  THICKNESS  AND  WIDTH  OF  INTEREST. 

(3)  PLASTIC  ZONE  ADJUSTMENTS  HAVE  BEEN  INCLUDED  IN  THESE  VALUES  AND 
NO  CRACK  BUCKLING  IS  ASSUMED. 

(4)  ENVIRONMENTAL  EFFECTS  HAVE  BEEN  NEGLECTED. 

(5)  BASED  ON  A  F.S.  OF  1.5. 


*  Although  the  enbr i t t lement  of  4340  alloy  at  the  high  <ru  is  such  that  this 
choice  would  be  for  comparison  purposes  only. 
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The  gross  fracture  stress  ( cr c)  for  the  1/2-inch  crack  in  the  titanium 
alloy  is  181  ksi.  A  calculated  critical  crack  length  at  the  maximum  limit 
design  stress  (89  ksi)  of  2.06  inches  can  be  tolerated  prior  to  fracture,  or 
a  10  percent  increase  in  critical  crack  length  over  the  7075-T6  alloy.  In 
addition,  the  strength-density  ratio  at  the  maximum  allowable  operating  stress 
is  14  percent  greater  for  the  titanium  alloy  than  for  either  7075-T6  or  4340. 
Therefore,  the  Ti-6A1-4V  would  be  a  good  fracture  resistance  choice  for  a 
substitute  material  for  this  application  (see  Note  below). 

IX . 1 ( b )  Other ^Methods  of  Material  Fracture  Comparison 

In  materials  trade-off  studies,  it  would  be  helpful  if  a  comparison  based 
on  weight  could  be  made.  For  example,  a  trade-off  between  weight,  yield 
strength,  and  fracture  toughness  (plane  strain),  first  presented  in  Reference 
25,  appears  as  shown  in  Figure  21(a).  Normalization  has  been  done  with  re¬ 
spect  to  a  materials  density.  An  upper  bound  curve  has  been  indicated,  setting 
the  upper  data  limit. 

On  a  linear  plot  of  Kjc  versus  cr yS  or  cru,  an  indication  can  be  obtained 
of  the  trade-offs  between  strength  and  toughness  (see,  for  example,  Reference 
45).  If  this  type  of  plot  is  density  normalized  (as  in  Figure  21  (a)  ),  a 
straight  line  drawn  from  the  origin  to  the  upper  bound  curve  at  a  particular 
strength  level  has  a  slope  of  Kxc/°"  ys  (see  Figure  21  (b)  ).  The  square  of 
this  slope  defines  a  test  specimen  of  suitable  dimension  for  fracture  tough¬ 
ness  determinations.  These  relationships  were  discussed  in  Section  VI.  It 
also  provides  an  indication  of  the  plastic  zone  size  (see,  e.g.,  Equation  VI-7). 

Any  valid  fracture  toughness  data  which  plot  below  the  upper  bound  curve 
indicate  that  the  material  is  not  as  crack  resistant  as  another  material  at 
the  same  strength  level.  In  this  manner,  the  density  normalized  plot  can  be 
used  to  compare  the  fracture  toughness  of  new  materials  with  other  materials 
at  a  given  strength  level. 

IX. 2  BENDING  LOADS 


An  arm  in  bending  is  shown  in  Figure  22.  It  has  developed  cracks  at  C, 
or  it  is  suspected  that  cracks  will  develop  at  C.  A  force  of  2,000  or  3,000 
lbs.  is  supplied  through  a  hydraulic  cylinder  at  point  B  and  the  bar  ends  are 
simply  supported. 

Question;  Which  material  would  be  a  better  fracture  resistant  choice  if 
a  1/4-inch  deep  throueh- the -thickness  crack  can  be  detected 
at  C;  7075-T6  or  2219-T851? 

This  geometry  has  significance  in  plane  strain  toughness  testing,  which 
will  be  discussed  in  Section  XII. 


NOTE:  Caution:  The  values  normally  reported  for  plane  stress  fracture  tough¬ 

ness  do  not  include  the  effects  of  environments  or  other  than  uniaxial 
loading.  Where  doubt  exists  as  to  the  component  stress  or  environment, 
simulated  tests  must  be  conducted  to  determine  or  fracture  stress 
values . 
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FIGURE  21  DENSITY  NORMALIZED  PLOT  OF  PLANE  STRAIN  FRACTURE  TOUGHNESS  AND 
YIELD  STRENGTH  FOR  DIFFERENT  MATERIALS  (REF.  25) 
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FIGURE  22  SIMPLY  SUPPORTED  BAR  OF  EXAMPLE  PROBLEM,  IN  BENDING 


STEP  1.  Equation  VII-8  gives  the  stress  intensity  factor  for  the  crack 

geometry  with  in-plane  bending.  Compute  the  stress  intensity  factor 
for  a  t-inch  deep  crack.  For  a  2,000  lb.  force, 

1  2  inc he<; 

M  =  2,000  lbs.  x  =  6,000  lb. -in. 


^  _  6M  JWa 

1  BW 


F  -  ^6,ppp.lb._-i.n,) - ^(0.25  inch)  =  1Q>6 

(0.75  inch) (2  inch)?  v 


for  a  3,000  lb.  force, 


M  =  3,000  lbs.  x  12  inches  =  9,000  lb. -in. 


KI  = 


_  6M  Jj a  _  6(9,000  lb. -in.)  ^77(0. 25  inch) 


BW2 


(0.75  inch) (2  inch) 2 


-  16.0  ksiy/inch 


STEP  2.  Compute  the  stress  intensity  factors  with  geometric  correction  using 
Figure  17  with  an  a/W  of  0.125,  A=  1.05  for  pure  bending.  For  a 
2,000  lb,  force  at  B 

K j  -  10.6  ks i^inch  (1.05) 

=  11.1  ks  flinch 

For  a  3,000  lb,  force  at  B 

=  16.0  ksiyinch  (1.05) 

-  16.8  ksiv^Lnch_ 


STEP  3,  Compare  these  values  of  stress  intensity  factor 
fracture  toughness  values. 


_  Kt 

MATERIAL  KIc  (ksi^nch)  (2,000  lbs.) 


7075-T6 


40 


0.28 


2219-T851 


32.6 


0.34 


to  the  plane  strain 
Kt 

(3,000  lbs.) 

KIc 

0.42 

0.52 


For  the  2,000  lb.  force,  either  material  would  be  suitable  for  this  design-- 
the  7075  alloy  being  a  slightly  better  choice.  The  2219  alloy  could  be 
considered  as  the  poorer  choice  for  the  3,000  lb.  force  and  1/4-inch  crack 
depth  condition,  if  environmental  considerations  are  not  a  factor. 


71 


IX. 3  OTHER  CRACK  GEOMETRIES 


In  service,  an  aircraft  can,  and  usually  does,  sustain  cracks  other  than 
the  center  crack  configuration.  There  are  cases  where  cracks  start  at  sheet 
edges,  on  the  surface,  or  at  angles  to  the  loading  path.  Fracture  mechanics 
can  handle  these  problems  with  reasonable  confidence  and  accuracy  for  areas 
where  plane  strain  mode  fracture  prevails  and  with  fair  success  for  the  mixed 
mode  and/or  plane  stress  situations  at  this  time. 

To  become  familiar  with  the  usage  of  the  fracture  mechanics  analysis, 
simple  crack  geometries  will  be  considered  as  they  may  occur  in  typical  de¬ 
signs.  In  this  manner,  one  can  become  familiar  with  the  steps  involved  and 
the  limitations  of  the  current  state-of-the-art  of  fracture  analysis. 

IX. 4  EFFECTIVE  OR  EQUIVALENT  CRACK  LENGTH 

In  cases  where  there  may  be  a  crack  traversing  from  the  edge  of  a  hole 
or  at  angles  to  the  loading  direction,  it  becomes  necessary  to  introduce  an 
equivalent  or  effective  "Griffith"  crack  length  rather  than  the  absolute  mea¬ 
sure  of  the  crack  length.  In  most  cases,  the  equivalent  crack  length  is  the 
same  as  the  absolute  crack  length  plus  the  discontinuity,  for  the  crack  and 
hole  situation.  However,  for  particular  applications,  estimates  should  be 
made  of  effective  length  magnitude,  and  compared  with  the  absolute  crack 
length.  The  requirements  for  calculating  equivalent  crack  length  will  be 
indicated  in  the  following  sections. 

IX . 4 ( a )  The  Crack  at  a  Hole 


A  design  for  a  proposed  feed-through  arm  is  shown  in  Figure  23. 

In  its  center  is  a  hole  with  a  tubular  insert  welded  at  the  top  and  bottom 
surface.  Fatigue  is  a  consideration,  particularly  if  a  brittle  weld  condi¬ 
tion  occurs.  Inspection  dictates  that  a  i-inch  crack  can  be  detected  at 
regular  inspection  periods  (we  will  assume  a  through- the-thickness  crack). 

The  arm  material  is  3/4  inches  thick  2219-T851  aluminum  wi th  a  yield  of 
51.2  ksi,  a  plane  strain*  fracture  toughness  of  32.6  ksi^inch,  and  an  operat¬ 
ing  stress  environment  of  35  ksi  prevails. 

Question:  Can  this  design  be  considered  safe  in  the  presence  of 

an  inspectable  crack,  at  a  257®  overload  condition  due 
to  pressure  pulses  at  P?  What  design  modifications 
would  be  recommended? 

STEP  1.  Assuming  the  worst  fatigue  crack  condition,  i.e.,  the  tube  has  split 
along  the  crack  axis,  it  can  be  seen  that  this  situation  then  becomes 
the  Bowie  crack  solution  (see  Section  VII  -  Eq.  VII-10)  as  indicated  at 
the  bottom  of  Figure  23. 

In  this  case,  to  solve  the  stress  intensity  equation,  an  equivalent  crack 
length  should  be  introduced.  In  other  words,  an  equivalent  "Griffith"  crack 
length  must  be  found  for  the  geometry  in  question. 

*  The  plane  strain  TKjc)  fracture  toughness  value  is  applicable  in  this  case 
with  3/4"  thick  plate  material.  It  would  be  a  safe  assumption  that  plane 
strain  fracture  behavior  predominates. 
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FIGURE  23  FEED  THROUGH  ARM  WITH  WELDED  SUPPORT  BASE 
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As  the  importance  of  the  effective  or  equivalent  crack  length  has  many 
implications  in  fracture  analysis,  we  will  present  this  typical  example  which 
will  have  direct  bearing  on  our  problem. 


The  stress  intensity  for  a  crack  or  cracks  at  a  hole  in  an  infinite 
plate  in  uniform  tension  is  from  Equation  VII- 10, 

Kj  -  <r  J(  f  ) 

and  the  stress  intensity  for  a  central  crack  in  an  infinite  plate  in  uniform 
tension  is  from  Equation  VII-1, 


K  =  <r  Jlfk 

N  e 


where  ae  is  the  equivalent  "Griffith"  crack.  By  equating  these  two  solutions, 
the  equivalent  crack  length  is  determined  by 


a 


=  a 


e 


2 


and  the  values 


are  shown  with  Equation  VII-10. 


At  this  point,  a  comparison  can  be  made  of  the  effective  crack  length 
and  the  actual  crack  length  of  the  Bowie  geometry.  For  a  1.0  inch  diameter 
hole,  the  equivalent  crack  lengths  are  shown  in  Figure  24.  For  the  half¬ 
crack  length  of  interest  (0.25  inch,  see  Figure  23)  ,  it  can  be  seen  that 
the  effective  crack  length  is  the  hole  radius  plus  actual  crack  length.  Be¬ 
yond  a  crack  length  of  0.12  inch,  there  is  a  little  difference  between  the 
equivalent  crack  lengths  computed  by  the  Bowie  solution  or  the  crack  length 
plus  radius  approximation. 


STEP  2.  The  equivalent  Grif f i th"crack  length  for  the  feed-through  arm  has 
been  found  to  be 


a  ^  a  +  r  =  0.25  inch  4-0.50  inch 
e 


or 


a  ^  0.75  inch 
e 

STEP  3.  Compute  the  finite  geometry  correction  factor  for  the  problem  using 
the  equivalent  crack  length.  The  crack  aspect  ratio  is  (from  Fig.  23) 

2a  _  1.5  inches  _  n  9c 
W  6  inches 

From  Figure  14,  X  =  1.04  for  a  crack  aspect  ratio  of  0.25. 
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10.0 


75 


STEP  4.  Solve  the  basic  stress  intensity  equation  for  critical  conditions; 

i.e.,  ae — ►aec,  cr - ►  crc,  and  K— *- Kxc  in  plane  strain.  Therefore, 

^Xc  32.6  ksijinch 

c  X  y 71(0 . 75  inch)  (1.04) 

o-  =  20.5  ksi 
c 


Clearly  this  design  will  not  tolerate  an  inspectable  crack  at  the 
anticipated  operating  conditions. 


What  would  be  the  least  tolerable  crack  for  this  design?  The  answer 
to  this  question  can  be  obtained  through  fracture  mechanics  by  assuming 
bounds  on  the  crack  length  parameter.  For  example,  manipulation  of  the  basic 
stress  intensity  equation,  again  assuming  critical  conditions  at  the  operat¬ 
ing  stress,  35  ksi, 


a 

e 


K 


Ic 


°cAj 


2 


gives  the  solution  for  equivalent  critical  crack  length.  However,  the  finite 
width  correction,  X  ,  is  also  a  function  of  crack  length.  For  this  crack 

geometry,  an  estimate  can  be  made  as  -fa - -0,  \  - -1.00,  so  the  smallest 

equivalent  critical  crack  length  will  be  0.276  inch.  Referring  to  Figure  24, 
an  equivalent  crack  length  of  0.276  inch  would  be  a  crack  less  than  0.05  inch 
long  from  the  edge  of  the  1-inch  diameter  hole.  This  perhaps  would  correspond 
to  a  slight  nick  caused  by  a  tool  at  the  hole  edge;  therefore,  a  re-evaluation 
of  the  design  is  inevitable. 

It  may  be  possible  to  introduce  four  reasonable  design  modifications 
which,  as  alternatives,  could  assure  a  fracture  safe  part.  They  are: 

1,  Reduce  allowable  operating  stress, 

2,  Reduce  hole  diameter, 

3,  Use  a  different  material  (higher  Kjc), 

4,  Reduce  plate  thickness  and  still  meet  stress  criteria. 

In  many  cases  it  is  difficult  to  change  the  configuration  to  satisfy 
new  operating  parameters.  Items  1  or  2  could  satisfy  the  design  requirement, 

°c  >  35  ksi,  for  the  inspectable  crack.  But,  we  have  already  seen  that  the 
smallest  tolerable  crack  for  this  geometry  is  very  small  (<0,05  inch),  which 
is  an  undue  restriction;  and  the  operating  stress  level  may  be  invariant. 

This  brings  us  to  the  third  or  fourth  alternatives.  Changing  to  a 
tougher  material  (higher  Kjc)  will  solve  the  problem.  However,  by  increasing 
the  material  plane  strain  fracture  toughness,  one  may  gain  some  of  the  necessary 
fracture  resistance,  but  not  enough  for  the  given  operating  stress  and  inspec¬ 
table  crack.  This  leaves  the  remaining  fourth  alternative  as  a  consideration. 


It  will  be  remembered  that  by  reducing  the  plate  thickness,  one  takes 
advantage  of  the  materials  ability  to  fracture  by  net  section  yielding.  In 
other  words,  a  brittle  fracture  behavior  had  been  assumed  and  by  reducing  the 
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thickness  a  mixed  mode  or  possibly  plane  stress  fracture  would  result.  This 
would  tend  to  result  in  higher  values  of  fracture  toughness.  Plane  stress 
(Kc)  fracture  toughness  values  of  %  113  ksi-s/inch^^ )  have  been  reported  for 
this  2219  material  in  0.10  inch  thick  sheet.  This  value  reflects  the  in¬ 
fluence  of  plastic  zone  size  (see  Eq.  VI-  2  ).  In  this  problem  we  could  safely 
assume  that  for  a  reduction  in  plate  thickness  to  3/8  inch,  we  would  gain 
%  507o  of  the  plane  stress  fracture  toughness,  or  56.5  ksi-Jinch.  (This  would 
be  a  conservative  estimate  and  would  imply  that  full  plastic  zone  development 
or  total  plane  stress  behavior  would  not  be  realized.) 

IX. 4(b)  Design  with  Reduced  Thickness 

STEP  5.  Compute  the  new  critical  fracture  stress  based  on  a  mixed  mode  type 

of  failure  with  a  reduced  plane  stress  fracture  toughness,  T^.^57  ksis/inch. 

*c  _  57  ksi^inch 

C  \Aaec  ^  y7T(0.75  inch)  (1.04) 

<r  =  36  ksi 
c 

Therefore,  this  re-design  (assuming  a  35  ksi  operating  stress)  would 
be  considered  fracture  safe  for  the  i;-inch  crack  condition,  but  not 
for  the  257,  overload  condition.  However,  compensation  must  be  made 
for  the  reduction  in  area  caused  by  the  reduction  in  thickness. 

The  designer  and  stress  analyst  must  now  decide  on  the  basis  of  stress 
analysis  of  this  re-design  if  a  reduction  in  plate  thickness  will  meet  static 
strength  requirements,  and  other  stress  criteria,  for  example.  This  is  where 
the  interplay  between  material/stress  and  fracture  parameters  enters  the  pic¬ 
ture,  and  trade-offs  are  selected. 


IX. 5  CRACK  AT  RIVETS 


Many  times  in  aircraft  design  a  part  may  develop  cracks  at  rivet  holes 
where  the  skin  is  attached  to  the  frame  or  stringer.  This  situation  is  de¬ 
picted  in  Figure  25,  and  will  be  analyzed  as  a  simple  case  in  which  the  sheet 
is  in  uniaxial  tension  and  the  rivets  above  and  below  the  crack  are  influen¬ 
tial  in  keeping  the  crack  closed.  (Tests  of  panels  with  concentrated  forces 
superimposed  on  uniform  tension  simulate  crack  growth  behavior  in  the  presence 
of  rivets.)  The  insert  of  Figure  25  shows  the  local  parameters  necessary  for 
determining  the  stress  intensities. 

Question:  Assuming  that  a  crack  grows  from  the  rivet  hole,  what 

is  the  local  stress  intensity  for  this  geometry? 

STEP  1.  This  problem  involves  the  use  of  stress  intensity  superposition  prin¬ 
ciples  which  were  discussed  earlier  in  this  report  (see  Section  V). 
Closer  examination  of  Figure  25  indicates  two  stress  intensities  as 
shown  in  Figure  26.  Geometry  (a)  and  associated  stress  intensity  is  for 
a  center  cracked  sheet  in  tension,  with  the  equivalent  crack  length 
approximately  equal  to  2a.  Geometry  (b)  and  associated  stress  intensity 
is  for  a  central  crack  subjected  to  concentrated  forces.  NOTE:  This 
stress  intensity  equation  is  applicable  only  if  the  concentrated  forces 
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/ 


FIGURE  25  CRACK  AT  RIVET  IN  A  RIVETED  SKIN-STRINGER  PANEL 
(NO  CRACK  BUCKLING) 
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UNIFORM  TENSION  CONCENTRATED  FORCE  SUPERIMPOSED 
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AND  CONCENTRATED  FORCE 


are  applied  at  the  center  line  of  the  sheet  and  at  a  distance  greater 
than  3  or  4  times  the  hole  diameter.  The  finite  width  corrections,  X, 
are  those  for  a  center  cracked  sheet  (see  Figure  14)  and  are  the  same 
in  both  cases. 

Inasmuch  as  the  concentrated  forces  are  in  an  opposite  direction  to 
the  uniform  stress,  and  tend  to  close  the  crack,  this  stress  intensity 
is  subtracted  from  the  uniform  extensional  stress  intensity  factor. 

STEP  2.  The  individual  stress  intensity  factors  can  now  be  superimposed  to 

compute  the  total  stress  intensity  factor  for  the  problem  geometry.  To 
better  understand  the  principles  of  superposition,  a  pictorial  repre¬ 
sentation  of  the  individual  and  superimposed  stress  intensities  is 
shown  in  Figure  27. 

With  knowledge  of  the  stress  intensity  solution  for  this  geometry, 
it  is  possible  to  determine  what  effect  the  rivet  closure  forces  have 
on  the  local  stress  field  for  similar  problems. 

Methods  exist  to  analyze  riveted,  cracked  structure  when  boundary  forces 
are  known.  These  will  be  demonstrated  later  in  this  section. 


IX. 6  ECCENTRICALLY  LOADED  CRACKS 


A  situation  which  occurs  quite  frequently  in  service  is  that  of  a 
crack  starting  at  a  rivet  hole  and  progressing  to  a  critical  length  by  fatigue 
crack  propagation.  Analyzing  this  type  of  problem  through  fracture  mechanics 
involves  both  equivalent  crack  length  and  combined  geometric  correction  prin¬ 
ciples. 


Take,  for  example,  the  geometry  of  Figure  28.  In  a  line  of  rivets, 
a  crack  has  started  at  rivet  A  and  has  progressed  to  a  length  of  3  inches  from 
the  hole  edge  prior  to  being  detected.  The  operating  stress  of  the  panel  is 
30  ksi,  and  it  is  fabricated  of  7075-T6  sheet  with  a  plane  stress  fracture 
toughness  of  85  ksiv/inch.  Overall  panel  width  is  2  feet.  (For  now  we  will 
not  consider  the  methods  available  to  analyze  riveted  structure,  but  only 
analyze  the  local  crack  problem  using  fracture  mechanics.)* 

Question:  What  are  the  critical  conditions  for  fracture  for  this 

crack  geometry  and  loading?  Should  the  panel  be  replaced? 

STEP  1.  In  examining  this  problem,  it  can  be  seen  that  it  is  the  crack  at 

hole  geometry  shown  as  an  insert  of  Figure  28.  (NOTE:  This  is  not  the 
simple  solution  for  a  centrally  stressed  crack  at  a  hole  as  in  Equation 
VII-10.)  This  problem  involves  an  eccentric  load  and  determination  of 
an  equivalent  crack  length  as  in  the  preceeding  problem. 

STEP  2.  Find  the  equivalent  Griffith  crack  length  and  compare  this  result  to 
the  hole  radius  plus  crack  length  approximation.  For  small  diameter 
rivet  holes,  the  function  $  (—) ,  from  the  corrections  of  Equation  VII- 
10,  for  a  single  crack  from  the  edge  of  a  hole,  is  0.707. 

*Two  opposing  point  rivet  forces  (above  and  below  the  crack)  are  assumed  in 
this  analysis. 
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STRESS  INTENSITY  FACTOR 


FIGURE  27  TRENDS  IN  STRESS  INTENSITY  „ FOR  LOADING  OF  FIGURE  26 
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FIGURE  28  EDGE  CRACKED  PANEL  GEOMETRY  FOR  EXAMPLE  EROBLEM 
(NO  CRACK  BUCKLING) 
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STEP 


STEP 


From  equivalent 
length) 


2a 


crack  length  computations  (hole  plus 
e  =  (a  +  r)  s  3.1875  inches 


absolute  crack 


for  a  3/16  inch  rivet  hole,  edge  crack  configuration.  This  is  compared 
to  3.0  inches  for  the  equivalent  crack  length  solution. 


There  is  a  stress  intensity  solution  available  for  eccentric  cracks 
in  tension^29)  which  includes  the  finite  width  correction.  The  stress 
intensity  factor  for  the  problem  is 

Ki  =  °V^7  ( x1>2) 


where  X  2  is  the  correction  factor  for  eccentric  loading  conditions 
and  a  function  of  the  eccentricity,  e,  and  X 2  is  the  finite  width 
correction  and  a  function  of  a/W1 .  (See  Figure  18). 


3.  Determine  the  crack  aspect  ratio  f  and  normalized  eccentricity 

(  6  )  for  the  problem. 


and 


£ 

W» 


1.5  inch 
(1.5  +  1  inch) 


0.60 


2e  _  2(11.6  inch) 
W  24  inch 


The  associated  correction  factor  for  eccentricity  and  finite  width  from 
Figure  18  is  X^,  X2  =  1.10  for  the  end  of  the  3  inch  long  crack  tip 
which  is  nearest  the  panel  center  line. 


4.  Assume  critical  conditions  for  the  stress  intensity  equation, 

Kt - -K  ,  <r - -cr  and  a  - -  a 

I  c*  c’  e  c 

and  solve  for  critical  stress. 

85  ksi^/Tnch" 

C  \/7lac  X L  2  v/71(1.5  inch)  (1.10) 

v  -  35.6  ksi 
c 

Hence,  it  can  be  assumed  the  panel  analyzed  as  a  simple  plate  would  be 
fracture  safe  under  the  30  ksi  operating  stress  with  this  length  crack. * 

The  second  part  of  the  question  may  be  re-phrased  for  our  purposes. 
When  should  the  panel  be  replaced?  Or,  in  particular,  at  what  crack 
length,  at  the  given  operating  stress,  should  the  panel  be  replaced? 

*  NOTE:  This  is  a  safe  assumption  as  long  as  the  other  end  of  the  crack 
remains  stopped  at  the  rivet  hole. 
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The  easiest  way  to  approach  this  problem  is  to  plot  the  representative 
stress-crack  length  history  of  this  eccentric  cracked  panel  as  shown  in  Figure 
29.  To  a  crack  length  of  4.0  inches,  this  panel  would  be  fracture  safe. 

This,  of  course,  depends  on  how  accurately  the  value  used  for  the  plane  stress 
fracture  toughness  reflected  the  material  behavior  in  service,  and  assumes 
that  the  crack  remains  stationary  at  the  rivet  hole.  NOTE:  There  are  addi¬ 
tional  ways  to  analyze  crack-riveted  structure  using  fracture  mechanics  prin¬ 
ciples.  These  will  be  presented  in  the  following  paragraphs. 

IX. 7  COMPLEX  THROUGH -CRACKED  STRUCTURE 


In  this  section  we  will  examine  the  analysis  of  built-up  structure 
(primarily  skin/stiffener  construction)  which  contain  a  crack  using  fracture 
mechanics  principles.  Skin-stringer  or  frame  construction  poses  a  particular 
fracture  analysis  problem  in  aircraft  design  when  a  crack  is  caused  by  sudden 
puncture  or  fatigue.  In  addition,  fail  safe  operation  specifies  load  carry¬ 
ing  requirements  for  particular  bay  sections  cracked.  To  meet  these  needs, 
it  is  possible  to  utilize  available  boundary  corrections  in  conjunction  with 
the  Griff ith-Irwin  fracture  toughness  parameters  to  compute  critical  fracture 
stresses . 


Examples  of  the  use  of  these  procedures  will  be  presented  as  a  series 
of  problems.  In  this  manner,  the  user  can  become  more  familiar  with  the  steps 
involved  in  the  analysis. 

NOTE:  Examples  presented  herein  are  for  illustrative  and  demonstrative 

purposes  only.  In  practice,  each  case  must  be  considered  as  indi¬ 
vidual  problems  with  unique  solutions. 


IX. 7 (a)  Damage  Configurations 

The  types  of  damaged  structure  which  can  be  analyzed  are: 

(1)  Elastic*  skin  with  central,  rivet  attached  stringer  and 
cracks  at  stringers, 

(2)  Elastic*  skin  with  fixed  edges  and/or  cracks  between  stringers, 

(3)  Combinations  of  1  and  2. 

Figures  30  and  31  show  the  various  combinations  of  damaged  structure 
for  which  the  analysis  is  valid  with  applicable  nomenclature. 

IX. 7(b)  Geometric  Corrections 


In  1964,  J.  Bloom  presented  an  analysis  of  the  boundary  corrections 
for  a  cracked,  riveted  sheet  stringer(^)  construction  under  tensile  loads 
(see  Figures  30(a)  and  (b)  and  31).  The  following  year  Grief  and  Sanders 


*  Elastic  behavior  is  assumed  in  the  derivation  of  all  the  analysis  methods, 
even  though  elas to-plas tic  action  is  known  to  occur  at  the  crack  tip.  If 
these  zones  are  small  compared  to  the  crack  length,  then  these  assumptions 
are  valid  and  fracture  mechanics  principles  can  be  used  to  provide  reason¬ 
able  results. 
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CRITICAL  FRACTURE  STRESS,  cr  (ksi) 


rr 


FIGURE  29  CRITICAL  FRACTURE  STRESS  FOR  PANEL  GEOMETRY 
OF  FIGURE  28  (NO  CRACK  BUCKLING) 
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(a) 


(b) 


CRACK  BETWEEN  RIVET 


CRACK  THROUGH  RIVET  HOLE 


CT 


I  t  t 


MULTIPLE  CRACKS 
(REINFORCED  EDGES) 


<r 


(d) 

CENTER  CRACK 
(REINFORCED  EDGES) 


FIGURE  30  SYMMETRIC  CRACK  CONFIGURATION  IN  UNIAXIAL  TENSION 
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cr 


FIGURE  31  NONSYMMETRIC  CRACK  IN  RIVETED  PANEL  -  UNIAXIAL  LOADING 
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presented  the  boundary  corrections  for  a  continuous  stringer-cracked  sheet  combi¬ 
nation^  .  Isida,  et  al,  (Reference  28)  have  calculated  the  crack  tip  stress 
intensities  for  various  combinations  of  extensional  and  bending  stiffness  of 
the  stiffeners  for  the  crack  geometry  of  Figures  30(c)  and  30(d). 

The  associated  correction  factors  for  the  crack  geometries  of  Figures 
30  and  31  are  shown  graphically  in  Figures  32  through  40.  The  crack  between 
rivets  (stringer  intact)  is  indicated  in  Figure  32.  For  the  crack  through  the 
rivet  hole  and  stringer  intact,  see  Figure  33. 

The  boundary  corrections  for  a  nonsymmetr ic  crack  near  a  riveted 
stringer  for  various  crack  aspect  ratios,  for  values  of  extensional  stiffness, 

S,  equal  to  0.5,  1  and  2  are  shown  in  Figures  34  and  35.  Boundary  corrections 
for  a  center  cracked  plate  in  tension,  with  reinforced  edges,  as  analyzed  by 
Isida  and  Itagahi^S)  are  shown  in  Figures  36  through  40  for  various  values 
of  extensional  stiffness,  S,  and  bending  stiffness  of  the  stringers,  CL  .  The 
CL  term  is  essentially  a  nondimens ional  inertia  parameter.  It  is  a  function 
of  the  moment  of  inertia  of  the  stiffener  about  its  neutral  axis  perpendicular 
to  the  plane  of  the  plate  or  sheet,  Iz,  divided  by  the  product  of  the  plate 
or  sheet  thickness,  B,  and  the  cube  of  the  half  stringer  spacing  or  sheet  width 
W 

— .  This  factor  is  multiplied  by  the  nondiraens ional  ratio  of  stringer  to  plate 
modulus . 


CL 


I 

z 


(IX-1) 


The  dashed  portion  of  the  correction  factor  curves,  in  all  cases, 
reflects  the  limiting  values  obtained  experimentally  by  Grief  and  Sanders (^) 

on  integrally  stiffened  panels  assuming  a  value  of  —  of  100. 

P 


IX. 8  GEOMETRIC  COMBINATIONS 


To  analyze  the  panel  section  of  Figure  41  for  fracture  strength,  we 
will  utilize  applicable  geometric  corrections  and  fracture  mechanics.  The 
following  conditions  apply  to  the  panel  of  Figure  41, 

.  center  cracked-perpendicular  to  loading  direction 
.  cracked  at  rivet  with  stringer  intact 

.  uniaxial  loading  (sheet  and  stringer  stressed  uniformly, 

remotely  from  the  crack) 

.  rivet  pitch  -  1.00  inch 

.  symmetrical  crack  growth  with  respect  to  panel  center  line 
.  repeative  structure  (analyze  only  portion  of  structure-3  bays). 

To  determine  the  influence  of  the  central  stringer  on  the  critical  fracture 
stress  for  this  geometry,  we  will  separate  the  problem  into  parts.  In  this 
way,  the  influence  of  the  center  stringer  will  become  more  apparent. 

Consider  the  problem  first  as  a  center  cracked,  edge  stiffened  sheet 
with  center  stiffener  removed.  The  geometry  of  Figure  30(d)  applies.  In 
Reference  28,  the  edge  stiffeners  are  considered  as  integral;  however,  it  is 
felt  that  negligible  error  is  introduced  for  crack  aspect  ratios  .<0.8. (28) 

As  the  crack  approaches  the  edges  stiffeners,  the  influence  of  these  stiffeners 
predominates  and  the  crack  geometry  of  Figure  31  applies. 
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FIGURE  32  CORRECTION  FACTORS  FOR  SYMMETRIC  CRACK  MIDWAY  BETWEEN  RIVETS 
RIVETED  STRINGER  -  STRINGER  INTACT  (REF.  47) 


half  crack  length 
rivet  pitch 
plate  thickness 
modulus  of  plate 
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FIGURE  33  CORRECTION  FACTORS  FOR  SYMMETRIC  CRACK  THROUGH  RIVET  HOLE 
RIVETED  STRINGER  -  STRINGER  INTACT  (REF.  47) 


CORRECTION  FACTOR, A  CORRECTION  FACTOR 


a  =  half  crack  length 
p  =  rivet  pitch 

B  =  plate  thickness 
E  =  modulus  of  plate 
EgiR  =  modulus  of  stringer 


STRINGER  £  ASTR  area  of  Wringer 


STRINGER  £ 


0.8 


0.4 


CRACK  ASPECT  RATIO, 


2a 

W 


FIGURE  34  CORRECTION  FACTORS  FOR  NONSYMMETR1C  CONDITIONS  -  CENTER 
CRACKED  PANEL  FOR  S'  =  0.5  AND  1.0  (REF.  47) 
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CORRECTION  FACTOR 


a  -  half  crack  length 
p  =  rivet  pitch 

B  =  plate  thickness 
w  E  =  modulus  of  plate 


FIGURE  35  CORRECTION  FACTORS  FOR  NONSYMMETRIC  CONDITIONS  -  CENTER  CRACKED 
PANEL  FOR  S'  =  2.0  (REF.  47) 
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CORRECTION  FACTOR 


W 

FIGURE  36  CORRECTION  FACTORS  FOR  EDGE  STIFFENED,  CENTER  CRACKED  PLATE  - 
BENDING  STIFFNESS,  Ot  =  0  (REF.  28) 
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EXTENSIONAL  STIFFNESS, 


CORRECTION  FACTOR 
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FIGURE  37  CORRECTION  FACTORS  FOR  EDGE  STIFFENED,  CENTER  CRACKED 
PLATE  -  BENDING  STIFFNESS,  0!  =  0.01  (REF.  28) 
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EXTENSIONAL  STIFFNESS 


CORRECTION  FACTOR 
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FIGURE  38  CORRECTION  FACTORS  FOR  EDGE  STIFFENED,  CENTER  CRACKED  PLATE  - 
BENDING  STIFFNESS,  a  =  0.03  (REF.  28) 
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EXTENSIONAL  STIFFNESS, 


CORRECTION  FACTOR 


1 .6"T 


FIGURE  39  CORRECTION  FACTORS  FOR  EDGE  STIFFENED,  CENTER  CRACKED  PLATE 
BENDING  STIFFNESS,  OC  =  0.05  (REF.  28) 
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FIGURE  40  CORRECTION  FACTORS  FOR  EDGE  STIFFENED,  CENTER  CRACKED  PLATE  - 
BENDING  STIFFNESS,  a  =  0.09  (REF.  28) 
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EXTENS10NAL  STIFFNESS 


(no  crack  buckling) 


FIGURE  41  CENTER  CRACKED  PANEL  GEOMETRY  OF  EXAMPLE  PROBLEM 
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STEP  1.  Determine  the  nondimensional  inertia  parameter,  from  Equation  IX-1. 
In  this  case  (Iz)  the  moment  of  inertia  of  the  edge  stiffener  (angle 
section  length  2  inch,  thickness  1/4  inch)  about  the  stringer  neutral 
axis-perpendicular  to  the  plane  of  the  plate  is 

I  =  0.322  inches4 
z 

and  the  ratio  of  stringer  to  sheet  modulus  is  1.0.  Therefore,  a  is 


a  = 


0.322  inches 


0.063  inch  ( 


18  inches 


r 


10 


psi 


10  psi 


S  0.01 


and  the  correction  factors  of  Figure  37  apply. 


STEP  2.  Compute  the 
strip.  This  is 

S 


extensional  stiffness, 
from  Figure  37. 


S,  for  the  finite  width,  cracked 


(IX- 2) 


c  -  inch  (2  inch  +  1.75  inch)  10^  psi 

.  / 18  inches  \  lrt7  .  # 

0.063  inch  ( - ^ - f  10  psi 

STEP  3.  Prepare  a  table  of  correction  factors  (X*)  for  various  crack  aspect 
ratios,  as  shown  in  Table  IV  for  S  =  1.65  from  Figure  37. 


TABLE  IV 

CORRECTION  FACTORS  FOR  EXAMPLE  PROBLEM 


TOTAL 

CRACK  LENGTH, 

2a  (inches) 

CRACK  ASPECT 

RATIO, 

2a/W 

A' 

(From  Fig.  37) 

1.8 

0.1 

1.1 

00 

3.6 

0.2 

5.4 

0.3 

7.2 

0.4 

9.0 

0.5 

1.00 

10.8 

0.6 

0.' 

98 

12.6 

0.7 

0.' 

97 

14.4 

0.8 

0.' 

91 

As  the  crack  approaches  the  edge  stiffeners  —  >  0.8  the  analysis  of  Reference 
28  is  not  valid.  However,  the  problem  becomes  one  of  a  non-symmetr ic  crack 
approaching  a  riveted  stiffener,  or  the  geometry  of  Figure  31.  (With  the 
center  stiffener  broken,  little  effect  will  be  noticed  on  the  correction  from 
this  stringer. ) 
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STEP  4.  Determine  the  correction  factors  as  the  crack  goes  through  the  edge 

stringer  areas  and  beyond  (see  Figure  34  and  35).  Compute  the  stiffness 
ratio  S'  and  ratio  of  crack  length  to  rivet  pitch  for  crack  aspect 
ratios  >0.8.  Where  Sf  is 


2  (aB)  E 

'  - F— 

str.  str. 

and  corresponding  correction  factors  are  shown  in  Table  V. 

TABLE  V 


CORRECTION  FACTORS  FOR  EXAMPLE  PROBLEM 


2a 

W 

S’ 

(?)* 

X"** 

0.9 

1.09 

4.05 

0.70 

1.0 

1.21 

4.50 

0.45 

1.1 

1.33 

4.95 

0.28 

1.2 

1.45 

5.40 

0.32 

*  The  crack  is  assumed  to  grow  through  a  rivet 
hole,  thus  effectively  doubling  the  rivet 
pitch,  p'  =  2p. 

**From  Figures  34  and  35 


(IX-3) 


The  correction  factors  X1  and  XM  for  the  two  crack  geometries  are 
shown  in  Figure  42  as  a  function  of  crack  aspect  ratio.  It  is  apparent 
that  the  maximum  reduction  in  correction  factor  will  occur  when  the 
crack  tip  is  approaching  the  riveted  stiffener  (0.9  >">  0.7).  The 
correction  factor  curves  for  a  nonsymmetric  crack  (Figures  34  and  35) 
also  indicate  this  trend  for  small  to  moderate  rivet  spacings.  This 
reduction  is  caused  by  the  remaining  rivet  forces,  near  the  crack  tip 
tending  to  pinch  the  crack  together  more  effectively  as  the  crack  is 
within  the  riveted  stringer  area.  Before  the  crack  reaches  this  area, 
this  effect  is  lost. 


STEP  5.  We  can  now  examine  the  effect  of  the  center  stiffener  on  the  correc¬ 
tion  factor  to  the  stress  intensity.  As  the  crack  is  assumed  to  start 
in  the  sheet  at  the  center  stiffener  and  progress  in  a  symmetric  manner, 
the  corrections  of  Figures  32  or  33  apply  depending  on  where  the  crack 
starts,  i.e.,  at  or  between  rivets.  The  author  of  Ref.  47  has  computed 
the  correction  curves  of  Figures  32  and  33  using  a  Poisson’s  ratio  of 
0.33  and  a  rivet  diameter  to  rivet  pitch  ratio  of  -jr  =  0.1.  However, 
he  did  indicate  relative  insensitivity  of  the  correction  factor  to 
variation  in  this  ratio.  It  should  be  noted  that 
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FIGURE  42  CORRECTION  FACTORS  FOR  SYMMETRIC  -  CRACKED  PLATE  WITH  REINFORCED 
EDGES  FOR  THE  GEOMETRY  OF  FIGURE  41 
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STRINGER  £ 


the  crack  running  through  the  rivet  hole  effectively  doubLes  the  pitch 
over  that  of  a  crack  midway  between  rivets  (i.e.,  compare  Figures  32 
and  33),  We  will  first  examine  the  most  plausible  situation,  a  crack 
starting  at  a  rivet  with  stringer  intact.  In  this  case,  the  geometric 
corrections  are  shown  in  Table  VI  from  Figure  33. 

TABLE  VI 

CORRECTION  FACTORS  FOR  EXAMPLE  PROBLEM 


2a 

W 

s 

a 

P 

X"  '  (**) 

0.1 

0.12 

0.90 

0.73* 

0.2 

0.24 

1.80 

0.68* 

0.3 

0.36 

2.70 

0.68* 

0.4 

0.48 

3.60 

0.675 

0.5 

0.60 

4.50 

0.68 

0.6 

0.73 

5.40 

0.68 

0.7 

0.85 

6.30 

0.69 

0.8 

0.97 

7.20 

0.69 

0.9 

1.09 

8.10 

0.69 

1.0 

1.21 

9.00 

0.695 

1.1 

1.33 

9.90 

0.695 

1.2 

1.45 

10.80 

0.70 

*  Obtained  by  Interpolation 
**From  Figure  33 

STEP  6.  Multiply  the  correction  factors  obtained  in  Table  IV  for  an  edge 

stiffened  plate  with  those  above  to  obtain  a  combined  correction  factor 
A*  and  A**«  These  factors  are  defined  and  indicated  in  Table  VII 
along  with  the  nonsymmetric  case  from  Table  V. 

TABLE  VII 

COMPOSITE  CORRECTION  FACTORS  FROM  TABLES  IV-VI 


2a 

W 

X*  =  X'  x  A"' 

X**  =  X"  x  X"' 

0.1 

0.73 

0.2 

0.68 

0.3 

0.68 

0.4 

0.68 

0.5 

0.68 

0.6 

0.67 

0.7 

0.66 

0.8 

0.63 

0.9 

0748 

1.0 

0.31 

1.1 

0.19 

1.2 

— 

0.22 
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IX. 8 (a)  Selection  of  Fracture  Toughness  Value 


With  the  use  of  fracture  mechanics  we  can  now  make  an  estimate  of 
critical  fracture  stress  for  each  crack  aspect  ratio.  Plane  stress,  Kc 
data***  on  0.063  inch  thick  7075-T6  aluminum  sheet  indicated  a  ksi>/inch  . 

This  average  plane  stress  fracture  toughness  value  includes  plastic  zone 
correction.  However,  when  any  correction  for  plasticity  is  employed  it  will 
lead  to  larger  values  of  fracture  toughness,  hence  fracture  stress.  As  a 
matter  of  choice,  it  is  usually  better  practice  to  base  preliminary  design 
and  fracture  analysis  on  a  fracture  toughness  value  which  does  not  include 
plasticity  corrections.  The  reason  is  that  any  calculations  of  fracture 
stress  will  tend  to  be  conservative  particularly  for  the  small  crack  situation 
(see  e.g.  Eq.  VI-5,  Section  VI).  For  the  tougher  materials  (e.g.  7075-T6  and 
2024-T3  aluminum)  this  may  become  unduly  conservative  and  lead  to  restrictive 
designs.  In  these  cases  it  would  be  better  to  include  plasticity  corrections 
in  the  analysis.  The  examples  presented  herein  assume  plasticity  corrections 
to  the  plane  stress  or  plane  strain  fracture  toughness. 

STEP  7.  Using  equation  VIII-1  compute  the  critical  fracture  stress  for  the 
crack  aspects  ratios  of  interest.  Equation  VIII-1  becomes,  for  this 
example 


K 


c 


cr 


(IX-4) 


for  critical  conditions  (i.e.  K — -  Kc ,  a - -ac,or  ^°^).  The  trend  in 

fracture  stress  with  increasing  crack  aspect  ratio  as  well  as  the  compo¬ 
site  correction  factor  is  shown  in  Figure  43. 

Thus,  for  the  geometry  of  Figure  41  the  fracture  stress  decreases  until 
the  crack  approaches  the  edge  stringers  (14. 5"  >  2a  >  0) .  The  edge  stringers 
start  to  pick  up  the  load  at  that  point  and  the  maximum  fracture  stress  is 
reached  once  the  crack  goes  through  the  stringer  area  (2a  =  20  inches).  It  is 
interesting  to  note  that  the  influence  of  the  center  stringer  on  the  correction 
factor  remains  fairly  constant  throughout  the  crack  range  of  interest.  (See 
Table  VI.) 

If  this  stringer  is  broken,  however,  the  effect  on  critical  stress  is 
quite  pronounced.  In  this  case  the  crack  geometry  becomes  one  of  Figure  30(d) 
until  a  crack  aspect  ratio  of  0.8  is  reached  and  the  geometry  of  Figure  31 
applies.  Critical  fracture  stresses  can  be  computed  as  before  using  the 
geometric  corrections  of  Table  IV  and  V.  Fracture  stresses  are  plotted  in 
Figure  44  for  the  broken  center  line  stringer  case.  Comparison  of  this 
analysis  with  that  for  the  stringer  intact  (dashed  curve)  indicate  an  approxi¬ 
mate  307o  reduction  in  fracture  stress  for  most  values  of/2ajfor  this  condition. 
With  the  center  line  stringer  broken,  there  will  be  some'^crack  buckling 
resulting  in  lower  fracture  stresses  at  the  larger  crack  lengths.  The 
analysis,  as  presented,  does  not  account  for  this  effect.  However,  it  is 
felt  that  this  reduction  will  be  minimal  due  to  the  influence  of  the 


***Valid  fracture  data  from  the  material  in  question  must  be  available  for  the 
thickness  of  interest  to  solve  this  problem.  This  implies  "infinite"  sheet 
data. 
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FIGURE  43  COMPOSITE  CORRECTION  FACTORS  &  CRITICAL  STRESS  AS  A 

FUNCTION  OF  CRACK  LENGTH  FOR  GEOMETRY  OF  FIGURE  41  WITH 
£  STRINGER  INTACT  (NO  CRACK  BUCKLING) 


Oi 


CRACK  ASPECT  RATIO,  % 


I  I  1  I  I  II 

10.8  9.0  4.5  o  4.5  9.0  10.8 

CRACK  LENGTH,  a  (INCHES) 

FIGURE  44  COMPARISON  OF  PREDICTED  FRACTURE  STRESS  FOR  CRACK  GEOMETRY  OF 

FIGURE  41  -  CENTER  STRINGER  INTACT  &  BROKEN  (NO  CRACK  BUCKLING) 
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2 

edge  stiffeners  within  the  important  crack  aspect  ratios  of  0.8>-^->0.2. 

For  cases  of  repetative  structure,  it  can  be  seen  that  this  problem 
then  becomes  similar  to  the  geometry  (many  skin/stringer  bays)  of  Figure 
30(c)  and  can  be  treated  in  the  same  manner.* 

IX. 9  NON-SYMMETRIC  CRACK- INTEGRAL  STRINGER 

In  1965,  Grief  and  Sanders^®)  performed  an  analysis  of  a  cracked  plate 
with  continuously  attached  (integral),  zero  width  stringer.  They  assumed  zero 
bending  stiffness  and  continuous  attachment  (integral)  of  the  stringer.  There¬ 
fore,  an  analysis  of  an  edge-stiffened  sheet  can  be  performed  for  the  geometry 
of  Figure  30(d)  using  the  Isida  corrections  for  crack  aspect  ratios  ^  <  0.8 
and  the  Grief -Sanders  results  for  >0.8.  Comparisons  can  then  be  made 
between  critical  stress  for  an  integral  stiffened  and  riveted  stringer.** 

IX. 9(a)  Example 

In  this  case,  we  have  the  structure  as  shown  in  Figure  45.  It  is  similar 
to  that  of  Figure  41  with  centerline  stringer  removed  and  0.125,!  thick  hat 
section  stringers  installed  instead  of  angles.  The  hat  sections  are  welded 
to  the  skin  or  could  be  riveted  with  very  small  pitch. 

STEP  1.  Compute  the  correction  factors  for  crack  aspect  ratios,  ^  <  0.8. 

In  this  case  from  equation  IX-1  and  Iz  =  0.385  inches^, 

IZ  . Estr.  _  0.385  inch4  ,  lO^psl 

®  ~  _/W  \  3  E  _ _ _  .  / 18  inches  \ 3  ,„7 

B^-)  0.063  inch( - ^ - )  10  psi 


fit  a  0.01 

As  in  the  previous  problem,  the  correction  curves  of  Figure  37 
apply.  The  extensional  stiffness  is  (Equation  IX-2), 

g  _  Astr .  _  Estr  . _ 0.594  inch^ _  lC^psi 

b(|)  E  0.063  inch ( 18  ^ph"8- )  lC7psi 

S  =  1.05 


and  the  correction  factors  are  shown  in  Table  VIII. 


*This  problem  has  recently  been  analized  by  C.  Poe,  Jr.,  ’’The  Effect  of 
Riveted  and  Uniformly  Spaced  Stringerson  the  Stress  Intensity  Factor  of 
a  Cracked  Sheet,”  M.S.  Thesis,  Va.  Polytechnic,  1969. 

**N0TE  Based  on  the  assumptions  in  this  and  previous  analyses,  satisfactory 

results  will  be  obtained  as  long  as  the  distance  of  the  crack  from  the 
stringer  is  on  the  order  of  the  rivet  spacing  or  stringer  width. 


FIGURE  45  SAMPLE  PROBLEM  PANEL  GEOMETRY 
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TABLE  VIII 

CORRECTION  FACTORS  FOR  EXAMPLE  PROBLEM 


2a 

W 

X  x 
(Fig.  37) 

0 

1.0 

0.1 

1.0 

0.2 

1.01 

0.3 

1.015 

0.4 

1.02 

0.5 

1.025 

0.6 

1.025 

0.7 

1.015 

0.8 

0.98 

STEP  2.  As  the  crack  approaches  the  hat  stiffeners,  the  correction  factors 
of  Figure  46  apply  and  the  values  for  S  and  X'  are  given  in  Table  IX. 
S  in  this  case  is 


S 


4EB  / 

. —  '  > 

Astr  Estr  (1+V)(3-V) 


where 


I  2a  +  W 
1  =  2 


for  a  center  crack  geometry. 


(IX -5) 


TABLE  IX 

CORRECTION  FACTORS  FOR  EXAMPLE  PROBLEM 


2a 

W 

l,  (inch) 

S 

X’ 

(Fig.  46) 

0.6 

14.4 

1.72 

0.925 

0.7 

15.3 

1.83 

0.885 

0.8 

16.2 

1.94 

0.810 

0.85* 

16.65 

1.99 

0.785 

0.90* 

17.1 

2.04 

0.690 

★Analysis  not  valid  for  i-2a  >  stringer  width. 


It  will  be  noticed  that  there  is  a  small  ( <  0.2)  difference  in  correction 
factors  for  crack  aspect  ratios,  between  0.6  and  0.8  in  Tables  VIII  and 

IX.  Therefore,  average  composite  correction  factors,  \  *,  between  Tables  VIII 
and  IX  are  indicated  in  Table  X  for  crack  aspect  ratios  (-^jof  0.6,  0.7,  and 
0.8  for  the  geometry  of  Figure  45. 
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FIGURE  46  CORRECTION  FACTORS  -  NONSYMMETRIC  CRACK  -  INTEGRAL  STRINGER  (REF.  48) 


TABLE  X 

COMPOSITE  CORRECTION  FACTORS  FROM  TABLES  VIII  &  IX 


2a 

W 

X 

x' 

* 

X 

0 

1.0 

- ■ 

1.0 

0.1 

1.0 

— 

1.0 

0.2 

1.01 

— 

1.01 

0.3 

1.015 

— 

1.015 

0.4 

1.02 

— 

1.02 

0.5 

1.025 

— 

1.025 

0.6 

1.025 

0.925 

0.975 

0.7 

1.015 

o.sas 

0.950 

0.8 

0.98 

0.810 

0.895 

0.85 

— 

0.785 

0.785 

0.90 

— 

0.690 

0.690 

A  plot  of  critical  stress  for  the  integral  stringer  of  this,  example  com¬ 
pared  with  the  riveted  stringer  of  the  previous  problem  is  shown  in  Figure  47, 
assuming  a  plane  stress  fracture  toughness,  Kc  =  85  ksis/inch. 

It  will  be  noticed  there  is  little  difference  between  the  riveted  edge 
stringer  and  welded  stringer  results.  In  the  case  of  the  riveted  stringer, 
the  rivet  pitch  is  small  enough  that  it  can  be  assumed  that  it  will  behave  in 
a  manner  similar  to  the  integral  stringer.  This  is  borne  out  by  the  results 
plotted  in  Figure  47.  Also  indicated  in  this  comparison  is  the  trend  in 
stress  pick-up  by  the  edge  stringer  at  a  crack  aspect  ratio  of  0.8,  which 
corresponds  to  the  riveted  structure  case. 


IX. 10  SUMMARY 

Thus  the  fracture  strength  can  be  computed  for  riveted  skin/stringer 
and  integral  stringer/skin  construction  using  fracture  mechanics  as  a  base. 

It  should  be  remembered  that  the  analysis  as  presented  here  has  included 
crack  tip  plasticity  which  can  be  quite  large  for  tougher  materials.  Neg¬ 
lecting  this  plasticity  can  lead  to  overly  conservative  results  for  the  higher 
fracture- tough  alloys.  Therefore,  this  type  of  analysis  is  limited  only  by 
the  selection  of  the  proper  value  of  plane  stress  fracture  toughness  (Kc). 

When  possible,  this  value  should  be  obtained  by  test,  from  cracked  sheets 
(preferably  of  the  geometry  of  interest)  of  the  thickness  which  is  under  con¬ 
sideration.  However,  lacking  positive  fracture  toughness  data  for  the  mate¬ 
rial  of  interest,  a  conservative  estimate  of  fracture  stress  can  be  computed. 
In  the  following  section,  it  will  be  shown  how  fracture  envelopes  can  be  use¬ 
ful  in  this  regard. 
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IX. 11  FRACTURE  STRESS  ESTIMATES  WITH  AVERAGE  TOUGHNESS  DATA 


In  many  cases  an  estimate  of  critical  stress  or  tolerable  crack  length 
must  be  made  for  designs  for  which  the  material's  fracture  data  is  known  from 
only  a  limited  number  of  tests.  Perhaps  average  Kc  data  is  not  available  for 
the  thickness  of  interest.  Or  dissimilar  materials  may  be  under  considera¬ 
tion  for  skin/stiffener  applications  in  a  particular  design.  Through  the  use 
of  the  composite  correction  principles  shown  thus  far,  it  is  possible  to  make 
preliminary  estimates  of  critical  fracture  stress  based  on  fracture  mechanics 
analysis.  One  such  problem  will  be  analyzed  to  indicate  the  applicability  of 
these  methods. 


A  design  shown  in  Figure  48  could  represent  a  portion  of  a  lower  center 
wing  panel  which  is  subjected  to  primary  tension  loads.  Width  between  the 
ZEE  stringers  is  9  inches.  In  this  problem,  consideration  is  being  given  to 
two  skin  materials,  2219-T87  and  Ti-6A1-4V.  If  the  titanium  alloy  is  used, 
the  ZEE  stringers  will  be  of  a  high  strength  steel  alloy.  For  the  aluminum 
construction,  the  ZEE  stringers  will  also  be  aluminum.  Average  plane  stress 
fracture  toughness  values  (Kc)  for  the  T1-6A1-4V  skin  range  from  155-165  ksi 
s/inch . *  For  the  2219-T87  aluminum  alloy,  Kc  ranges  from  115-125  ksis/inch.** 

To  provide  the  maximum  data,  two  adjacent  bays  should  be  analyzed,  assum¬ 
ing  a  crack  would  start  at  a  rivet  hole  in  the  sheet  at  the  center  ZEE  strin¬ 
ger.  It  can  be  seen  from  Figure  48  that  the  wing  panel  is,  for  the  crack 
problem,  similar  to  the  geometry  of  Figure  41.  Therefore,  the  crack  geometry 
is  the  composite  of  Figures  30  (b  and  d)  and  31,  and  the  procedures  followed 
in  the  first  reinforced  panel  problem  of  this  section  are  applicable. 

IX. 11 (a)  Computation  of  Geometric  Corrections 

STEP  1.  The  moment  of  inertia  of  the  ZEE  stringer  about  the  stringer  neutral 
axis,  perpendicular  to  the  plane  of  the  plate,  Iz  is  0.0511  inches^,  and 
the  Modulus  of  the  steel  stringer  is  28  x  10^  psi,  and  the  titanium 
sheet  16  x  10^  psi. 


For  the  all  2219-T87  aluminum  alloy  construction,  the  inertia  para¬ 
meter  (Eq.  IX-1)  is  calculated  as 

a  .  . , „  = - — inch  r - — P'—  =  o.oooi  2  o 


A //A/ 


0.100  inch  (• 


18  inches  \  3 


lO^ps i 


Assuming  the  same  ZEE  section  properties  for  the  steel  stringers  (i.e 
I7  =  0.0511  inches^)  for  simplicity,  the  titanium  panel  with  alloy 
steel  stringers  gives  a  nondimensional  inertia  parameter  of 


Oi  / 
Ti/St 


0,0 


_ 0.051 1  inch^ 

,  / 1 8  inches 
q0  inch^ - - - 


2.8  x  lO^psi  a  0.002  ^  0 
1.6  x  10  psi 


Therefore,  the  corrections  of  Figure  36  apply  in  this  analysis  to  both 
skin  materials. 


The  extensional  stiffnesses  are  from  Figure  36,  again  assuming  equal 
ZEE  stringer  areas. 

*  Approximate  range  in  values  for  B  =  0.100  inch.  Wide  panel  tests. 

**  Approximate  range  in  values  for  B  =  0.05  inch.  Wide  panel  tests. 


112 


ZEE  STRINGER: 

LEG  WIDTH  =  1.00  INCH 
HEIGHT  ~  1.50  INCHES 
AREA  -=  0.314  INCH2 
RIVET  PITCH  1.0  INCH 


FIGURE  48  LOWER  CENTER  WING  PANEL  FOR  EXAMPLE  PROBLEM 
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2  7 

-  =  0.314  inches _  .  10  psi  ^  n  .c 

A^/A^  0.100  inch  ( 18  1^ches)  I07psi 

and 

2  7 

_ 0,314  inches _  #  2,8  x  10  psi  ^  .  , 

STi/St  =  „  f  v  ( 18  inches  \  ,  *  ir,7  J  ”  * 

0,050  inch  ^ ^ - J  1,6  x  10  psi 

Values  of  the  correction  factors  for  these  material  configurations 
are  listed  in  Table  XI, 


TABLE  XI 

CORRECTION  FACTORS  FOR  EXAMPLE  PROBLEM 
ALUMINUM- ALUMINUM  AND  TITANIUM- STEEL  CONSTRUCTION 


2a/W 

2a 

(inches ) 

Aluminum- Aluminum 

Ai/Ai* 

Titanium-Steel 

A  Ti/St 

0.1 

1.8 

1.01 

1.01 

0.2 

3.6 

1.015 

1.01 

0.3 

5.4 

1.03 

1.015 

0.4 

7.2 

1.06 

1.02 

0.5 

9.0 

1.10 

1.025 

0.6 

10.8 

1.15 

1.03 

0.7 

12.6 

1.20 

1.025 

0.8 

14.4 

1.28 

1.02 

*  From  Figure  36, 


STEP  2,  For  crack  aspect  ratios  0.8>-^j->0.6,  as  the  crack  approaches  the 

edge  ZEE  stiffeners,  for  small  rivet  pitch  (assume  crack  grows  through 
rivet  and  center  stringer  remains  intact),  the  correction  factors  from 
Figure  46  apply.  Tabular  results  are  given  in  Table  XII, 


TABLE  XII 

CORRECTION  FACTORS  FOR  EXAMPLE  PROBLEM 
ALUMINUM- ALUMINUM  AND  TITANIUM- STEEL  CONSTRUCTION 


2a/W 

2a 

Aluminum- Aluminum 

Titanium-Steel 

(inches ) 

SAi/Ai** *** 

yt  .  ,  .irk* 

A  A/ /A/ 

c  ** 

5Ti/St 

%  it  irk* 

A  Ti/St 

0.6 

10.8 

5.2 

0.995 

1.55 

0.93 

0.7 

12.6 

5.5 

0.97 

1.6 

0.875 

0.8 

14.4 

5.8 

0.94 

1.7 

0.80 

0.85* 

15.3 

6.0 

0.92 

1.74 

0.76 

0.90* 

16.2 

6.1 

0.91 

1.8 

0.69 

0.95* 

17.1 

6.3 

0.80 

1.83 

~ .  6 

*  Analysis  is  invalid  for  i-2a  > stringer  width,  however,  proper  trend  is 
noted. 

**  From  Equation  IX-5. 

***  From  Figure  46. 
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STEP  3.  The  problem  can  now  be  analyzed  for  the  geometry  of  Figure  30(b) 

with  the  crack  emanating  from  the  rivet  hole  of  the  central  ZEE  strin 
ger.  The  correction  factors  from  Figure  33  for  the  two  cases  are 
shown  in  Table  XIII. 


TABLE  XIII 

CORRECTION  FACTORS  FOR  EXAMPLE  PROBLEM 
ALUMINUM-ALUMINUM  AND  TITANIUM-STEEL  CONSTRUCTION 


2a 

W 

2a 

(inches ) 

Aluminum- 

Aluminum 

Titanium-Steel 

SAi/A/* 

STi/St 

\  ii  irk'k 

A  Ti/St 

0.1 

1.8 

0.57 

0.75 

0.16 

0.71** 

0.2 

3.6 

1.15 

0.725 

0.33 

0.675** 

0.3 

5.4 

1.72 

0.715 

0.49 

0.675 

0.4 

7.2 

2.29 

0.73 

0.655 

0.67 

0.5 

9.0 

2.87 

0.73 

0.82 

0.67 

0.6 

10.8 

3.44 

0.73 

0.98 

0.67 

0.7 

12.6 

4.01 

0.74 

1.15 

0.675 

0.8 

14.4 

4.59 

0.74 

1.31 

0.675 

0.85 

15.3 

4.87 

0.745 

1.39 

0.68 

0.9 

16.2 

5.16 

0.75 

1.47 

0.68 

0.95 

17.1 

5.45 

0.75 

1.56 

0.68 

*  From  Equation  IX-3. 

**  Interpolated  from  Figure  33. 
***  From  Figure  33. 


STEP  4.  The  composite  correction  factors  for  the  aluminum/aluminum  and 

titanium/steel  construction  can  now  be  computed  by  superposition.  The 
results  for  the  crack  history  are  shown  in  Table  XIV. 

TABLE  XIV 

COMPOSITE  CORRECTION  FACTORS  FOR  EXAMPLE  PROBLEM 
ALUMINUM-ALUMINUM  AND  TITANIUM- STEEL  CONSTRUCTION  FROM  TABLES  XI-XIII 


2a 

2a 

( inches ) 

Aluminum- Aluminum 

Titanium-Steel 

W 

x*  =  X'xA"xA"' 

X*  =  A'xA"xA'" 

0.1 

1.8 

0.77 

0.72 

0.2 

3.6 

0.74 

0.68 

0.3 

5.4 

0.74 

0.685 

0.4 

7.2 

0.77 

0.68 

0.5 

9.0 

0.80 

0.69 

0.6 

10.8 

0.835 

0.629 

0.7 

12.6 

0.861 

0.596 

0.8 

14.4 

0.890 

0.548 

0.85 

15.3 

0.685 

0.517 

0.9 

16.2 

0.682 

0.469 

0.95 

17.1 

0.60 

0.408 
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STEP  5.  Using  the  composite  correction  factors  of  Table  XIV  and  the  range 
in  average  Kc  values  for  the  2219-T87  (115-125  ksiv/inch)  and  T1-6A1-4V 
(155-165  ksi^inch )  and  Equation  VIII-1,  compute  the  associated  critical 
stress  ranges  for  this  panel  geometry.  Comparison  of  the  two  proposed 
panel  fracture  stress  envelopes  is  shown  in  Figure  49. 

In  this  comparison  of  fracture  strength  to  yield  strength  ratio,  (this 
is  only  one  of  the  many  ways  of  presenting  fracture  data  comparisons)  the  all 
aluminum  construction  shows  as  the  better  fracture  safe  configuration  for 
•^<0.7.  However,  in  Figure  50,  a  comparison  of  fracture  stress  alone  indi¬ 
cates  the  superiority  of  the  2219-T87  construction  for  inspectable  size  cracks 
For  example,  the  dashed  line  could  be  the  anticipated  limit  design  stress  for 
the  2219-T87  and  the  solid  line  the  limit  design  stress  for  the  Ti-6A1-4V  con¬ 
struction.  Clearly  the  aluminum  construction  will  tolerate  a  7.2-inch  long 
crack  at  the  limit  design  stress  whereas  the  titanium/steel  construction  will 
permit  a  4  V2-inch  long  crack  for  a  90  percent  increase  in  limit  design  stress 
Both  crack  lengths  are  within  inspectable  limits. 

IX. 12  SUMMARY 


Thus  far  we  have  seen  how  fracture  mechanics  can  be  used  to  make  pre¬ 
liminary  estimates  of  fracture  stress.  Composite  correction  factors  which 
account  for  various  boundary  conditions  for  riveted  and  integral  stringer/ 
skin  construction  are  used  to  correct  for  crack  geometry.  All  of  these  ana¬ 
lyses  have  been  performed  with  the  assumption  of  limited  crack  tip  plasticity 
to  provide  conservative  results.  Therefore,  actual  material  fracture  data  for 
the  problem  or  design  of  interest  would  be  required  to  present  the  total  frac¬ 
ture  envelope  (i.e.,  same  thickness,  alloy  and  environment).  In  all  cases  the 
analysis  deals  with  through- the- thickness  crack  situations.  However,  it  is 
known  that  in  many  cases  (i.e.,  pressure  vessels,  thick  sections,  etc.)  the 
part  through  crack  situation  is  prevalent  in  service.  This  problem  will  be 
treated  in  the  next  section. 

IX. 13  PART -THROUGH  CRACKS* 

Normally  part-through  cracks  can  occur  in  pressure  vessels  and  initiate 
from  small  flaws  or  discontinuities  which  may  or  may  not  manifest  themselves 
during  proof  testing.  Flaw  growth  can  occur  during  subsequent  pressurization 
cycles  and  lead  to  catastrophic  failure.  One  of  the  prime  difficulties  in 
analyzing  for  fracture  is  the  presence  of  environmental  factors  (corrosive 
media)  which  may  enter  into  the  problem.  Fortunately,  fracture  mechanics  can 
be  used  with  good  success  to  analyze  these  structures.  A  great  amount  of 
research  has  been  directed  at  this  problem  due  to  its  historical  permanence 
(see  Section  II). 

It  is  necessary  to  show  that  the  largest  inspectable  initial  flaw  will 
not  become  critical  during  its  life  span.  The  crack  size  requirements  are 
dictated  by  nondestructive  testing  (NDT)  and  the  parameters  by  fracture  mecha¬ 
nics  (stress  levels,  fracture  toughness,  material  thickness  and  location  and 
orientation  of  flaws. 


*  Although  discussion  of  this  crack  problem  is  usually  oriented  towards  pres¬ 
sure  vessel  applications,  the  part-through  crack  can  also  occur,  and  the 
same  analysis  applies  to  any  structure  with  these  crack  conditions. 
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9.0 


2219-T87  SKIN/STRINGER 
^ys  =  59  ksi 
Ti-6A1-4V  SKIN/STEEL 
STRINGER  OyS  =  125  ksi 


CRACK  ASPECT  RATIO,  ^ 

-JL. 


+ 

4.5  0  4.5 

CRACK  LENGTH,  a  (INCHES) 


9.0 


FIGURE  49  ESTIMATES  OF  FRACTURE  ENVELOPES  FOR  THE  PANEL  GEOMETRY  OF 
FIGURE  48  FOR  TWO  SKIN  MATERIALS  (NO  CRACK  BUCKLING) 
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Ti-6A1-4V  SKIN/ 
STEEL  STRINGER 


,  r~  i  1  u  I I I 

1.0  0.5  0.1  0.1  ,  0.5  1.0 

CRACK  ASPECT  RATIO, 

W 


j- | | | 1 

9.0  4.5  0  4.5  9.0 

CRACK  LENGTH,  a  (INCHES) 

FIGURE  50  COMPARISON  OF  FRACTURE  STRESS  L  CRACK  LENGTHS  FOR  PANEL 
GEOMETRY  OF  FIGURE  48  FOR  TWO  SKIN  MATERIALS  (NO  CRACK 
BUCKLING) 
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IX . 1 3 ( a )  Determination  of  Critical  Discontinuity  Dimensions  and  Kj  for 
Part-Through  Cracks 


Part-through  cracks  usually  manifest  themselves  as  surface  or  embedded 
flaws.  For  thick  wall  structures,  pressure  vessels  for  example,  with  surface 
discontinuities  or  small  embedded  flaws,  the  following  equation  applies  for 
critical  conditions. 


K 


Ic 


(IX -6) 


Irwin  (41),  first  described  the  opening  mode  stress  intensity  solution  for  a 
part-through  elliptical  crack  as 


K 


I 


A  (TjJr 


(1X-7  ) 


where(§)in  both  cases  is  used  to  describe  the  flaw  size.  The  relationship 
between  flaw  shape  (Q)  and  discontinuity  aspect  ratio  (depth  length,--  )  is 
given  in  Figure  51*  In  equation  IX-6,  the  constant  C  =  =  1.21  c  for 

surface  discontinuities  and  unity  for  small  embedded  flaws.  Figure  51.  shows 
the  flaw  shape  relationship  for  various  operating  to  yield  stress  ratios. 

Notice  that  the  critical  dimension  for  part-through  cracks  is  always  through- 
the- thickness ,  therefore,  a,  is  the  crack  depth  dimension  for  this  class  of 
crack  situation. 

Equation  IX-6  has  been  found  quite  accurate  for  *§<  50%  inhere  B  is  material 
thickness);  however,  at  greater  crack  depths  the  stress  intensity  solution  of 
Equation  IX-7  must  be  modified  due  to  the  stress  intensity  magnification  as 
the  flaw  reaches  a  free  surface  (i.e.  ^>50%).  Therefore,  an  elastic  magnifica¬ 
tion  factor  (M^)  is  applied  to  K  which  varies  with  a/B  as  shown  in  Figure  52 
and  is  described  functionally  by  Equation  IX-8. 


K 


I 


Acrj7r 


M  [\ 
KV  Q 


(IX-8) 


This  implies,  for  example,  that  in  “thin”  walled  pressure  vessels  the  stress 
intensity  becomes  large  (see  Equation  IX-8  and  Figure  52) for  small  flaw  sizes. 

Or  Ki — ^^critical  faster  for  a  given  flaw  size  when  ^  >0.5. 

The  elastic  magnification  factors  for  semi-elliptical  surface  flaws  in  a  section 
in  bending  at  the  crack  surface  (  /3  =  90°)  and  root  ( (3  =  0)  for  <  0.5  are 
given  in  Figure  53.  And  the  magnification  factors  for  embedded  flaws  in  finite 
thickness  plate  as  a  function  of  eccentricity  and  crack  depth  aspect  ratio  is 
shown  in  Figure  54.  The  combined  elastic-plastic  magnification  factors  for 
deep  surface  flaws  are  given  in  Figure  55,  (i.e.,  "thin”  wall  cylinders). 

IX. 13(b)  Arm  In  Bending  as  a  Part-Through  Crack  Problem 

In  analyzing  the  arm  of  Figure  22,  it  is  possible  that  a  par t- through 
(elliptical)  crack  occurs  at  point  C  (shown  as  insert  in  Figure  22).  This 
part  has  been  previously  treated  as  a  plane  bending  problem.  We  will  now  use 
part-through  crack  principles  and  fracture  mechanics  and  compare  these  results 
with  the  through- the- thickness  crack  situation. 
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DISCONTINUITY  ASPECT  RATIO 


0.5-r 


i  ° 

HJ  |  CN 


0.4- 


0.3- 


0.2- 


0.1- 


°4- 

0.6 


FLAW  SHAPE  PARAMETER,  Q 


b- 2c-H  ra 


SURFACE  FLAW 


EMBEDDED  FLAW 


FIGURE  51  FLAW  SHAPE  PARAMETER  CURVES  (REF.  45) 
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ELASTIC  BENDING  MAGNIFICATION  FACTOR, 


-'Original  Sources  -  Smith,  F.  W.  ,  ’’Stresses  Near  Semi-Circular 
Edge  Crack,”  Ph.D.  Dissertation,  Univ.  of  Washington,  1966. 

Smith,  F.  W.  ,  ’’Stress  Intensity  Factors  for  a  Semi- 
Elliptical  Surface  Flaw  in  Plate,”  Structural  Development 
Research  Memo,  The  Boeing  Company,  1966. 

Gross,  B.  J.  &  Srawley ,  J.  E. ,  NASA  TND-2603,  1965. 


FIGURE  53  STRESS  MAGNIFICATION  FACTORS  FOR  SEMI -ELLIPTICAL  SURFACE  FLAW  IN 
SECTION  IN  BENDING  AT  y8  =  0°  &  90°  (REF.  25)* 
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ELASTIC  MAGNIFICATION  FACTOR 


T 

B 

L 


x 


0.35 


0.25 


1.0 


0 


I  ~  I  '  -  I 

0.1  0.2  0.3 


T  \ 

0  v4  0.5 


FLAW  ECCENTRICITY,  (1  -  -^) 

B 


FIGURE  54  ELASTIC  MAGNIFICATION  FACTORS  AT  POINTS  Pt  &  Pt  FOR  AN  EMBEDDED, 
ECCENTRIC  LOADED  FLAW  IN  TENSION  (REF.  49 )L  2' 
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COMBINED  ELASTIC  6.  PLASTIC  MAGNIFICATION  FACTORS 


FIGURE  55  COMBINED  ELASTIC  &  PLASTIC  MAGNIFICATION  FACTORS  FOR  A  DEEP 
SURFACE  DISCONTINUITY  -  IN  TENSION  <  0.3)  (REF-  25)* 
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Assume  that  a  surface  crack  starts  at  point  C  (Figure  22),^  and  has  a 
semi-elliptical  shape  through  the  thickness.  In  this  case,  we  will  be  analyz¬ 
ing  a  part-throu2h  crack  situation  with  the  shape  of  Figure  52  with  tensile 
loading.  Equation  IX-8  applies  and  the  ratio  of  crack  depth  to  thickness,  ^ 

is  ^^ineh0^  =  ^•^5.  Assume  that  the  discontinuity  aspect  ratio  -j—  ^  3^4  In* 

^0.33,  the  elastic  bending  magnification  factor,  Mg,  is  from  Figure  53  equal 
to  0.9.  This  value  is  not  unexpected  as  it  is  approaching  a  through-crack 
geometry  and  the  correction  factor  will  tend  toward  unity. 

STEP  1.  Determine  the  flaw  shape  parameter,  £.  The  yield  strength  for  the 
7075-T6  alloy  is  65  ksi  and  the  2219  alloy,  50  ksi.  The  nominal 
(tensile)  stress  for  the  arm  at  the  crack  position  is  by  simple  stress 
analysis  12  ksi  and  18  ksi  for  the  2,000  and  3,000  lb.  force,  respec¬ 
tively.  The  values  of  the  flaw  shape  parameters  are  given  in  Table  XV 
for  the  operating  to  yield  stress  ratios  and  a  discontinuity  aspect 
ratio  (a/2c)  of  0.33. 


TABLE  XV 

FLAW  SHAPE  PARAMETERS  FOR  EXAMPLE  PROBLEM 


MATERIAL 

-2-  (2000  lb  force) 
cr 
ys 

Q* 

(3000  lb  force) 
ys 

Q* 

7075-T6 

0.18 

1.75 

0.28 

1.73 

2219-T851 

0.24 

1.73 

0.36 

1.70 

*  From  Figure  51. 


STEP  2.  From  Equation  IX-8,  A  is  1.1  for  a  surface  crack.  The  value  of  Kj 
can  now  be  determined  for  the  1/4"  deep  surface  flaw.  For  a  2,000  lb. 
force  and  7075-T6  material 


K  =  1.1  (12  ksi)  s/F  (0.9) 


Vk  inch 
1.75 


8.0  ksi^/Tnch 


and  for  a  3,000  lb.  force  and  7075-T6  material 


Kj  =  1.1  (18  ksi)  JW  (0.9)  yj  =  12.0  ksi^inch 

The  stress  intensities  for  the  2219-T851  material  are  8.0  ksi^/inch 
(2,000  lb.  force)  and  12.1  ksi^/inch  (3,000  lb.  force).  The  stress 
intensities  for  a  2,000  and  3,000  lb.  force  and  through  the  thickness 
crack  obtained  previously  are  11.1  ksi^/inch  (2,000  lb.  force)  and 
16.8  ksi^/inch  (3,000  lb.  force). 


Thus,  results  obtained  for  the  semi-elliptical  part- through  crack  com¬ 
pared  to  the  through  crack  indicate  slightly  lower  (<  3  ksiVinch)  stress 
intensities.  However,  the  geometric  corrections  approach  unity  in  both  cases. 


125 


IX. 14  APPLICATION  TO  DESIGN  -  PART-THROUGH  CRACKS 


An  excellent  summary  of  the  application  of  fracture  mechanics  to  the 
design  of  pressure  vessels  is  given  in  Reference  42.  Those  concerned  with  the 
analysis  of  part-through  surface  or  embedded  flaws  should  refer  to  that  docu¬ 
ment  for  the  limitations  of  the  analysis  as  well  as  complete  illustrative 
examples . 


For  a  given  environment  and  loading  profile,  the  time  or  cycles  to 
failure  of  a  pressure  vessel  depends  primarily  on  the  initial  stress  intensity 
(Kii)  for  an  assumed  discontinuity  compared  to  the  material  plane  strain  frac¬ 
ture  toughness.  That  is  the  number  of  cycles  to  failure  or  time  to  failure 
is  a  function  of  ^li  .  Through  the  use  of  emperical  data,  curves  of  ^li  can 

be  established  as  ^cshown  in  Figure  56  for  static  loaded  and  cyclic  lol&ed 
structure.  For  the  static  loaded  case,  it  can  be  seen  that  below  some  thresh¬ 
old  value  or  £th  that  the  crack  or  flaw  will  not  grow  as  load  is  sustained 

klc  Klc 

indefinitely;  whereas,  if  the  initial  flaw  (and  corresponding  K^)  ab°ve 
the  threshold  value,  the  crack  will  grow  with  time  as  constant  load  (stress) 
sustained  and  failure  will  ensure.  For  the  cyclic  load  case,  no  threshold 
value  is  evident  and  some  crack  growth  and  eventual  failure  can  be  expected 
unless  a  non-propagating  crack  situation  develops.  Much  variation  in  thresh¬ 
old  value  has  been  experimentally  detected  and  is  caused  by  alloy,  heat  treat, 
temperature,  etc.  The  techniques  used  in  employing  threshold  stress  intensity 
data  in  design  and  analysis  are  reported  in  Reference  42  in  detail. 
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IX. 15  EFFECT  OF  BULGING  DUE  TO  INTERNAL  PRESSURE 


In  thin  gage,  pressurized  cylinders  it  is  known  that  the  stress  inten¬ 
sity  of  a  through  crack  is  influenced  not  only  by  the  membrane  stresses,  but 
also  by  bulging  of  the  material  on  either  side  of  the  crack. (50)  The  crack 
geometry  of  this  problem  is  shown  in  Figure  57.  This  problem  has  been  treated 
in  Reference  50  as  a  stress  intensity  superposition.  When  the  crack  is 
oriented  parallel  to  the  cylinder  axis,  the  superimposed  stress  intensity  is 
the  sum  of  the  factors  for  hoop  membrane  (K^)  and  bulging  (Kg),  or,  (see  also 
Section  V.  ) 

K  -  *h  +  \  ( IX-9 ) 


The  complete  solution,  from  Reference  50,  is 


K  = 


(IX-10) 


where  C  and  ft  are  empirically  determined  constants  obtained  from  pressurized 
burst  tests  on  precracked  cylinders,  and  R  is  the  cylinder  radius.  Equation 
IX-10  does  not  include  correction  for  crack  tip  plasticity.  The  plasticity 
corrected  form  of  Equation  IX-10  is  (see  also  Section  VI). 


1  +  c  |  (f)"l  (IX-11) 

m 

In  equation  IX-11,  cr  g  is  the  material  yield  strength  in  a  2:1  biaxial 
stress  field.  The  constants  C  and  jl  for  the  materials  tested  in  Reference  50 
were  found  to  be  equal  and  can  be  computed  from  pressure  cylinder  fracture  data. 
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RATIO  OF 


CYCLES  TO  FAILURE 


FIGURE  56  SCHEMATIC  OF  STATIC  AND  CYCLIC  FLAW  GROWTH  FOR  A 
GIVEN  ENVIRONMENT  (REF.  42) 


127 


<7h  MEMBRANE  HOOP  STRESS 
( =  MEMBRANE  LONGITUDINAL  STRESS 
P  INTERNAL  PRESSURE 


FIGURE  57  CRACK  GEOMETRY  &  STRESS  SYSTEM  FOR  PRESSURIZED  CRACK  WITH 
PRESSURE  BULGING  (REF.  50) 
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Thus,  for  critical  conditions,  a - *"ac>  K - ►  Kc  and  O’  — -,cr^  the  bulge 

_  c 

coefficient  C  is  f“  ^ 

c 


R 

c  -  M  =  r 


K  2 
c 


—  1 


"V  2^72 


yB 


(IX-12) 


Rearranging  terms  and  solving  for  the  critical  hoop  stress  (  CT  )for  critical 

he 

conditions 


K 


c 


(IX -13) 


or 


(IX -14) 


based  on  knowledge  of  the  original  crack  length,  aQ,  only.  In  this  case, 
Kc  is  the  nominal  plane  stress  fracture  toughness  of  the  material  obtained 
from  flat  sheet  fracture  tests. 


In  equations  IX-13  and  IX-14  the  geometric  correction  is  A  ==(1  +  C  . 

R 

Values  of  the  bulge  coefficient  C,  based  on  original  crack  length  a  ,  as 
given  in  Reference  50  are  repeated  in  Table  XVI.  It  should  be  pointed 
out  that  C  was  empirically  determined  in  Reference  50  and  its  dependence  on 
temperature  or  material  is  not  clearly  understood. 
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TABLE  XVI 

VALUES  OF  BULGE  CORRECTION  AT  TEMPERATURE  FOR  VARIOUS  MATERIALS 

(DATA  FROM  REF.  50) 


Material 

Test 

Temperature 

Bulge  Correction 
(C) 

2014-T6 

Room 

6.4 

2014-T6 

-320°  F 

3.5 

2014-T6 

-423°F 

2.9 

T1-5A1-2. 5  Sn  (ELI) 

-320°  F 

4.3 

T1-5A1-2. 5  Sn  (ELI) 

-423°F 

1.6 

2024-T3 

Room 

9.5 

7075-T6 

Room 

9.5 

IX. 15(a)  Example  Problem 

Compare  the  critical  hoop  tension  fracture  stress  (0"hc)  for  a  20-inch 
diameter  pressure  bottle  which  can  withstand  a  V2-inch  long  crack  without 
leaking  catastrophically  or  bursting  for  two  materials  which  are  being  consi¬ 
dered,  2014-T6  ( O"  =  68  ksi)  and  7075-T6  (°~yS  =  ksi)  aluminum.  The 
material  thickness  is  0.060  inch  for  the  2014-T6  and  7075-T6  cylinders ,  with 
nominal  YLq  of  58.4  ksi^inch  for  2014-T6  aluminum  and  Kc  of  85  ksi-Jinch  for  the 
7075. 


Using  the  bulge  corrections  from  Table  XVI  and  Equation  (IX- L4),  deter¬ 
mine  the  critical  hoop  tension  stresses.  (The  biaxial  yield  strength  from 
Von  Mises  yield  criteria  can  be  assumed  to  be  ~1.15  cr yS  --  see  Ref.  50.) 


2014-T6 


cr 


58.4  ksiy/inch 


he 


/  0.5  inch  .  (58.4  ksi*/inch)  ,,  t  ,  , 
f 7T  5 -  +  - * - t  (  1  +  6.4 

2  [l .  1 5  (68  ksi)]Z 


0.5  in.' 


2SL Jju. 


=  48.8  ksi 


7075-T6 


cr 


he 


85  ksiy/inch 


0>5  inch  (85  ksi^inch)  [  l  +  9  5 
2  2 [l .  15  (77  ksi)]2 


0.5  in. 


20  in. 


=  61.5  ksi 
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The  7075-T6  material  would  be  a  better  choice  on  a  fracture  stress  basis,  or 
a.  fracture  stress  to  biaxial  ratio,  ^hc  ~  0.7  compared  to  0.6  for  the  2014- 
T6  material.  ayB 

Recently,  the  authors  of  Reference  51  have  compared  the  various  methods 
available  to  fracture  analyze  pressurized,  through- the- thickness  cracks.  This 
analysis  included  the  data  obtained  in  Reference  50.  Better  correlation  was 
obtained  if  a  criteria  was  set  on  the  crack  tip  plasticity.  The  general  form 
of  the  equation  presented  in  Reference  51  for  long  cracks  is, 

K 

a-  =  - - -  (IX- 15) 

c  s/iriy  s/~£~ 

where  7  is  a  plasticity  correction  and  C  is  the  magnification  factor  applied 
to  lioop  stress  which  provides  nominal  stress  values  so  comparisons  can  be 
made  with  flat  panel  fracture  behavior.  In  this  case,  the  Kc  value  is  ob¬ 
tained  from  a  pressure  vessel  burst  test.  Hahn,  et.  al.,  (51)  have  proposed 
that  for  "thin"  wall,  low  to  intermediate  toughness  vessels,  with  relatively 
long  cracks  (i.e.,  rocket  tankage  and  pressure  fuselages)  f  can  be  expressed 
as 

£  -  J  1  +  l.«  ^  (SO  tanh^) 

The  plasticity  correction  y  is  given  by 


7 


17 ‘VS 

2  cr 


In 


sec 


^hcS 

2  cr 


(IX-16) 


In  Equation  IX-16,  cr  is  the  average  flow  stress  and  is  not  precisely  defined. 
The  upper  and  lower  bounds  are  given  as  <XyS  <  cr  <  c ru. 

Other  fracture  equations  for  intermediate  thickness  vessels  with  long 
and  short  crack  situations  (e.g.,  steel  pipe  lines  and  pressure  vessels  operat¬ 
ing  at  low  temperatures)  are  presented  in  Reference  51. 


IX. 16  APPLICATION  TO  REINFORCED  PRESSURE  CYLINDERS 


In  pressurized  fuselage  structure,  the  analysis  of  References  50  and  51 
which  were  just  presented,  can  be  seen  to  have  direct  application.  The  pri¬ 
mary  loads  on  these  structures  are  the  membrane  hoop  tension  due  to  cabin 
pressurization,  and  bending  loads  from  maneuver,  gust,  landing  and  take-off. 
From  Equation  IX-10,  (Reference  50)  the  crack  bulging  correction,  [hc  ^ 

is  seen  as  a  crack  geometry  correction,  X ,  which  can  be  treated  in  the  same 
manner  as  other  geometric  corrections.  For  example,  it  would  be  possible  to 
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analyze  riveted  skin-stringer  fuselage  section  designs  where  reduction  due  to 
pressure  bulging  would  occur  due  to  sudden  puncture.  Methods  employed  earlier 
in  this  section  to  obtain  composite  correction  factors  for  skin-stringer  de¬ 
signs  could  then  be  used  in  combination  with  the  bulge  factor.  Thus,  pressuri¬ 
zation  effects  can  be  accounted  for  in  calculating  fracture  stress  envelopes. 
For  large  fuselage  diameters  (i.e.,  R»a),  it  can  be  seen  that  the  bulge 
correction  factor  (  1  +  C  t)  approaches  1.0  or  infinite  plate  behavior.  How¬ 
ever,  in  failsafe  analysis,  long  crack  lengths  are  common  and  the  reduction  in 
fracture  stress  due  to  crack  bulging  can  be  extensive  (depending  on  material) 
and  should  enter  into  the  overall  fracture  mechanics  analysis. 

IX. 17  CRACK  ARREST 


A  frequent  objective  in  structural  design  is  to  provide  some  technique 
for  limiting  or  arresting  c rack  growth  so  that  catastrophic  failure  can  be 
prevented.  Requirements  for  fail-safe  structure  necessitate  the  development 
and  application  of  such  techniques. 

There  are  two  ways  in  which  it  can  be  attempted  to  arrest  a  running 
crack;  one  material  oriented,  the  other  by  structural  design  considerations. 
Both  of  these  approaches  have  one  tiling  in  common,  that  is  reducing  the  stress 
(and,  therefore,  stress  intensity)  near  the  crack  tip.  Material  changes  which 
lead  to  higher  toughness  materials  can  assist  in  making  crack  arrest  easier. 
However,  when  a  material’s  higher  toughness  is  combined  with  designs  which 
incorporate  structural  members  in  a  redundant  load  carrying  structure,  the 
optimum  possibility  of  crack  arrest  is  achieved  (see  Section  III). 

There  is  a  popular  and  effective  concept  frequently  employed  in  commer¬ 
cial  aircraft  design,  which  utilizes  crack  stoppers  or  tear  straps  at  strate¬ 
gic  points  in  a  load  bearing  structure.  These  straps  of  thin  gage  aluminum 
or  titanium  are  normally  adhesively  bonded  to  the  main  structure.  In  princi¬ 
ple,  they  act  as  a  point  of  stress  redistribution  for  an  advancing  crack.  The 
increase  in  effective  thickness  due  to  the  layering  effect  is  accomplished 
with  no  loss  in  fracture  toughness.  In  other  words,  the  fracture  toughness  of 
the  layered  structure  is  on  the  order  of  the  one  layer  itself  and  does  not 
change  with  layering  thickness  (see,  i.e.,  Reference  52).  Geometrical  boun¬ 
daries  (splices,  etc.)  also  aid  in  arresting  a  running  crack  and  are  often 
employed  as  crack  arrest  barriers. 

Reference  53  contains  a  description  of  the  c rack  arrest  problem  and 
summarizes  the  research  which  has  been  performed  in  this  area.  The  crack 
arrest  problem  in  riveted  sheet-stringers  is  analyzed  in  Reference  54.  Addi¬ 
tional  sources  can  be  located  in  the  Bibliography  of  this  report. 


IX. 18  FRACTURE  OF  COMPOSITES 


Research  into  the  fracture  of  composite  materials  has  been  limited  until 
recently  (see,  for  example,  References  55  and  56).  However,  with  the  greater 
employment  of  these  materials  in  aerospace  structural  applications,  additional 
data  will  become  available  on  both  crack  propagation  and  fracture  strength 
which  can  be  used  in  design.  In  general,  it  is  known  that  the  fracture  process 
in  fiber  composites  is  discontinuous.  The  ductile  matrix  material  tends  to 
aid  in  crack  branching  along  the  fiber  at  the  fiber/matrix  interface  and  leads 
to  increasing  the  overall  fracture  strength  in  uniaxial  tension  loading  parallel 
to  the  fibers. 
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X  FATIGUE  AND  CRACK  PROPAGATION 


X#1  INFLUENCE  OF  FATIGUE  ON  FRACTURE 


As  mentioned  in  earlier  sections  of  this  report,  fracture  mechanics 
assumes  the  preexistence  of  flaws,  inhomogeni ties  and  discontinuities  in  a 
material.  It  is  from  these  flaws  that  cracks  often  nucleate  and  grow  under 
cyclic  load.  The  prediction  of  crack  nucleation  is  beyond  the  scope  of  this 
report;  however,  it  must  be  recognized  that  any  local  geometric,  physical  or 
metallurgical  feature  which  intensifies  the  stress  is  a  prospective  site  for 
fatigue. 

During  stress  or  strain  cycling  of  a  material,  there  is  a  "shake-down" 
period  of  a  few  number  of  cycles  (dependent  on  magnitude)  during  which  the 
material  will  cyclic  harden  or  soften.  (57)  For  a  period  of  time,  the  crack 
will  nucleate  and  then  grow  by  stages.  (58)  In  fracture  mechanics  we  are 
interested  in  the  propagation  of  engineering  cracks  of  detectable  size.  Ihe 
prediction  of  crack  growth  is  discussed  later  in  this  section. 

X.2  MICRO-MECHANICS 

It  is  possible,  f rom  microf ractography ,  to  determine  not  only  the  site 
of  fatigue  but  the  subsequent  stress  history  as  long  as  the  plastic  zone  is 
small  at  the  crack  tip.  This  is  usually  the  case  in  high  cycle  fatigue  and/or 
brittle  material  which  produce  flat  (plane  strain)  type  fracture  surfaces.  An 
example  of  an  electron  microscope  fractograph  taken  from  an  area  on  a  7075-T6 
aluminum  fracture  surface  where  a  crack  has  grown  is  shown  in  Figure  58.  The 
direction  of  crack  propagation  is  indicated  by  the  arrow.  The  rippled  appear¬ 
ance  (striations)  is  the  result  of  successive  fatigue  cycles,  and  the  distance 
between  successive  striations  can  give  a  measure  of  the  crack  growth  rate  (per 
cycle).  Through  the  use  of  fracture  mechanics,  it  is  then  possible  to  deter¬ 
mine  operating  stress  and  cycles  to  fracture.  The  Bibliography  of  this  report 
lists  documents  which  should  be  consulted  for  additional  information. 

X . 3  FRACTURE  MECHANICS  AND  CYCLIC  CRACK  GROWTH 

In  order  to  establish  a  method  of  inspection  intervals  for  aircraft 
structure,  it  would  be  desirable  if  a  means  were  available  to  predict  succes¬ 
sive  crack  lengths  during  service  history.  During  preliminary  design,  it 
would  be  helpful  if  a  method  of  material  comparison  (in  addition  to  conven¬ 
tional  S-N  data)  were  used  to  rank  material  on  the  basis  of  cyclic  crack 
growth  resistance.  For  example,  it  has  been  shown  (Refs.  59-62)  that  the 
stress  intensity  factor  is  a  good  correlating  factor  for  crack  propagation 
under  fluctuating  as  well  as  static  loading. 

It  was  previously  indicated  that  microscopic  crack  propagation  takes 
place  on  a  cyclic  scale  in  tune  with  the  fluctuating  load  or  stress.  On  this 
basis,  it  is  reasonable  to  assume  that  incremental  crack  extension  (A  a)  in  a 
small  number  of  cycles  (AN)  is  related  to  the  slope  (derivative)  of  the  con¬ 
tinuous  crack  length  (a)  vs.  cycle  (N)  curve,  =  iifL  . 

AN  dN 
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FIGURE  58  ELECTRON  FRACTOGRAPHIC  OF  FATIGUE 
AREA  ON  7075-T6  ALUMINUM  PLATE 
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The  growth  rate  J  is  a  function  of  (1)  Material ,  (2)  Environment,  (3) 
Load-Time  History,  and  (4)  Geometry,  including  crack.  Stress  intensity  can 
describe  items  3  and  4  with  a  single  parameter,  K.  The  two  remaining  items 
will  be  discussed  later.  Paris(62)  showed  that  the  stress  intensity  factor 
(K)  describes  the  effect  of  external  load  and  crack  geometry  for  cyclic  crack 
growth.  For  this  reason,  he  showed  that  =  £  (K).  Since  K  is  a  linear  fac¬ 
tor  in  the  elastic  stress  equations,  it  depends  on  the  load  or  stress  magni¬ 
tude  and  crack  length  in  a  linear  fashion.  If  K(t)  and  cr(t)  are  time  history 
of  load  (stress)  and  stress  intensity,  and  if  the  crack  length  varies  slowly 
with  time  (t),  the  local  crack  tip  stress  intensity  history  is  K(t)  =  P(t) 

.  f(a)  (62).  pn  terms  of  stress  for  a  center  crack  (infinite)  plate  (see  Eqs. 
V- 25 ) ,  K(t)  =  cr(t)  •  s/tFT. 

Therefore,  the  time  dependent  character  of  fatigue  crack  growth  can  be 

expressed  as 


f-K(t)  (X-l) 

The  reason  that  fracture  mechanics  analysis  has  shown  greater  applicabi¬ 
lity  to  fatigue  crack  propagation  is  due  to  the  conditions  for  fatigue  being 
less  than  critical.  In  other  words,  the  fatigue  cracking  process  involves 
stress  levels  and  crack  lengths  (hence,  Kfs)  below  the  static,  critical  values. 
Work  hardening  during  cyclic  loading  raises  local  (crack  tip)  yield  strengths 
and  produces  small  yield  zones  (see,  e.g.,  Equation  VI-2).  Therefore,  elastic 
stress  field  solutions  are  more  directly  applicable  due  to  smaller  plastic 
zones . 


X . 4  SINUSOIDAL  LOADING 


In  fatigue  testing,  sinusoidal  loading  is  quite  common.  K(t)  local  to 
the  crack  is  then  a  sine  wave  about  some  mean  value  of  K  and  variations  in  K 
can  then  be  generally  described  by  a  sinusoidal  amplitude,  a  mean,  and  fre¬ 
quency. 

Many  equations  have  been  proposed  for  the  prediction  of  fatigue  crack 
growth  based  on  fits  to  empirical  data  (Bibliography).  These  equations 
appear  to  be  valid  as  long  as  they  are  applied  within  the  test  limits  and 
conditions  from  which  they  are  derived.  In  1963,  Paris  and  Erdogan^bu^ 
reviewed  existing  crack  propagation  equations  and  from  a  wide  range  of  test 
data  proposed  that  a  best  fit  might  be: 


da 

dN 


£  (AK)n  , 


(X-  2 ) 
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where  C  is  a  material  constant  and  AK  is  the  stress  intensity  factor  range 
(Kmax.  —  Kmin.)'  Equation  X-2  suggests  a  straight  line  data  fit,  with  n  =  4 
for  log  AK  vs.  log  da/dN.  With  this  relationship,  numerous  researchers 
attempted  correlations  based  on  Equation  X-2  (see,  e.g.,  Reference  63),  but  it 
became  obvious  that  a  single  equation  could  not  account  for  large  differences 
in  fatigue  stress  ratio.  To  account  for  this  parameter,  the  authors  of  Refer¬ 
ence  64  presented  equation  X-3  as  a  better  data  fit  which  would  account  for 
stress  ratio  effects  and  apply  to  the  point  of  fracture. 

da  C(A  K)n  ,  . 

dN  (l-R)K-AK 
c 

In  1967,  Hudson^^  empirically  compared  equation  X-3  with  other  equations  for 
crack  growth  and  concluded  that  it  represented  the  data  quite  well.  Some  data 
from  Reference  65  has  been  abstracted  and  is  shown  in  Figure  59  for  7075-T6 
and  2024-T3  aluminum  alloy  sheet  tested  at  two  cyclic  stress  ratios. 

However,  it  must  be  pointed  out  that  environmental,  metallurgical,  multi 
axial  stress  and  frequency  effects  appear  only  as  lumped  parameters  in  the 
material  constant  (C)  and  numerical  exponent  (n)  of  equations  X-2  and  X-3. 
(Typical  values  of  C  and  n  for  7075-T6  and  2024-T3  are  shown  in  Figure  59.) 

The  individual  effect  of  these  parameters  on  fatigue  crack  propagation  have 
not  as  yet  been  separated. 

X. 5  RANDOM  LOADING 


Crack  propagation  data  from  random  load  tests  have  indicated  that  corre¬ 
lation  on  a  stress  intensity  basis  is  not  only  feasible  but  adequately  des¬ 
cribes  the  crack  growth  for  such  loading  (see,  for  example,  Reference  62). 

The  data  of  Swanson' ^ '  on  7079-T6  aluminum  indicates  that  fracture  mechanics 
was  the  most  effective  approach  to  analyzing  crack  propagation  for  both  con¬ 
stant  amplitude  and  Rayleigh,  random  amplitude  loading. 


X. 6  PREDICTION  OF  CRACK  GROWTH  TO  FRACTURE  IN  DESIGN 


Estimates  can  be  made  of  crack  growth  or  extension  using  Equation  X-2  or 
X-3.  For  uniform,  sinusoidal  loading  in  a  plate  (small  crack  lengths)  integra 
tion  of  Equation  X-3  was  suggested  by  Forman,  et.  al.,  to  compute  cycles 

to  failure  for  a  given  fatigue  crack, 


N 

c 


2 

7TC(1-R)2o-2 


K 

1  '^n  r 

o 


(X-4 ) 


where  N  is  the  critical  value  of  cycles  at  crack  instability,  N0  is  the  ini¬ 
tial  value  of  cycles,  and  KQ  is  the  initial  crack  tip  stress  intensity,  i.e., 

j  (ao) • 


We  will  examine  a  problem  in  crack  growth  in  two  ways,  growth  to  frac¬ 
ture  and  cycles  to  crack  detection. 
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CRACK  GROWTH  RATE  (INCHES/CYCLE)  CRACK  GROWTH  RATE  (INCHES/CYCLE) 
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FIGURE  59  TYPICAL  CRACK  GROWTH  VS  STRESS  INTENSITY  FACTOR  RANGE 

DATA  FOR  7075-T6  &  2024-T3  ALUMINUM  (DATA  FROM  REF.  65) 
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X.6(a)  Example  Problem 


A  corner  crack  has  developed  at  the  top  surface  of  a  bolt  hole  in  the  part 
shown  in  Figure  60.  The  plate  thickness,  B,  is  1  inch  and  the  material  is  a 
100  ksi  yield  strength  alloy  steel.  The  plane  strain  fracture  toughness  (Kjc) 
is  assumed  to  be  130  ksiVinch.  The  operating  cyclic  stress  of  the  part  in 
question  is  sinusoidal  with  a  maximum  stress  of  70  ksi  and  a  minimum  of  1.4 
ksi  (R  of  0.02). 

Question:  Assume  that  a  0.3  inch  deep,  semi-circular  crack  grows  under  the 

bolt  heat,  how  many  cycles  will  it  take  for  the  crack  to  grow 
through  to  the  back  surface  (through  the  plate  thickness)?  Also, 
how  many  cycles  will  be  allowed  before  critical  conditions  prevail, 
assuming  that  the  crack  changes  from  a  semi-circular  to  a  straight 
front- through  the  thickness  crack  after  penetrating  the  thickness. 

STEP  1.  To  answer  the  first  part  of  the  question,  it  is  necessary  to  assume 
that  the  crack  depth  (ac)  is  reached  when  the  crack  grows  through  the 
plate  thickness.  Thus,  when  the  crack  length  ac  is  1  inch,  compute  a 
critical  or  Kq.  From  Equation  VII-7  when  Kj — ►Kq  with  tjqtt  =  0: 

=  0 . 705rrV7ra  , 

Q  c 

=  0.705  (70  ksi)  V  77(1  inch)  =  87.5  ksiVinch' 

Using  Equation  X-4,  solve  for  the  number  of  cycles  for  the  initial  crack 
(0.3  inch)  to  grow  through  the  thickness  (1  inch),  where  Kq  is  87.5  ksi 
Vinch.  Values  of  the  constant,  C,  have  not  been  reported  for  this  con¬ 
structional  steel,  but  we  may  assume  that  it  is  on  the  order  of  10 
(see,  e.g.,  Figure  59  and  Refs.  64  and  65  for  aluminums). 


N 

c 


_ 2 _ 

77 10'13(l-0.02)2(70,000psi)2 


_ 87,500  psiv/inch _ 

0.705  (70,000psi)s/7r  (0.3  inch) 


1 


-  in 


_ 87,500  psi\/inch _ 

0.705  (70,000psi)2s/77(0.3  inch) 


N  =  303  cycles. 

STEP  2.  The  second  question  can  be  answered  by  using  Equation  X-4  and  the 
applicable  stress  intensity  solution  and  elastic  magnification  factor 
for  a  through  crack  at  hole  (Equation  VII-10)  where  a  =  ac  =  1  inch  and 
r  =  0.25  inch.  In  this  case,  j{^)  =  0.75  (see  Equation  VII-10). 

Compute  Nc : 
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FIGURE  60  EXAMPLE  PROBLEM  -  CORNER  CRACK  AT  BOLT  HOLE  LOADED  IN 
TENSION 
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2 


1 


N 

c 


7Tl0"13(l-0.02)2(70,000psi)2 


_ 130,000  psiVin6h _ 

0.75  (70,000  psi)  V7r(1.0  inch) 


m  130,000  psiVinch 

-  In  - 1  * -  '  ~ 

0.75  (70,000  psi)  Vtt(1.0  inch! 

N  =  84.8  cycles.* 
c 


Thus,  for  a  sinusoidal  loading  condition  of  70  ksi  maximum  stress,  the 
0.3  inch  corner  crack  at  a  bolt  hole  will  grow  through  the  plate  in  «  300 
cycles,  and  will  approach  critical  conditions  in  385  cycles.  In  an  aircraft 
application,  this  would  mean  that  approximately  85  cycles  of  life  remain 
prior  to  fracture,  once  the  crack  completely  penetrates  the  plate. 

Equation  X-2  can  be  used  in  conjunction  with  graphical  integration  to 
obtain  load  (stress)  or  K  vs.  crack  growth  curves  for  particular  problems. 

A  good  review  of  this  procedure  is  contained  in  Reference  67,  where  predictions 
of  crack  growth  rate  using  fracture  mechanics  principles  are  outlined  and  com¬ 
pared  with  data  from  simulated  tests. 


X.7  ENVIRONMENTAL  AND  LOAD  HISTORY  EFFECTS 

If  environmental  conditions  prevail,  it  can  be  expected  that  the  values 
of  cyclic  life  will  be  reduced  on  the  order  of  20-30  percent  for  certain  mate¬ 
rials.  There  are  data  available  on  environmental  crack  growth  which  can  be 
used  to  modify  the  predicted  curves  based  on  Equations  X-2,  X-3  (see  Biblio¬ 
graphy).  It  must  be  remembered  that  for  varying  load  history,  the  order  of 
integration  is  extremely  important.  To  better  illustrate  the  effect  of  load 
history  on  crack  growth,  refer  to  Figure  61.  This  is  a  typical  representation 
of  the  delays  in  crack  growth  rate  caused  by  the  application  of  three  success¬ 
ive  cycle  overloads  during  a  sinusoidal  load  spectrum.  Data  of  this  type  may 
be  found  in  Reference  63. 

Until  the  parameters  which  influence  crack  growth  are  fully  understood, 
any  estimates  of  service  crack  growth  conditions  must  be  considered  as  pre¬ 
liminary  estimates  only. 


X.8  CRACK  GROWTH  DURING  PRESSURE  VESSEL  PROOF  TESTING 

This  special  class  of  problems  will  not  be  explored  in  this  report.  For 
those  who  are  concerned  with  this  problem,  an  excellent  summary  of  crack 
growth  analysis  for  part-through  cracks  in  pressure  vessels  is  contained  in 
Reference  42. 


*  It  can  be  seen  that  critical  conditions  are  being  rapidly  reached  (87.5 
ksi^inch  — -Kic),  and  this  value  is  not  unexpected  at  the  high  crm  . 
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FIGURE  61  SCHEMATIC  OF  DELAY  IK  CRACK  GROWTH  DUE  TO  STRESS  OVERLOAD 
IK  ALUMINUM  ALLOYS  (REF.  63) 
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XI  ENVIRONMENTAL  EFFECTS 


XI. 1  STRESS  CORROSION  CRACKING 

The  traditional  measure  of  stress  corrosion  is  based  on  the  time-to-failure 
of  a  smooth  specimen  at  different  stress  levels  in  the  environment  of  interest. 
However,  the  effect  of  environment  on  crack  initiation  (production  of  pitting) 
and  the  slow  growth  of  subsequent  cracks  cannot  be  separated.  In  design  prac¬ 
tice,  the  presence  of  crack  or  crack-like  defects  must  be  considered  to  occur 
in  the  structure.  Thus  far  wa  have  touched  lightly  on  the  effect  of  environ¬ 
ments  on  static  or  fatigue  crack  growth  and  the  role  of  stress  intensity  as  a 
correlating  parameter.*  An  example  of  the  trend  of  dependence  of  crack  growth 
on  K  is  shown  in  Figure  62  for  wedge  open  loaded  specimens  (K  decreases  with 
increasing  a)  from  the  data  of  Smith^6®).  Such,  crack  growth  is  often  referred 
to  as  sub-critical  crack  growth  or  delayed  failure  under  environmental  conditions. 

Experimental  measurement  of  sub-critical  stress  corrosion  susceptibility 
involves  the  measurement  of  time-to-failure  for  precracked  specimens  at  vari¬ 
ous  load  levels.  By  reducing  the  load  (stress)  in  steps,  a  threshold  level 
of  Kj  is  reached  below  which  no  failure  should  occur.  This  level  is  called 
KIscc*  T^e  relationship  between  this  environmentally  determined  threshold 
stress  intensity  (KIscc)  and  the  plane  strain  fracture  toughness  in  air  pro¬ 
vides  a  measure  of  material  stress  corrosion  susceptibility  which  can  be  used 
in  material  selection  and  design. 

Data  presentation  for  statically  loaded  specimens  can  take  one  or  two 
forms  as  shown  in  Figure  63:  the  applied  stress  intensity  or  normalized 
(to  Kjj.)  stress  intensity  vs.  time  to  failure.  Comparisons  of  a  "wet"  and  "dry" 
environment  can  be  made  in  this  manner. 

XI. 2  ENVIRONMENTAL  FATIGUE  CRACK  GROWTH 


The  effect  of  environment  on  fatigue  crack  growth  is  best  summarized  by 
examining  data  reported  by  Hartman^9)  on  2024-T3  clad  material.  In  these 
data,  the  effect  of  a  moist  environment  had  a  10-fold  influence  on  crack  growth 
at  low  values  of  K  over  a  "dry  argon"  reference  environment.  Trends  for 
2024-T3  shown  in  Figure  64  for  two  values  of  the  stress  ratio  are  typical  for 
7075-T6  also.  In  the  presence  of  a  saline  or  other  more  corrosive  media,  more 
drastic  changes  can  be  expected.  As  with  delayed  failure,  the  prediction  of 
environmental  fatigue  crack  growth  is  not  a  simple  task,  due  to  the  many 
parameters  which  influence  fatigue  and  compounded  by  the  presence  of  environ¬ 
mental  effects. 


*  It  bears  repeating  that  the  crack  tip  stress  intensity  will  represent  the 
driving  force  with  reasonable  accuracy  in  materials  which  exhibit  limited 
plasticity.  Actual  cases  must  be  verified  by  experimental  results. 
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FIGURE  62  TYPICAL  BEHAVIOR  OF  NET  STRESS  AND  K_  FOR  A  WEDGE  FORCE- 
LOADED  SPECIMEN  1 
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FIGURE  63  DELAYED  FAILURE  DUE  TO  ENVIRONMENTAL  CRACK  GROWTH 
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FIGURE  64  THE  EFFECT  OF  MOISTURE  ON  FATIGUE -CRACK  GROWTH  IN  A  CLAD 
2024-T3  ALUMINUM  ALLOY  (REF.  69) 


XI. 3  TEMPERATURE  AND  STRAIN  RATE  EFFECTS  ON  CRITICAL  K 


The  available  data  on  temperature  and  strain  rate  effects  indicates  in 
general  that  the  lower  strength  meterials  such  as  structural  steels  can  be 
expected  to  exhibit  a  sensitivity  to  these  parameters.  These  materials  have 
received  more  definitive  examination  due  to  their  pronounced  changes  in 
strength  and  toughness  with  temperature  and  strain  rate. 

If  it  is  suspected  in  design  that  temperature  and/or  high  strain  rates 
will  be  operative,  it  would  be  expedient  to  determine  the  critical  stress 
intensities  for  maximum  or  minimum  temperature  and  strain  rates  of  interest. 

In  some  applications,  thermal  cycling  may  occur  during  service  (e.g., 
supersonic  aircraft),  which  also  should  be  considered  in  the  evaluation  of 
fracture  criteria. 


XI. 4  SUMMARY 


To  insure  satisfactory  results  from  a  fracture  mechanics  evaluation  of 
a  specific  design,  the  critical  stress  intensity  for  the  environment  of  inter¬ 
est  should  be  determined.  Environment  enhanced  crack  growth  can  be  correlated 
by  the  use  of  stress  intensity  for  sub-critical  crack  growth  and  fatigue  crack 
propagation.  For  those  interested  in  the  application  of  fracture  mechanics 
to  stress  corrosion  cracking  and  the  current  limitations,  a  good  review  is 
presented  in  Reference  70. 

In  the  overall  view  of  design  against  fracture,  the  role  of  environment 
cannot  be  avoided  and  must  be  represented  in  the  analysis  by  modification  to 
the  critical  stress  intensity. 
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XII  FRACTURE  TESTING  PROCEDURES 


XII. I  FRACTURE  TOUGHNESS  TESTING 


The  objective  of  the  fracture  toughness  test  is  to  determine  an  accurate 
value  of  the  critical  stress  intensity  for  the  material  in  a  given  environment. 
Most  of  the  historic  effort  has  been  expended  in  determination  of  the  critical 
plane  strain  fracture  toughness  (Kjc).  The  following  general  comments  apply 
to  all  specific  types. 

.  Specimen  must  contain  a  sharp  fatigue  induced  crack. 

.  Specimen  must  exceed  minimum  geometric  dimensions  (see  e.g.  Section  VI). 

.  Crack  extension  and  loading  must  be  accurately  determined. 

The  details  which  are  required  to  obtain  Kjc  data  will  not  be  covered  in 
this  report  as  documents  are  available  which  cover  these  items  in  great  detail. 
(See  Refs.  7,  8,  and  current  ASTM  proposed  test  method,  Ref.  71). 

At  present,  there  is  no  standard  test  method  for  fracture  toughness  test¬ 
ing  of  thin  gage  materials  which  do  not  exhibit  brittle  like  (plane  strain) 
behavior,  A  special  Task  Force  within  Committee  E-24  of  the  American  Society 
for  Testing  and  Materials  has  been  established  to  study  and  eventually  recom¬ 
mend  standards  for  thin  section  fracture  testing. 

XII. 2  PLANE  STRAIN  TOUGHNESS  TESTING 


Starting  from  the  center  cracked  through  the  three  point  bend  and  com¬ 
pact  tension,  the  evolution  of  plane  strain  fracture  toughness  specimens  (Kjc) 
is  shown  in  Figure  65,  From  this  progression,  two  specimens,  the  bend  and 
compact  tension,  have  been  proposed  by  ASTM  (Ref.  71)  as  they  have  been  shown 
to  fulfill  the  requirements  for  valid  and  repeatable  data. 

Accurate  measurements  are  made  of  both  load  and  displacement  through  the 
use  of  autographic  recording  of  test  machine  load  and  the  displacement  of  two 
precisely  located  points  above  and  below  the  crack.  With  the  aid  of  suitable 
calibration  curves  of  specimen  displacement  vs.  elastic  compliance  established 
with  the  displacement  gage  in  place,  accurate  measurements  of  crack  extension 
are  possible.  A  typical  calibration  curve  is  shown  in  Figure  66  in  terms  of 
relative  crack  extension  for  a  center  crack  geometry.  The  specimen  compliance 
(vEBj  determined  from  the  fracture  test  where  v  is  the  crack  displacement, 
P  is  the  load  associated  with  displacement,  and  E  and  B  are  Modulus  and  speci¬ 
men  thickness,  respectively.  Any  particular  calibration  applies  to  any  other 
geometrically  similar  specimen  with  different  Modulus  or  thickness.  Thus,  the 
information  obtained  from  the  autographic  recording  during  the  fracture  test 
(see,  for  example,  Figure  67)  can  be  translated  into  stress/crack  length  and 
critical  points  established.  With  these  critical  points,  the  critical  plane 
strain  fracture  toughness  (Kic)  can  be  determined  from  the  suitable  stress 
intensity  equation  for  the  geometry  and  loading  configuration.  Many  different 
types  of  load  displacement  curves  are  possible;  for  a  further  discussion  see 
Reference  8. 
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FIGURE  66  TYPICAL  CALIBRATION  CURVE  FOR  A  CENTER  CRACK  GEOMETRY 
WITH  0.2  INCH  GAGE  POINT  DISTANCE 
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FIGURE  67  REPRESENTATION  OF  TYPICAL  DISPLACEMENT  GAGE-LOAD  RECORD  IN 
FRACTURE  TOUGHNESS  TEST 
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To  obtain  a  valid  value  of  plane  strain  fracture  toughness,  it  is  re¬ 
quired  that  both  the  specimen  thickness,  B,  and  crack  length,  a,  exceed  2.5 
times  the  square  of  the  ratio  of  plane  strain  fracture  toughness  to  0.2  percent 
offset  yield  strength  (see  equation  VI-8).  As  shown  in  Section  VI,  this  re¬ 
quirement  sets  plastic  zone  limits  so  that  elastic  fracture  mechanics  principles 
apply  directly.  One  advantage  of  the  displacement  measurement  of  crack  exten¬ 
sion  is  that  the  value  of  crack  length  obtained  includes  both  crack  length  and 
the  increment  due  to  plasticity.  Therefore,  approximate  expressions  for  plas¬ 
ticity  corrections  need  not  be  used. 

With  knowledge  of  the  critical  crack  length  and  associated  critical  load, 
the  critical  stress  intensity  factor  can  be  calculated  for  the  crack  geometry. 
This  stress  intensity  is  designated  Kq  until  such  time  as  it  is  shown  that  all 
criteria  of  Reference  71  are  met.  Then,  and  only  then,  can  this  Kq  value  be 
designated  a  valid  plane  strain  fracture  toughness  value  (Kjc).  Reference  71 
describes  the  limitations,  specimen  design,  and  testing  criteria  required  in 
plane  strain  fracture  toughness  testing  and  should  be  consulted  for  additional 
information. 

XII, 3  SOURCES  OF  PLANE  STRAIN  AND  PLANE  STRESS  DATA 
XII, 3(a)  Plane  Strain  (Kjr) 

In  evaluating  reported  plane  strain  fracture  toughness  data,  it  is  imper¬ 
ative  that  a  prime  criterion,  the  relationship  of  Equation  VI-8,  be  met.  As  a 
general  rule,  any  data  which  has  been  taken  using  the  proposed  method  (Ref.  71) 
can  be  considered  as  valid  data  as  long  as  the  limited  plasticity  of  Equation 
VI-8  is  fulfilled.  Those  organizations  which  are  members  of  ASTM  Committee 
E-24  or  who  follow  the  criteria  of  the  proposed  testing  method  can  be  used  as 
sources  for  valid  Kjc  data.  Older  data  which  does  not  meet  this  criteria  must 
be  considered  as  Kq  values  and  used  only  for  the  thickness  and  specimen  geometry 
for  the  reported  data. 

XII. 3(b)  Plane  Stress  (Kq) 

In  general,  data  on  Kc  must  come  from  specific  tests  employing  the  mate¬ 
rial,  thickness,  and  specimen  width  for  the  intended  application.  The  net 
section  stress,  at  fracture,  must  not  exceed  0.8  times  the  yield  stress 
(  crnet  <  0.8  crys)  and  crack  buckling  must  be  restrained  or  accounted  for  in 
the  analysis.  The  preceeding  are  general  criteria  and  must  be  considered  as 
such  until  specific  testing  methods  are  available  from  ASTM  Committee  E-24. 

It  is  recommended  that  those  who  are  interested  in  further  discussion  of 
the  fracture  testing  of  metals  consult  References  7,  8,  and  71. 

XII, 4  CRACK  GROWTH  RESISTANCE  (R  CURVE  CONCEPT) 


In  Section  II  (Figure  1)  the  crack  growth  resistance  concept  of  Irwin  is 
depicted.  Further  discussions  of  this  concept  are  contained  in  Reference  72 
and  73.  Krafft^22),  et  al,  hypothesized  that,  neglecting  environmental  effects 
and  with  controlled  testing,  the  crack  growth  resistance  (resistance  to  crack 
extension)  is  primarily  a  function  of  absolute  crack  extension  (ac  “  &o)  an< ^ 
independent  of  initial  crack  length.  This  concept  leads  to  many  implications, 
such  as  uniqueness  of  R  curves  for  materials,  etc. 
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Recent  research  has  shown  that  crack  growth  resistance  offers  a  means  of 
determining  a  critical  plane  stress  fracture  toughness,  Kc  for  (74)  those 
thin  materials  whose  fracture  properties  are  highly  geometric  (width  and 
original  crack  length)  dependent.  Thus,  it  is  hoped  that  the  crack  growth 
resistance  curve  for  a  given  material  and  thickness  (again  neglecting  loading 
and  environmental  conditions)  can  be  used  to  determine  critical  crack  driving 
forces  or  stress  intensities  for  specimens  of  other  geometries.  The  results 
of  Heyer  and  McCabe'^'  are  encouraging  in  this  respect  and  bear  further  study 
as  additional  crack  growth  resistance  data  for  various  crack  geometries  becomes 
available . 
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